Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,866,488 G→C 5.4% intergenic (+23/+71) ygbN → / ← rpoS putative transporter/RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,866,4880GC 5.4% 258.3 / 7.0 92intergenic (+23/+71)ygbN/rpoSputative transporter/RNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base G (45/40);  new base C (2/3);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

TGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTG  >  NC_000913/2866343‑2866636
                                                                                                                                                 |                                                                                                                                                    
tGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCAt                                                                                                                                                   >  1:689146/1‑149 (MQ=255)
tGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCAt                                                                                                                                                   >  1:142179/1‑149 (MQ=255)
tGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCAt                                                                                                                                                   >  1:1054081/1‑149 (MQ=255)
 cgATGGTGACCAAATATCTGGGGGTGTCGGTATCAGACGGGCTGAAAACCTGGACAGTGTTAACGACCATTCTCCGGTGTACCGGATGCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATc                                                                                                                                                  <  2:485579/148‑1 (MQ=255)
  gATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGCCCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCg                                                                                                                                                 >  1:1186776/1‑149 (MQ=255)
   tttGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCgg                                                                                                                                                <  1:983439/148‑1 (MQ=255)
     tgtgACCAAAAATCGGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTCTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGAc                                                                                                                                              <  1:608226/149‑1 (MQ=255)
     ggtgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTAACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGAc                                                                                                                                              <  2:206484/148‑1 (MQ=255)
       tgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCt                                                                                                                                            >  1:462369/1‑149 (MQ=255)
        gACCAAAGATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTATTTACCAGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAGCtt                                                                                                                                           >  2:710082/1‑149 (MQ=255)
         ccccaaTATGTGGGGGTGTCGGTAGCAGTCGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTCCCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCttt                                                                                                                                          <  1:1118923/145‑1 (MQ=255)
         aCCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTGAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGACCttt                                                                                                                                          <  2:1024534/149‑1 (MQ=255)
          ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtttt                                                                                                                                         >  1:329366/1‑149 (MQ=255)
           cAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTa                                                                                                                                        >  1:845924/1‑149 (MQ=255)
             aaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTAtt                                                                                                                                      >  1:842974/1‑149 (MQ=255)
             aaTATCTGGGGTGGTCGGTAGCAGACGGTCTGCAAACCTGGACAGTGTTAACGACCATTCTCGGTATTACCGGATTCTTAATTCCCTGGGGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTAtt                                                                                                                                      <  1:569492/149‑1 (MQ=255)
              atatGTGGGGTTGTCGGTAGCAGACGGGCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTTATAAAAAGGTCCGATGGGCATCGGACCTTTTATtg                                                                                                                                     <  2:352797/149‑1 (MQ=255)
               tctCTGGGGTTGTCGGTACCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATtgt                                                                                                                                    <  1:1017792/147‑1 (MQ=255)
               tatCTGGGGGTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATtgt                                                                                                                                    <  1:998691/149‑1 (MQ=255)
                atCTGGGGTTGGCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATtgtg                                                                                                                                   <  1:908722/149‑1 (MQ=255)
                atCTGGGGTCGCCGGTACCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGCTAATTTGACCTTAATAACAAGGGCCGATGGGCAGCGGACCTTTTATtgtg                                                                                                                                   <  1:1078524/149‑1 (MQ=255)
                 tCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACATTTTATTGTGc                                                                                                                                  >  1:1160786/1‑149 (MQ=255)
                  cTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACAGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGca                                                                                                                                 <  2:904334/149‑1 (MQ=255)
                          gTCGGTAGCAGACGGTCTGAAACCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTCCCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAg                                                                                                                         <  1:229727/149‑1 (MQ=255)
                          gTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAg                                                                                                                         >  1:496580/1‑149 (MQ=255)
                           tCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAgg                                                                                                                        >  1:912830/1‑149 (MQ=255)
                              gTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCa                                                                                                                     <  1:529805/149‑1 (MQ=255)
                                aGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGc                                                                                                                   <  2:525767/149‑1 (MQ=255)
                                aGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGc                                                                                                                   <  1:1047034/149‑1 (MQ=255)
                                 gCAGACGGGCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAgtc                                                                                                                  <  1:1116136/149‑3 (MQ=255)
                                   aGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCTGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCta                                                                                                                >  2:827189/1‑148 (MQ=255)
                                   aGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGCATTCTTAATTACAAGGTGCGTATGTGCGGAAATTTGAACATAATAAAAAGGTCCCATGGGCATCTGACCTATTATTTTGCACAGAAAAGTACACCCTc                                                                                                                >  2:663999/1‑149 (MQ=255)
                                    gACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCg                                                                                                               >  2:913609/1‑149 (MQ=255)
                                       ggTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCtt                                                                                                            <  1:914707/149‑1 (MQ=255)
                                       ggTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCtt                                                                                                            >  1:333409/1‑149 (MQ=255)
                                       ggTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGAGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCtt                                                                                                            >  2:841960/1‑149 (MQ=255)
                                        ttCTGAATAAGTGGACATTGTTAACTACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTAAGGGCGGTAATTTGACCTTAATAAAAAGGGCCGATGGTCATCGGACCTTTTATTGTGCACAGTAAAGGACAGCCTCGCTTg                                                                                                           <  2:341857/148‑1 (MQ=255)
                                         tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                          >  1:681241/1‑149 (MQ=255)
                                         tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                          >  1:394267/1‑149 (MQ=255)
                                         tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                          >  1:1168280/1‑149 (MQ=255)
                                           tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgaga                                                                                                        >  1:1122348/1‑149 (MQ=255)
                                            gAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGAc                                                                                                       >  2:502204/1‑149 (MQ=255)
                                            gAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGAc                                                                                                       >  1:405424/1‑149 (MQ=255)
                                             aaaaCCTGGACAGTGGTAACGACCATTCTCGGGTTTCCCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACt                                                                                                      <  2:656725/149‑1 (MQ=255)
                                             aaaaCCTGGAAAGGGAGAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGAGCGTATGGGCGGTAATTTGACCTAAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACt                                                                                                      <  2:653908/149‑1 (MQ=255)
                                              aaaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGTATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTg                                                                                                     >  1:1053311/1‑149 (MQ=255)
                                                aCCTGGACAGTGTTAACGACCTTTCGCGGTTTTCCCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGATACTGGc                                                                                                   <  1:646244/149‑1 (MQ=255)
                                                   tGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGTCTCAAGCGAGGCTGGCCtt                                                                                                <  1:1023790/149‑1 (MQ=255)
                                                    gggcAGGGTTAACGACCATTCTCGGCTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCttt                                                                                               <  2:980500/146‑1 (MQ=255)
                                                    ggACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGGATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCttt                                                                                               <  2:539503/149‑1 (MQ=255)
                                                     gACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                              <  2:878643/149‑1 (MQ=255)
                                                     gACAGGGTTAACGACAAATCTCGGTTTTACCGGATGGTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCATAAAAAAAAGGTCCGAAGGTCATCGTAAATTTTATTGTGCACAGAAAAGGACAGCCTCTCTTGAGACTGGCCTTTc                                                                                              <  2:979935/149‑1 (MQ=255)
                                                      aCAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCt                                                                                             <  2:1160786/149‑1 (MQ=255)
                                                          tgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCTGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGATTGAGACTGGCCTTTCTGACa                                                                                         >  1:517282/1‑149 (MQ=255)
                                                          tgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACa                                                                                         >  1:75378/1‑149 (MQ=255)
                                                             tAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGAt                                                                                      >  1:96191/1‑149 (MQ=255)
                                                             tAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAACTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGAt                                                                                      <  1:171407/149‑1 (MQ=255)
                                                             tAACGACCATTATCTGTTTTACCGGATTATTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTAGCTTGAGACTGGCCTTTCTGACAGAt                                                                                      <  2:776167/149‑1 (MQ=255)
                                                              aaCGCCCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATg                                                                                     >  1:581023/1‑149 (MQ=255)
                                                               aCGACCATTCTCGGTTTTACCGGATTCTAAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGc                                                                                    <  1:871355/149‑1 (MQ=255)
                                                                 gACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGctt                                                                                  >  2:691593/1‑149 (MQ=255)
                                                                 gACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGctt                                                                                  >  1:325746/1‑149 (MQ=255)
                                                                     aTTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttactt                                                                              >  2:168229/1‑149 (MQ=255)
                                                                       tctcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCAGCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCttacttac                                                                            <  2:937380/149‑1 (MQ=255)
                                                                        ctcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGTGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCCGACCTTTTATTGTGCACAGAAAAGGCCATCCTCGCTTGAGCATGGCCGTTCTGACAGATGCTtacttact                                                                           >  1:388109/1‑149 (MQ=255)
                                                                             ttttACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCgg                                                                      <  1:853480/149‑1 (MQ=255)
                                                                                   cGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAgc                                                                >  1:184839/1‑149 (MQ=255)
                                                                                       tgctaaaTTACCTGGAGCGTATGGGCTGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTACTGACAGATGCTTACTTACTCGCGGAACAGCGCtt                                                            <  2:500272/144‑1 (MQ=255)
                                                                                        tCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTc                                                           <  2:1168280/149‑1 (MQ=255)
                                                                                               tACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTc                                                    <  1:28089/149‑1 (MQ=255)
                                                                                                  cTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGc                                                 <  1:841960/149‑1 (MQ=255)
                                                                                                   tGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGcc                                                >  2:47183/1‑149 (MQ=255)
                                                                                                          tATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCTCTTCGATATTCAGCCCCTGCGt                                         >  2:348454/1‑149 (MQ=255)
                                                                                                          tATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGt                                         >  2:965993/1‑149 (MQ=255)
                                                                                                            tGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGttt                                       >  2:945870/1‑149 (MQ=255)
                                                                                                               gCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCGTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                    <  2:689146/149‑1 (MQ=255)
                                                                                                               gCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACATATGCTTACTTACTCGCGGAACAGCGCGTCGATATTCAGCCCCTGCGTTTGCa                                    >  2:1100754/1‑149 (MQ=255)
                                                                                                               gCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                    <  2:1171767/149‑1 (MQ=255)
                                                                                                                cGGTAATTTGACCTTAATAAAAATGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAg                                   <  2:1053311/149‑1 (MQ=255)
                                                                                                                 ggTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAgg                                  >  1:1034950/1‑149 (MQ=255)
                                                                                                                       ttGACCTTAATAAAAAGGTCCGATGTGCATCGGACCTTTTATTGTGCACAGAAAAGTCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTcg                            <  1:817069/149‑1 (MQ=255)
                                                                                                                           ccTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCa                        <  1:913609/149‑1 (MQ=255)
                                                                                                                            cTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGACTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaa                       >  1:707953/1‑149 (MQ=255)
                                                                                                                              taataaAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAAc                     <  2:329366/149‑1 (MQ=255)
                                                                                                                              acaataaaaGGTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTAAGCCCCTGCGTTTGCAGGATTTCGCGCAAAc                     >  1:209700/6‑149 (MQ=255)
                                                                                                                                  aaaaaGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGgcg                 >  2:632520/1‑149 (MQ=255)
                                                                                                                                    aaaGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa               <  2:142179/149‑1 (MQ=255)
                                                                                                                                      aGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGGAGGATTTCGCGCAAACGGCGCAgg             >  1:448603/1‑149 (MQ=255)
                                                                                                                                       ggTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc            >  2:552554/1‑149 (MQ=255)
                                                                                                                                        gTCCGATGGGCATCGGACCTTTTATTGTTCACAGAAAAGGCCAGCCTCGCTTGAGACTGTCCTTTCTGACAGATGCTTACTTACTAGCGGAACAGCGCATCGATATTCAGCCCCTGCGTTTGCAGGATTTCGAGCAAACGGCGCGGGcc           >  2:476824/1‑149 (MQ=255)
                                                                                                                                           cGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGGACAGCGCGTCGATATTCAGCCCCTGCGGTTGCAGGATTTCGCGCAACCGGCGCAGGCGTTc        >  2:94712/1‑149 (MQ=255)
                                                                                                                                           cGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc        >  1:756048/1‑149 (MQ=255)
                                                                                                                                                 gCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg  >  2:129760/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTG  >  NC_000913/2866343‑2866636

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: