Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,542,072 G→C 12.3% intergenic (+52/+51) nadE → / ← PP_4870 NH(3)‑dependent NAD(+) synthetase/azurin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,542,0720GC12.3% 463.8 / 39.3 191intergenic (+52/+51)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base G (70/94);  new base C (15/10);  total (85/106)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

GCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCG  >  NC_002947/5541927‑5542219
                                                                                                                                                  |                                                                                                                                                   
      gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGAACCGCGAATAAGATGGTGAACCCACCGAAGCATTagc                                                                                                                                           >  2:4463301/1‑150 (MQ=255)
      gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGccatagc                                                                                                                                           >  1:6552026/1‑144 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                         <  1:7097956/150‑1 (MQ=25)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                         <  1:932690/150‑1 (MQ=25)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                         <  2:148373/150‑1 (MQ=25)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                         <  2:8338593/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                         <  2:7983939/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                         <  1:4420434/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                         <  2:6531424/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                         <  2:5403176/150‑1 (MQ=255)
         gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCggg                                                                                                                                        <  1:4477474/150‑1 (MQ=255)
          ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGt                                                                                                                                       <  1:5117890/150‑1 (MQ=25)
          ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGt                                                                                                                                       <  2:398643/150‑1 (MQ=25)
            ggTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGa                                                                                                                                     <  1:2374991/150‑1 (MQ=25)
                gACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCGCCCCCCAAAACGAGGGGGTACCCACCCCCGCAGTAGCGGGTGAAccc                                                                                                                                 >  2:2190563/1‑150 (MQ=255)
                    cggCAGCGTTCGACAACATCTTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTc                                                                                                                             <  1:7944844/149‑1 (MQ=255)
                    aGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCCGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTc                                                                                                                             <  2:7777743/150‑1 (MQ=21)
                     ggCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAAAACGAGGGTGAAACCACCGACGCAATAGCGGGTGGAACCGCGcc                                                                                                                            >  2:6958904/1‑150 (MQ=255)
                            ttCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTAACCCGCGGATACGGGGGTGCACCCACCGCCGGCTTCGGGGGTGAAGCCGCTCCTGCAggg                                                                                                                     >  2:7092570/1‑150 (MQ=255)
                            tccGAAATCATCGAGGCCACCTACCGCACGCCCCACCACACCGGGGAACTGCGGTTCGCTCCTTGATATGTGGCCTGAAGGAGGGGTTTAACCCTGTATTACTTGGGTCAACCCACAGACGCTTTCGCGGGTGAACCCGCTCCTGCattt                                                                                                                     <  2:6736596/148‑4 (MQ=255)
                                    catcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGCATACGAGGGTGAACCCACCGACGCAGTCGCGGGTGAACCCGCTCCGGCAGGGGTCGCCgg                                                                                                             >  2:4314900/1‑149 (MQ=255)
                                      tggtcgCCACCCACCGCACGACCCCCCACACGCGCGACCTGCCGTGCGCTCCTTTATACGGGGCCGGCAGGAGCGGGTTCACCCGCGAATCCGCGGGTGACCCCACCGACGCCTTCGGGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCa                                                                                                           <  2:140031/148‑1 (MQ=255)
                                      tcgtccCCACCCCCCCCACGACCCACCCCACGCGCGAACTGCGGTTCGCCCCTTTATACGTGGCCGGCCGGTGGGGTTTCACCGGCTAATACGCGGGTTAACCCCCCGACGCCTTCGCGGTTGAACCCGCTCCTGCAGGGTTCGCCGTCa                                                                                                           <  2:3915323/150‑1 (MQ=255)
                                        gtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCAGCGAAGACGAGGGTGAACCCACCGACGCAATCGCGGGTGAACCCCCTCCTGCAGGGTTCGCCGGCAGc                                                                                                         >  2:92140/1‑150 (MQ=255)
                                            ccACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGCCGCCTTCGCGGGTGACCCCGCGCCTGCCGGGTTCGCCGTCAGCCttt                                                                                                     >  2:5862935/1‑150 (MQ=255)
                                             cACCTACCGCAAGACCCCGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTg                                                                                                    <  2:6387194/150‑1 (MQ=255)
                                                  aCCGCAAGACCCAGCACAAGCGCGAACTGCGGTTCGCTCCTTGATACGTGGCCTGCCGGTGGGGGTTCACCCGCTAATCCGCGGGTCAACCCACCGACGCATTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTAc                                                                                               <  1:2800553/150‑1 (MQ=255)
                                                    cccaaccccccGCACACGCGCGACCTGCCGTCCGCTCCTGGATACGGGGCCTGCAGGCGCGGTTCCCCCCGCTAATACGAGGGTGACCCCACCGACGCATCCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                                             <  2:2983102/143‑1 (MQ=255)
                                                    cGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCAGCGCCAGCAGGGTTCGCCGGCAGCCGATGATTaagt                                                                                             >  2:4963554/1‑147 (MQ=38)
                                                    cGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCCCGAATACGAGGGTGAACCCACCCACGCAGTCGCGGGGGACCCCCCGCCGGCCGGGTTCGCCGGCACCCGTTGATTAcat                                                                                             >  2:140616/1‑148 (MQ=255)
                                                      cAACACCCACCACACGCGCTAACTGCCTTTCGCTCCTTTATAAGGGGCCTGCAGGAGCGGTTTCCCCCGCTAATACGCGGGTGAACCCCCCGACGCTTGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGa                                                                                           <  1:5683578/150‑1 (MQ=255)
                                                          aCCCCGCACAAGCGCGCCCGGCGGTCCGCACCTTGATACGGGGCCTGCAGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                                       <  2:7522606/150‑1 (MQ=25)
                                                              aGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGAATACGAGGGTGAACCCACCGACGCAATCGCGGGTGAACCCGCGCCGGCAGGGTTCGCCGGCCGCCGTTGATTACTTGACGATGAcc                                                                                   >  2:4334550/1‑150 (MQ=255)
                                                                cacaAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCtt                                                                                 <  2:7893357/150‑1 (MQ=255)
                                                                  cacgcgcCAACTGCCGTCCGCTCGTTTATATGGGGCCTGCTGGCGCGGGTTCACCCGGTAAACCGTGGGTGAACCCACCGCCGCTTTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGc                                                                               <  2:3170115/147‑1 (MQ=255)
                                                                       gcgACCTGCCGTTCGCTCGTTGATCCGGGGCCTGCTGGAGGGGGTGCCCCCGCGAATAGGGGGGTGCCCCCACGGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                                          <  2:2275259/150‑1 (MQ=255)
                                                                        cgAACTGCGGTTCGCTCCTTTATACGTGCCGTGCTGGCGCGGGTTCACCCGCTAATAGGGGGGTGAACCCACCGACGCTGTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                         <  2:3671456/150‑1 (MQ=255)
                                                                        cgAACTGCCGTTCGCTCCTTTATATGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                         <  2:4866169/150‑1 (MQ=255)
                                                                        cgAACTGCCGTTCGCTCCTTTATAAGTGGCCTGCAGGAGCGGGTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                         <  2:1330933/150‑1 (MQ=255)
                                                                        cgAACTGCCGTTCGCTCCTTGTTAAGTGGCCTGCTGGCGCGGGTTCCCCCGCGAATACGAGGGTGAACCCACCTACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                         <  2:1565416/150‑1 (MQ=39)
                                                                        cgAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                         <  2:8692609/150‑1 (MQ=255)
                                                                          accggCCGTTCGCTCCTTGACTAGTGGCCTGCAGGAGCGGGTTCACCCGCGATTAGGTGGGTGAACCCACGGACGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtcatc                                                                       <  2:5013263/146‑1 (MQ=38)
                                                                               ccGTTCGCTCCTTGATATGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGTGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCga                                                                  <  2:1226772/150‑1 (MQ=255)
                                                                                cGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgag                                                                 <  2:8606895/150‑1 (MQ=255)
                                                                                       tCCTTGATAAGTGGCCTGCAGGAGCGGGTTAAACCCCCAAAAACAGGGTGAAACCAACCACGCCATAGCGGGTGAACCCCCGCCTGCAGGGTTCGCCGTCAGCCGGTGATTACGTGTCGAGGACCGTGCCGTTCCTCCACGCGATGGGGc                                                          >  2:3848030/1‑150 (MQ=255)
                                                                                            gATAAGTGGCCTGCAGGAGCGGGTTCAACCCCCCAAAACAGGGTGAACCCCCCCACGCATTAGCGGGTGAAACCCCGCCTGCAGGGGTAGCCCTCAGCCGTTTATTAAGTGACCATGACCTTGGCGTTAAGCATCGAGATGGGGGCCggg                                                     >  2:4515695/1‑150 (MQ=255)
                                                                                            gATAAGTGGCCTGCAGGAGCGGGTTCACCCCCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAAACCGCTCCTGCAGGGTTCGCCGTCAGCCTTTTATTACTTGACGATGACCATGCCTTTCATCATCGAGATGTGGCCCggg                                                     >  2:5668116/1‑150 (MQ=255)
                                                                                              aacgTGGCCTGCTGGCGCGGTTTCACCCGCCAATACGCGGGTCAACCCACCCACGCCTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGaa                                                   <  2:2194906/147‑1 (MQ=255)
                                                                                                  tGGCCTGCAGGAGGGGGTCCACCCGCGAATACGCGGGTGAACCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTg                                               <  2:3279590/150‑1 (MQ=255)
                                                                                                  gggCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTg                                               <  2:4607262/149‑1 (MQ=255)
                                                                                                   ggccggccGGCGGGGGTCCCCCCGCTAATCGGCGGGTCACCCCCCCGCCGCATGCGGGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGc                                              <  1:449048/150‑1 (MQ=255)
                                                                                                   ggCCTGCCGGTGCGGTTCCCCCCGCTAATCCGGGGGGCACCCCACCGCCGCATTCGGGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCACCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGc                                              <  1:6727642/150‑1 (MQ=255)
                                                                                                        gCAGGAGCGGGTTGAACCCCGCAAACGAGGGTGTAACCACCCACCCCGTCGCGGGTGAACCCCCTCCTGCCGGGGTCGCCCTCAGCCATTTTTTACTTTACGAAGAACGAGGCATTCATCAACGAGAAGTGGCCCGGGAAGGTGCagaag                                         >  2:3367124/1‑150 (MQ=255)
                                                                                                             cgggggTCCCCCCGCAAAACCGCGGGTCACCCCCCCGCCGCATGCGCGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCt                                    <  1:2941298/147‑1 (MQ=255)
                                                                                                              gCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg                                   >  2:4406797/1‑150 (MQ=255)
                                                                                                               cGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGt                                  >  2:5931037/1‑150 (MQ=21)
                                                                                                                gggTTCACCCGCGAATACGTGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                                 <  2:1052574/150‑1 (MQ=255)
                                                                                                                gggTTCACCCGCGAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                                 <  2:3643978/150‑1 (MQ=255)
                                                                                                                gggTTCACCCGATAAAACGTGGGCCAACCCACCGACGCATTGGCGGGTCAACCCGCTCCTGCAGGGTTCGCCGACAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                                 <  1:2517961/150‑1 (MQ=255)
                                                                                                                 ggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAg                                <  2:8817010/150‑1 (MQ=255)
                                                                                                                  gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                               >  1:884616/1‑150 (MQ=255)
                                                                                                                  gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                               >  1:481584/1‑150 (MQ=255)
                                                                                                                   ttCACCCGCTATTACGTGGTTGAACCCACCGACGCTTTCGGGGGTGAACCCGCTCCTGCGGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc                              <  1:1307747/150‑1 (MQ=255)
                                                                                                                    tacccccGCCAATACGCGGGTCACCCCACCGCCGATTTCGGGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc                              <  1:464642/146‑1 (MQ=255)
                                                                                                                    tCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTca                             >  1:3643595/1‑150 (MQ=255)
                                                                                                                    tCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTca                             >  1:3863070/1‑150 (MQ=255)
                                                                                                                    ctcccccGCGAATACGCGGGTCACCCCACCGCCGCCTCCGGGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc                              <  1:100768/146‑1 (MQ=255)
                                                                                                                     cACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcac                            <  2:5551238/150‑1 (MQ=255)
                                                                                                                         ccGCGAATGCGTCGGTGGATTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                         >  1:7362138/1‑150 (MQ=37)
                                                                                                                         ccGCGAATGCGTCGGTGAGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTTGCCGTCAGCCGTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                         >  1:7260398/1‑150 (MQ=25)
                                                                                                                         ccGCGAATACGCGGGTGAACCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                         <  2:7403689/150‑1 (MQ=255)
                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                         >  1:7959346/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                         >  1:2033453/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                         >  2:6131921/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                         >  2:3587516/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                         >  1:8537599/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                         >  1:8800361/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                         >  1:3989995/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                         >  1:3313097/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                         >  2:1705786/1‑150 (MQ=255)
                                                                                                                          cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGc                                              <  2:5930010/128‑1 (MQ=255)
                                                                                                                          cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                        >  2:8519528/1‑150 (MQ=255)
                                                                                                                          cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                        >  1:1480735/1‑150 (MQ=255)
                                                                                                                          cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                        >  1:3615901/1‑150 (MQ=255)
                                                                                                                          cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                        >  1:3615859/1‑150 (MQ=255)
                                                                                                                           gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCg                       >  2:8631683/1‑150 (MQ=255)
                                                                                                                                tACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTc                  >  2:2461375/1‑150 (MQ=255)
                                                                                                                                 ccgcgggTCAACCCCCCGACGCAGGCGGGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTcc                 <  1:1062807/146‑1 (MQ=255)
                                                                                                                                      ggTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAGGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCtt            >  1:8803263/1‑150 (MQ=255)
                                                                                                                                       gttaaCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGAGTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg           <  2:6281814/147‑1 (MQ=255)
                                                                                                                                       gTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg           <  2:1783959/150‑1 (MQ=255)
                                                                                                                                       gTGAACCCACCGTCGCTTTCGCGGGTGCCCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg           <  1:8398156/150‑1 (MQ=255)
                                                                                                                                       gTGAACCCACCGACGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg           <  1:1314463/150‑1 (MQ=255)
                                                                                                                                       gTGAACCCACCCCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg           <  1:4951029/150‑1 (MQ=255)
                                                                                                                                          aaCCCACCTCCGATTTCGGGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa        <  1:6176517/150‑1 (MQ=255)
                                                                                                                                          aaCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa        >  1:264623/1‑150 (MQ=255)
                                                                                                                                             ccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTTGAAgtg     >  1:6270726/1‑150 (MQ=255)
                                                                                                                                              cACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtgt    >  1:717951/1‑150 (MQ=255)
                                                                                                                                               aCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGCTTACTTGACGATGACCTTGCCTTTCATCCTCGAGCTGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTc   >  2:6187574/1‑150 (MQ=255)
                                                                                                                                               aCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTc   >  2:7271764/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCG  >  NC_002947/5541927‑5542219

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: