Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,542,059 G→C 13.8% intergenic (+39/+64) 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,0590GC13.8% 423.3 / 56.4 153intergenic (+39/+64)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base G (58/73);  new base C (7/14);  total (66/87)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-02

Only 100 of 161 total aligned reads displayed.

CAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAG  >  NC_002947/5541910‑5542206
                                                                                                                                                     |                                                                                                                                                   
cAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAg                                                                                                                                                     <  2:4144673/150‑1 (MQ=255)
 aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                                    <  2:3261561/150‑1 (MQ=255)
 aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                    <  1:5016761/150‑1 (MQ=255)
 aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                    <  1:4695624/150‑1 (MQ=255)
  gATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGAGAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAggg                                                                                                                                                   >  2:6756091/1‑150 (MQ=255)
      gATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGaa                                                                                                                                               <  2:6004581/150‑1 (MQ=255)
      gATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGaa                                                                                                                                               <  1:3893897/150‑1 (MQ=255)
       aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCTCGTATTCGCGGGTGAAc                                                                                                                                              <  1:2134784/150‑1 (MQ=255)
       aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAc                                                                                                                                              <  2:5549213/150‑1 (MQ=255)
       aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAc                                                                                                                                              <  1:3205221/150‑1 (MQ=255)
       aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAc                                                                                                                                              <  1:3564079/150‑1 (MQ=255)
           cTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCCCGACAACCAGGGTGAACCCAc                                                                                                                                          >  2:4619169/1‑150 (MQ=255)
            ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGTGGGTTCACCCTCGTATTCGCGGGTGAAcccgct                                                                                                                                         <  1:6985110/150‑4 (MQ=255)
            ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                                  <  1:6240137/141‑1 (MQ=255)
            ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                                  >  2:6240137/1‑141 (MQ=255)
             tCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGACCCGCCAAAAACGGGGTGGACCCCCCg                                                                                                                                        >  2:2566368/1‑150 (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)
                         aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTCCCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  2:4038140/150‑1 (MQ=255)
                         aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  2:1717864/150‑1 (MQ=25)
                         aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  2:5334414/150‑1 (MQ=25)
                         aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  1:7301201/150‑1 (MQ=25)
                         aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  2:2294925/150‑1 (MQ=25)
                         aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  2:542946/150‑1 (MQ=25)
                         aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  2:655004/150‑1 (MQ=255)
                         aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  2:4370029/150‑1 (MQ=255)
                         aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  1:4956411/150‑1 (MQ=255)
                         aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATACGCgg                                                                                                                            <  2:1662057/150‑1 (MQ=255)
                           ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGt                                                                                                                          <  2:4528426/150‑1 (MQ=255)
                           ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGt                                                                                                                          <  2:7233694/150‑1 (MQ=255)
                              gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGACCTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGACTGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGaa                                                                                                                       <  1:4052597/150‑1 (MQ=21)
                              gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGaa                                                                                                                       <  1:5476027/150‑1 (MQ=25)
                              gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  2:4270111/145‑1 (MQ=255)
                              gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            >  1:4270111/1‑145 (MQ=255)
                              gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGaa                                                                                                                       >  2:4648403/1‑150 (MQ=255)
                              gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGaa                                                                                                                       <  1:4648403/150‑1 (MQ=255)
                                  aCCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCCACGCATTCGCGGGTGAACCCg                                                                                                                   >  2:488346/1‑150 (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)
                                            gTTCGACACCATCGCCGCCACCCCCCGAACCACCCAGCACACGCGCGAACTGCGGTTCGCTCCTTTTTACGGGGCCTGCTGGAGGGGTTTCACCCGCTAATACTAGGTTGAACCCACCGACGAATTCGCGGGTGAACCCGCTCCTGCAgg                                                                                                         <  2:4436251/150‑1 (MQ=255)
                                                gACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGTACCCACCGACGCCGTCGCGGGTGAACCCGCTCCGGCAGGGttg                                                                                                     >  2:3635366/1‑149 (MQ=255)
                                                cacacCATCGCCGCCACCAACCGCACGACCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGCCGCGGGTTCACCCGCTAATACGCGGGTGAACCCACCGACGCATCGGCGGGTGACCCCGCTCCTGCAGGGTTc                                                                                                     <  2:1645715/149‑1 (MQ=255)
                                                  catcatCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGAATACGAGGGTGGACCCACCGACGCAATCGCGGGTGAACCCCCGCCGGCCggggtggc                                                                                                   >  2:5256925/1‑145 (MQ=255)
                                                     cagcggcgCCACCTCCCCCACGCCCCCGCCCCCGCGCGACCGGCCGTGCGCTCCTTGATATGGGGCCTGCAGGAGGGGGTGCACCCGCTAAGACGCGGGGGACCCCACCGCCGCATGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGt                                                                                                <  2:5119424/144‑1 (MQ=255)
                                                      atcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCCCGCATACGGGGGTGAAACCACCGGCGCATTCGCGGGTGATCCCGCGCCTGCCGGGTTGGCCggc                                                                                               >  1:4171951/1‑148 (MQ=255)
                                                       tcgtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGCGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCa                                                                                              >  2:3214287/1‑150 (MQ=255)
                                                         gtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCCCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  1:254552/118‑1 (MQ=18)
                                                         gtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCCCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            >  2:254552/1‑118 (MQ=18)
                                                            gCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            >  1:5087222/1‑115 (MQ=18)
                                                             ccACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCGTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGCAAACGAGGGGGGACCCCCCCCCGCAATCGCGGGTGGACCCCCTCCGGGAGGGGTCGCCGGCAGCCatt                                                                                        >  2:3985134/1‑150 (MQ=255)
                                                                 cTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  2:3117100/110‑1 (MQ=255)
                                                                 cTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            >  1:3117100/1‑110 (MQ=255)
                                                                 cTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCAACCGCCAATACCAGGGTGTAACCACCGACGCAATCCCGGGTGAACCCGCTCCAGCAGGGGTCGGCGTCAGCCGTTGAtt                                                                                    >  2:5636896/1‑150 (MQ=255)
                                                                  tACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGACCCCCGAATAACAGGGTGAACCCACCGACGCATTCGCGGGTGTACCCGCGCCTGCAGGGGTCGCCCGCAGCCGTTGATTa                                                                                   >  2:3445653/1‑150 (MQ=37)
                                                                  tACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  2:3980811/109‑1 (MQ=17)
                                                                       cAAGACCCTGCACACGCGCGAACTGCCGTTCGATCCTTGATACGGGGCCTGCTGGAGCGGGTTCACCCGCGAATACTTGGGTTAACCCACCTACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGa                                                                              <  2:2345496/150‑1 (MQ=38)
                                                                        caccccccGCACACGCGCGACCGGCCGTTCGCTCCTTTATATGTGGCCTGCAGGAGGGGGTTCACCCGCGAATACGTGGGTGAACCCACCGACGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGAc                                                                             <  2:728528/146‑1 (MQ=255)
                                                                              ccgCACACGCGCGACGTGCGTTTCGCTCCTTGTTACGTGGCCTGCCGGCGGGGTTCCACCCGCGAATACGCGGGTGAACCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGAc                                                                       <  1:2375413/148‑1 (MQ=25)
                                                                                gCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGCCCCCCGAATCAGCGGGTGGACCCACCGACGCAGTCGCGGGTGGACCCGCTGCGGCCGGGTGCGCCGTCAGCCGTGGCTTACTTGTCGCTGGccg                                                                     >  2:7160398/1‑149 (MQ=255)
                                                                                   caAGCGCGAACGGCCGTTCGCTCCTTTATATGTGGCCTGCAGGAGGGGGTTCACCCGCGAATACGAGGGTTAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGc                                                                  <  2:6366829/150‑1 (MQ=38)
                                                                                      gcccccaCTGCCGTTCCCTCATTTATATGGGGCCGGCTGGAGCGGTTCCACCCGCCAATACGCGGGTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCtt                                                               <  1:1129412/144‑1 (MQ=25)
                                                                                        gcgAACTGCCGTTCGCTCCTTTATATGTGGCCTGCTGGAGCGGGTTCACCCGCTAATACGCGGGTTAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                             <  2:3393853/150‑1 (MQ=39)
                                                                                         cgAACTGCCGTTCGCTCCTTGATACGGGGCCGGCAGGAGCGGGTCCACCCGCGAATACGTGGGTGACCCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                            <  2:2259685/150‑1 (MQ=37)
                                                                                               gCCGTTCGCTCCTTGACACGTGGCCTGCAGGAGCGGGTACACCCGCGACTCCGAGGGTGAACCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCATTGATTACTTGACGATGACCTTGCCTTTCATCATCg                                                      <  2:4270051/150‑1 (MQ=38)
                                                                                                ccGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCCCCAAAAAGATGGTTAAACCAACCACGCAATAGCGGGTGAACCCCCACCGGCAGGGTTCGCCGGCACCCATTTATTAAATTAAGATGAACATTCCATTAATCATAga                                                     >  2:3495908/1‑150 (MQ=255)
                                                                                                 cGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCTAATACGCGGGTGAACCCACCGACGCATGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgag                                                    <  2:2269096/150‑1 (MQ=255)
                                                                                                            ggATAAGTGGCCGGCAGTAGCGGGGTCAGCCGCTATTACTAGGGTGAACCCACCGAGGCATGGGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGGCTTTGATTAGTTGACGATGACCTTGCCTTTCATCTTCGAGATGTGGCCCgg                                         <  2:4759394/149‑1 (MQ=255)
                                                                                                                 cggggccggccGGCGCGGGTTCCCCCGCGAATCCGCGGGTGCCCCCACCGACGCCGGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAgg                                    <  2:4678995/147‑1 (MQ=255)
                                                                                                                 cgccgccggccgGCGCGGGGCCCCCCGCGACTCCGCGGGTGCACCCACCGCCGCAGGCGCGGGGGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCCTCATCGAGAGGTGGCCCGGGAAgg                                    <  2:5724388/146‑1 (MQ=255)
                                                                                                                 cgTGGCCGGCCGGCGGGGGTTCACCCGCTAATCCGTGGGAGCCCCCCCCGCCGATTTCGGGGGTGCACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAgg                                    <  1:1607902/149‑1 (MQ=255)
                                                                                                                             gcgggggTCCCCCCGCGAATACGCGGGTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGc                        <  2:3762592/146‑1 (MQ=39)
                                                                                                                             gAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCCACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCGTTGATTACTTGTCGATGACCTTGCCTTTCATCATAGAGATGTGGCCCGGGAAGGTGCAGAAGAAGc                        >  2:2757502/1‑150 (MQ=255)
                                                                                                                              tgggggTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCt                       <  1:2391734/147‑1 (MQ=255)
                                                                                                                               gCGGGTACACCCGCGAATACGAGGGTCAACCCACCGACGAATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg                      <  2:1835592/150‑1 (MQ=255)
                                                                                                                                   gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                  >  2:929464/1‑150 (MQ=255)
                                                                                                                                       aCCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                         >  1:1102443/1‑139 (MQ=21)
                                                                                                                                       aCCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                         <  2:1102443/139‑1 (MQ=21)
                                                                                                                                       aCCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc              >  2:5134414/1‑150 (MQ=25)
                                                                                                                                       aCCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc              >  2:1435426/1‑150 (MQ=25)
                                                                                                                                       aCCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc              >  2:6134975/1‑150 (MQ=255)
                                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCa                                                                                                           >  2:5174832/1‑55 (MQ=255)
                                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCCTCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  2:1085133/1‑150 (MQ=255)
                                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                    >  2:3141033/1‑142 (MQ=255)
                                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:7329282/1‑150 (MQ=255)
                                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:7050781/1‑150 (MQ=255)
                                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:3566633/1‑150 (MQ=255)
                                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:1824776/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)
                                                                                                                                                   gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGACGGTGCAGGAGAAGCTGTAGTCACCACCGGCTTCCAg  >  2:7032291/1‑150 (MQ=255)
                                                                                                                                                   gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg  >  1:5508061/1‑150 (MQ=255)
                                                                                                                                                   gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg  >  2:4144575/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAG  >  NC_002947/5541910‑5542206

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