Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,542,067 C→T 14.2% intergenic (+47/+56) 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,0670CT14.2% 382.9 / 48.2 148intergenic (+47/+56)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base C (59/68);  new base T (7/14);  total (66/82)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

Only 100 of 157 total aligned reads displayed.

CTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTG  >  NC_002947/5541921‑5542216
                                                                                                                                                  |                                                                                                                                                     
cTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCCCGACAACCAGGGTGAACCCAc                                                                                                                                                    >  2:4619169/1‑150 (MQ=255)
 ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGTGGGTTCACCCTCGTATTCGCGGGTGAAcccgct                                                                                                                                                   <  1:6985110/150‑4 (MQ=255)
 ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATAAGAGGGTGAAcccccc                                                                                                                                                   >  2:1476290/1‑150 (MQ=255)
 ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGGCCCGCGCCTAGGTGGGGGGACCCAcc                                                                                                                                                   >  1:4162532/1‑150 (MQ=255)
  tCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGACCCGCCAAAAACGGGGTGGACCCCCCg                                                                                                                                                  >  2:2566368/1‑150 (MQ=255)
    cTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGAc                                                                                                                                                <  2:4018580/150‑1 (MQ=255)
      gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCCACGc                                                                                                                                              >  2:1374341/1‑150 (MQ=255)
          gggCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTc                                                                                                                                          <  2:7021840/150‑1 (MQ=255)
          gggCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACCCaatg                                                                                                                                          >  2:1571746/1‑147 (MQ=255)
              cgCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                      <  1:795059/149‑1 (MQ=255)
              aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                      <  1:7301201/150‑1 (MQ=25)
              aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                      <  2:5334277/150‑1 (MQ=25)
              aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                      <  2:522313/150‑1 (MQ=25)
              aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                      <  2:2295221/150‑1 (MQ=25)
              aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                      <  2:542946/150‑1 (MQ=25)
              aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                      <  2:2294925/150‑1 (MQ=25)
              aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                      <  2:655004/150‑1 (MQ=255)
              aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                      <  1:4956474/150‑1 (MQ=255)
              aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                      <  2:1689240/150‑1 (MQ=255)
              aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATACGCgg                                                                                                                                      <  2:1662057/150‑1 (MQ=255)
               gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCggg                                                                                                                                     <  2:6491413/150‑1 (MQ=255)
                ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGt                                                                                                                                    <  2:7233694/150‑1 (MQ=255)
                   gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGACCTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGACTGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGaa                                                                                                                                 <  1:4052597/150‑1 (MQ=21)
                   gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                      >  1:4270111/1‑145 (MQ=255)
                   gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                      <  2:4270111/145‑1 (MQ=255)
                   gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGaa                                                                                                                                 <  1:4648403/150‑1 (MQ=255)
                         cAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCCAATAAGAGGGTGGACCCACCCACGCAATAGCGGGTGAAacccca                                                                                                                           >  1:5968311/1‑147 (MQ=255)
                          aGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCCAATACGAGGGTGAAACCACCCACCCCATCGCGGGGGAccccccac                                                                                                                          >  2:5288232/1‑146 (MQ=255)
                          aGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCGAAAACGAGGGTGAAACCAACCACGCAATCGCGGGTGAAcccccgc                                                                                                                          >  2:184499/1‑146 (MQ=255)
                           ggCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACCAGGGTGAACCCACCGACGCATTCGCGGGTGAAACCGCTcc                                                                                                                         >  2:2374291/1‑150 (MQ=255)
                            gCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCt                                                                                                                        <  2:6059133/150‑1 (MQ=255)
                                 gTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGCAAACGAGGGTGAAACCACCCACGCAATAGCGGGTGAACCCGCACCTGCacg                                                                                                                   >  2:2761148/1‑148 (MQ=255)
                                 gTTCGACACCATCGCCGCCACCCCCCGAACCACCCAGCACACGCGCGAACTGCGGTTCGCTCCTTTTTACGGGGCCTGCTGGAGGGGTTTCACCCGCTAATACTAGGTTGAACCCACCGACGAATTCGCGGGTGAACCCGCTCCTGCAgg                                                                                                                   <  2:4436251/150‑1 (MQ=255)
                                     gACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGTACCCACCGACGCCGTCGCGGGTGAACCCGCTCCGGCAGGGttg                                                                                                               >  2:3635366/1‑149 (MQ=255)
                                        atcatcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGACAAGTGGCCTGCAGGAGCGGGGTCCCCCGCGCATAAGGGGGTGGACCCACCCACGCAATCGCGGGTGGACCCGCTCCGGGCGGGTTGGcc                                                                                                            >  2:1578040/1‑150 (MQ=255)
                                            tggtcgCCACCTACCCCACGACCCCCCACACGCGCTACGTGCGGTTGGCTCCTTTTTATGGGGCCTGCCGGCGGGGGTTCCCCGGCGAATCCGCGGGTCAACCCCCCGCCGCTTGGGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCa                                                                                                        <  2:123181/148‑1 (MQ=255)
                                              gtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCCCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                      <  1:254552/118‑1 (MQ=18)
                                               tcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGACCCCCGAAAAAGAGGGTGTAACCACCGACCCAATCGCGGGGGGACCCCCGCCAGCAGGGTTCGCCCTCCGcc                                                                                                     >  2:338790/1‑150 (MQ=255)
                                                 gCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                      >  1:5087222/1‑115 (MQ=18)
                                                 gCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                      <  2:5087222/115‑1 (MQ=18)
                                                  ccaccaACCCCAACACCCACCACAACCGCCAACTGCCGTTCGCTCCTTTATAAGTGGCCTGCAGCAGCGGTTTCACCCGCGAAAACGAGGGTCAACCAACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCatt                                                                                                  <  2:5361157/150‑1 (MQ=37)
                                                  ccACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCGTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGCAAACGAGGGGGGACCCCCCCCCGCAATCGCGGGTGGACCCCCTCCGGGAGGGGTCGCCGGCAGCCatt                                                                                                  >  2:3985134/1‑150 (MQ=255)
                                                      cTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                      >  1:3117100/1‑110 (MQ=255)
                                                      cTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCAACCGCCAATACCAGGGTGTAACCACCGACGCAATCCCGGGTGAACCCGCTCCAGCAGGGGTCGGCGTCAGCCGTTGAtt                                                                                              >  2:5636896/1‑150 (MQ=255)
                                                       tACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGACCCCCGAATAACAGGGTGAACCCACCGACGCATTCGCGGGTGTACCCGCGCCTGCAGGGGTCGCCCGCAGCCGTTGATTa                                                                                             >  2:3445653/1‑150 (MQ=37)
                                                       tACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGACCCCCCAAAACGAGGGTTTACCCACCGCCCCAGTCGCGGGGGTACCCCCGCCAGCCGGGGTGGCCGGCCGCCGTTGGTTa                                                                                             >  2:346898/1‑150 (MQ=255)
                                                        aCCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGGCAGCCGTTGATTAc                                                                                            >  2:2452103/1‑150 (MQ=255)
                                                             caccccccGCACACGCGCGACCGGCCGTTCGCTCCTTTATATGTGGCCTGCAGGAGGGGGTTCACCCGCGAATACGTGGGTGAACCCACCGACGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGAc                                                                                       <  2:728528/146‑1 (MQ=255)
                                                             caGACCCAGCACCCGCGCGAACTGCCGTTCGCTCCTTGTTAAGTGGCCTGCAGGAGCGGGTTCCCCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGAc                                                                                       <  2:4614210/149‑1 (MQ=255)
                                                                   ccgCACACGCGCGACGTGCGTTTCGCTCCTTGTTACGTGGCCTGCCGGCGGGGTTCCACCCGCGAATACGCGGGTGAACCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGAc                                                                                 <  1:2375413/148‑1 (MQ=25)
                                                                     gCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGCCCCCCGAATCAGCGGGTGGACCCACCGACGCAGTCGCGGGTGGACCCGCTGCGGCCGGGTGCGCCGTCAGCCGTGGCTTACTTGTCGCTGGccg                                                                               >  2:7160398/1‑149 (MQ=255)
                                                                      cacaAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGAATACGAGGGTGTACCCACCGACGCAGTCGCGGGTGTACCCGCTCCTGCAGGGTTGGCCCTCAGCCTTTGGTTTCTTGGCGCTgtgctt                                                                              >  2:3192185/1‑145 (MQ=255)
                                                                           gcccccaCTGCCGTTCCCTCATTTATATGGGGCCGGCTGGAGCGGTTCCACCCGCCAATACGCGGGTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCtt                                                                         <  1:1129412/144‑1 (MQ=25)
                                                                             gcgAACTGCCTTTCGCTCCTTGATAAGTGGCCTGCTGGAGGGGGTTCACCCGCGAATACGTGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                                       <  1:7347620/150‑1 (MQ=255)
                                                                                      cGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCTAATACGCGGGTGAACCCACCGACGCATGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgag                                                              <  2:2269096/150‑1 (MQ=255)
                                                                                       gTTCGCTCCTTTATACGTGGCCTGCTGGCGGGGGTACCCCCGCGAATACGCGGGTGAACCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgaga                                                             <  1:3635366/150‑1 (MQ=38)
                                                                                                 ggATAAGTGGCCGGCAGTAGCGGGGTCAGCCGCTATTACTAGGGTGAACCCACCGAGGCATGGGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGGCTTTGATTAGTTGACGATGACCTTGCCTTTCATCTTCGAGATGTGGCCCgg                                                   <  2:4759394/149‑1 (MQ=255)
                                                                                                     ccggggCCTGCTGGCGCGGGTTCACCCGGTAATCCGGGGGGGACCCCCCGGCCGCATTGGGGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAg                                               <  1:4913335/146‑1 (MQ=255)
                                                                                                     aaGTGGCCTGCCGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAg                                               <  2:4692962/150‑1 (MQ=255)
                                                                                                      cgccgccggccgGCGCGGGGCCCCCCGCGACTCCGCGGGTGCACCCACCGCCGCAGGCGCGGGGGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCCTCATCGAGAGGTGGCCCGGGAAgg                                              <  2:5724388/146‑1 (MQ=255)
                                                                                                          gCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCAt                                                                                                                                            >  2:6299842/1‑52 (MQ=255)
                                                                                                            cTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCCCCGCCGCAATCGCGGGTGAACCCCCTACTGCAGGGGTCGCCGTCGGCCTGTTATTACTATACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAATGTGCaga                                        >  2:4310640/1‑150 (MQ=255)
                                                                                                                  gcgggggTCCCCCCGCGAATACGCGGGTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGc                                  <  2:3762592/146‑1 (MQ=39)
                                                                                                                   tgggggTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCt                                 <  1:2391734/147‑1 (MQ=255)
                                                                                                                   aGCGGGTACACCCCCCAAAACGCGGTACAACCCACCTACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCt                                 <  2:7039242/150‑1 (MQ=255)
                                                                                                                    gCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg                                <  2:5164822/150‑1 (MQ=255)
                                                                                                                    gCGGGTACACCCGCGAATACGAGGGTCAACCCACCGACGAATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg                                <  2:1835592/150‑1 (MQ=255)
                                                                                                                     cGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGt                               >  2:2724962/1‑150 (MQ=21)
                                                                                                                        gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                            >  1:2127910/1‑150 (MQ=255)
                                                                                                                        gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                            >  2:4657484/1‑150 (MQ=255)
                                                                                                                        gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                            >  2:929464/1‑150 (MQ=255)
                                                                                                                         ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc                           <  1:5450829/150‑1 (MQ=255)
                                                                                                                            aCCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                   <  2:1102443/139‑1 (MQ=21)
                                                                                                                            aCCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                   >  1:1102443/1‑139 (MQ=21)
                                                                                                                            aCCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc                        >  2:5134414/1‑150 (MQ=25)
                                                                                                                            aCCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc                        >  2:6134975/1‑150 (MQ=255)
                                                                                                                              ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGATGCTGTa                              <  1:3141033/142‑1 (MQ=255)
                                                                                                                              ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                              >  2:3141033/1‑142 (MQ=255)
                                                                                                                              ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  1:1824776/1‑150 (MQ=255)
                                                                                                                              ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  2:4068844/1‑150 (MQ=255)
                                                                                                                              ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  2:4251675/1‑150 (MQ=255)
                                                                                                                              ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  1:7050781/1‑150 (MQ=255)
                                                                                                                              ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCATTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  2:3350257/1‑150 (MQ=255)
                                                                                                                                gctaaTCCGTGGGTGCCCCCCCCGCCGCATGGGCGGGGGCACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCg                    <  1:3812529/147‑1 (MQ=255)
                                                                                                                                gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCg                    >  1:6424253/1‑150 (MQ=255)
                                                                                                                                gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCg                    >  1:5278640/1‑150 (MQ=255)
                                                                                                                                    aTACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCtt                <  2:5568306/150‑1 (MQ=255)
                                                                                                                                     tccGAGGGTGAACCCACCGACGCATTCGCGGGGGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTc               <  2:4429485/148‑1 (MQ=255)
                                                                                                                                        gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg            >  2:4144575/1‑150 (MQ=255)
                                                                                                                                        gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg            >  1:5508061/1‑150 (MQ=255)
                                                                                                                                        gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg            >  2:3918958/1‑150 (MQ=255)
                                                                                                                                           gtttaaCCCACCGACGCTTTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCtt         <  1:6586406/146‑1 (MQ=255)
                                                                                                                                           ggtgccccccccGCCGCATGCGGGGGTCAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCtt         <  1:1571746/144‑1 (MQ=255)
                                                                                                                                           ggTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCtt         <  2:5968311/150‑1 (MQ=255)
                                                                                                                                            gTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg        >  2:5017221/1‑150 (MQ=255)
                                                                                                                                               aaCCCACCCACGCATTCGCGGTTTACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa     <  1:1476290/150‑1 (MQ=255)
                                                                                                                                                 cccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGATGTAGTCACCACCGGCTTCCAGCTTGGAAgt   >  2:844128/1‑150 (MQ=255)
                                                                                                                                                 cccACCGACGCATTCGCGGGTGAACCAGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgt   >  2:844692/1‑150 (MQ=255)
                                                                                                                                                  ccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtg  >  2:1362728/1‑150 (MQ=255)
                                                                                                                                                  ccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtg  >  1:579750/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTG  >  NC_002947/5541921‑5542216

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: