Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,542,058 A→T 20.9% intergenic (+38/+65) 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,0580AT20.9% 283.4 / 61.5 160intergenic (+38/+65)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base A (55/51);  new base T (8/21);  total (67/93)
Fisher's exact test for biased strand distribution p-value = 2.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

Only 100 of 166 total aligned reads displayed.

GCAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAG  >  NC_002947/5541909‑5542206
                                                                                                                                                     |                                                                                                                                                    
gcAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGa                                                                                                                                                      <  2:1511907/150‑1 (MQ=255)
gcAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGa                                                                                                                                                      <  2:1601160/150‑1 (MQ=255)
gcAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGa                                                                                                                                                      <  1:4558372/150‑1 (MQ=255)
gcAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGa                                                                                                                                                      <  2:229714/150‑1 (MQ=255)
 cAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAg                                                                                                                                                     <  2:4144673/150‑1 (MQ=255)
  aGATCGCTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                    <  1:968060/150‑1 (MQ=255)
  aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                                    <  2:3262030/150‑1 (MQ=255)
  aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                                    <  2:3261561/150‑1 (MQ=255)
  aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                    <  1:5016761/150‑1 (MQ=255)
   gATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGAGAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAggg                                                                                                                                                   >  2:6756091/1‑150 (MQ=255)
    aTCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                                  <  2:1806753/150‑1 (MQ=255)
       gATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGaa                                                                                                                                               <  1:3893897/150‑1 (MQ=255)
        aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCTCGTATTCGCGGGTGAAc                                                                                                                                              <  1:2134784/150‑1 (MQ=255)
        aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAc                                                                                                                                              <  1:3563965/150‑1 (MQ=255)
            cTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCCCGACAACCAGGGTGAACCCAc                                                                                                                                          >  2:4619169/1‑150 (MQ=255)
             ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                                  >  2:6240137/1‑141 (MQ=255)
             ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                                  <  1:6240137/141‑1 (MQ=255)
             ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGGCCCGCGCCTAGGTGGGGGGACCCAcc                                                                                                                                         >  1:4162532/1‑150 (MQ=255)
              tCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGACCCGCCAAAAACGGGGTGGACCCCCCg                                                                                                                                        >  2:2566368/1‑150 (MQ=255)
                  gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCCACGc                                                                                                                                    >  2:1374341/1‑150 (MQ=255)
                      gggCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTc                                                                                                                                <  2:7021840/150‑1 (MQ=255)
                        gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                              <  2:5549956/150‑1 (MQ=255)
                        gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                              <  2:5889380/150‑1 (MQ=255)
                        gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                              >  2:612642/1‑150 (MQ=255)
                         cAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgcg                                                                                                                             <  2:5678685/150‑1 (MQ=255)
                          cgCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  1:795059/149‑1 (MQ=255)
                          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTCCCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  2:4038140/150‑1 (MQ=255)
                          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  2:5334277/150‑1 (MQ=25)
                          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  2:1717864/150‑1 (MQ=25)
                          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  2:522313/150‑1 (MQ=25)
                          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  2:542946/150‑1 (MQ=25)
                          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  2:2294925/150‑1 (MQ=25)
                          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  2:4370029/150‑1 (MQ=255)
                          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  1:5512352/150‑1 (MQ=255)
                          aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  1:4956474/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)
                            ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGt                                                                                                                          <  2:4528426/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                                                                                                                       <  1:4648403/150‑1 (MQ=255)
                                      aGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCGAAAACGAGGGTGAAACCAACCACGCAATCGCGGGTGAAcccccgc                                                                                                                >  2:184499/1‑146 (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)
                                                   catcatCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGAATACGAGGGTGGACCCACCGACGCAATCGCGGGTGAACCCCCGCCGGCCggggtggc                                                                                                   >  2:5256925/1‑145 (MQ=255)
                                                    atcatcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGACAAGTGGCCTGCAGGAGCGGGGTCCCCCGCGCATAAGGGGGTGGACCCACCCACGCAATCGCGGGTGGACCCGCTCCGGGCGGGTTGGcc                                                                                                  >  2:1578040/1‑150 (MQ=255)
                                                      cagcggcgCCACCTCCCCCACGCCCCCGCCCCCGCGCGACCGGCCGTGCGCTCCTTGATATGGGGCCTGCAGGAGGGGGTGCACCCGCTAAGACGCGGGGGACCCCACCGCCGCATGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGt                                                                                                <  2:5119424/144‑1 (MQ=255)
                                                          gtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCCCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  1:254552/118‑1 (MQ=18)
                                                           tcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGACCCCCGAAAAAGAGGGTGTAACCACCGACCCAATCGCGGGGGGACCCCCGCCAGCAGGGTTCGCCCTCCGcc                                                                                           >  2:338790/1‑150 (MQ=255)
                                                             gCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  2:5087222/115‑1 (MQ=18)
                                                              ccaccaACCCCAACACCCACCACAACCGCCAACTGCCGTTCGCTCCTTTATAAGTGGCCTGCAGCAGCGGTTTCACCCGCGAAAACGAGGGTCAACCAACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCatt                                                                                        <  2:5361157/150‑1 (MQ=37)
                                                              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)
                                                                   tACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGACCCCCCAAAACGAGGGTTTACCCACCGCCCCAGTCGCGGGGGTACCCCCGCCAGCCGGGGTGGCCGGCCGCCGTTGGTTa                                                                                   >  2:346898/1‑150 (MQ=255)
                                                                    aCCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGGCAGCCGTTGATTAc                                                                                  >  2:2452103/1‑150 (MQ=255)
                                                                        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)
                                                                                  cacaAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGAATACGAGGGTGTACCCACCGACGCAGTCGCGGGTGTACCCGCTCCTGCAGGGTTGGCCCTCAGCCTTTGGTTTCTTGGCGCTgtgctt                                                                    >  2:3192185/1‑145 (MQ=255)
                                                                                    caAGCGCGAACGGCCGTTCGCTCCTTTATATGTGGCCTGCAGGAGGGGGTTCACCCGCGAATACGAGGGTTAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGc                                                                  <  2:6366829/150‑1 (MQ=38)
                                                                                       gcccccaCTGCCGTTCCCTCATTTATATGGGGCCGGCTGGAGCGGTTCCACCCGCCAATACGCGGGTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCtt                                                               <  1:1129412/144‑1 (MQ=25)
                                                                                          cgAACTGCCGTTCGCTCCTTGATACGGGGCCGGCAGGAGCGGGTCCACCCGCGAATACGTGGGTGACCCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                            <  2:2259685/150‑1 (MQ=37)
                                                                                                gCCGTTCGCTCCTTGACACGTGGCCTGCAGGAGCGGGTACACCCGCGACTCCGAGGGTGAACCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCATTGATTACTTGACGATGACCTTGCCTTTCATCATCg                                                      <  2:4270051/150‑1 (MQ=38)
                                                                                                  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)
                                                                                                                  cggggccggccGGCGCGGGTTCCCCCGCGAATCCGCGGGTGCCCCCACCGACGCCGGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAgg                                    <  2:4678995/147‑1 (MQ=255)
                                                                                                                  cgTGGCCGGCCGGCGGGGGTTCACCCGCTAATCCGTGGGAGCCCCCCCCGCCGATTTCGGGGGTGCACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAgg                                    <  1:1607902/149‑1 (MQ=255)
                                                                                                                      gCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCAt                                                                                                                                  <  1:6299842/52‑1 (MQ=255)
                                                                                                                              gcgggggTCCCCCCGCGAATACGCGGGTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGc                        <  2:3762592/146‑1 (MQ=39)
                                                                                                                               tgggggTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCt                       <  1:2391734/147‑1 (MQ=255)
                                                                                                                                gggggTTCACCCGCGAATCCGCGGGTGACCCCACCGACGCATGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg                      <  1:6756091/148‑1 (MQ=255)
                                                                                                                                gCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg                      <  2:5164822/150‑1 (MQ=255)
                                                                                                                                 cGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGt                     >  2:2724962/1‑150 (MQ=21)
                                                                                                                                    gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                  >  2:929464/1‑150 (MQ=255)
                                                                                                                                    gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                  >  2:4657484/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:1435426/1‑150 (MQ=25)
                                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCa                                                                                                           <  1:5174832/55‑1 (MQ=255)
                                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAggc                                                                                                        <  1:7184210/58‑2 (MQ=255)
                                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGATGCTGTa                    <  1:3141033/142‑1 (MQ=255)
                                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:1824776/1‑150 (MQ=255)
                                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:7050781/1‑150 (MQ=255)
                                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:3566633/1‑150 (MQ=255)
                                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:7329282/1‑150 (MQ=255)
                                                                                                                                          ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  2:4068844/1‑150 (MQ=255)
                                                                                                                                                    gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGACGGTGCAGGAGAAGCTGTAGTCACCACCGGCTTCCAg  >  2:7032291/1‑150 (MQ=255)
                                                                                                                                                    gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg  >  2:3918958/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAG  >  NC_002947/5541909‑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: