Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,385,160 G→T 5.5% intergenic (+74/+34) recN → / ← fur DNA repair/recombination protein/ferric iron uptake global transcriptional repressor Fur

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,385,1600GT 5.5% 871.8 / 26.1 256intergenic (+74/+34)recN/furDNA repair/recombination protein/ferric iron uptake global transcriptional repressor Fur
Reads supporting (aligned to +/- strand):  ref base G (124/118);  new base T (5/9);  total (129/127)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 258 total aligned reads displayed.

GATGCTCGGCGGGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGACTT  >  NC_002947/5385011‑5385309
                                                                                                                                                     |                                                                                                                                                     
gATGCTCGGCGGGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAACGCACAAAGGCGACCCg                                                                                                                                                       >  2:5296569/1‑150 (MQ=255)
 aagCTCGTCGCGATCGACCTGACCAACGAATCCCTGGCCCATGCACGCAAGATGGTGCTGATCGGCAAGGCCTGAGCCTTGCATTTCCAGCACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGa                                                                                                                                                      <  2:3952660/148‑1 (MQ=255)
      cggcggGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTc                                                                                                                                                 <  1:4951849/150‑1 (MQ=255)
      cggcggGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTc                                                                                                                                                 <  2:4099726/150‑1 (MQ=255)
      cggcggGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTc                                                                                                                                                 <  2:1244957/150‑1 (MQ=255)
       ggcggGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCg                                                                                                                                                >  2:5204494/1‑150 (MQ=255)
          gggATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCt                                                                                                                                             >  2:5447782/1‑150 (MQ=255)
          gggATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCt                                                                                                                                             <  2:2391588/150‑1 (MQ=255)
           ggATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCtt                                                                                                                                            <  1:419751/150‑1 (MQ=255)
            gATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCTTCGGGTCGCCttt                                                                                                                                           <  1:3497754/150‑1 (MQ=255)
            gATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCTTCGGGTCGCCttt                                                                                                                                           <  1:817623/150‑1 (MQ=255)
            gATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCttt                                                                                                                                           <  1:3779245/150‑1 (MQ=255)
              tCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGc                                                                                                                                         <  2:2567522/150‑1 (MQ=255)
              tCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGAGGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCTTCGGGTCGCCTTTgt                                                                                                                                         <  1:2791855/150‑2 (MQ=255)
                gACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTa                                                                                                                                       >  2:660442/1‑150 (MQ=255)
                 aCCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTAt                                                                                                                                      >  1:681444/1‑150 (MQ=255)
                  ccTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATc                                                                                                                                     >  2:5073328/1‑150 (MQ=255)
                        cAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAAc                                                                                                                               >  1:1299904/1‑150 (MQ=255)
                            gAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGaaa                                                                                                                           <  1:4587810/150‑1 (MQ=255)
                                cccTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAAcc                                                                                                                       <  1:1739783/150‑1 (MQ=255)
                                 ccTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCg                                                                                                                      >  2:5612990/1‑150 (MQ=255)
                                  cTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGt                                                                                                                     <  1:5269161/150‑1 (MQ=255)
                                  cTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGt                                                                                                                     <  1:5495624/150‑1 (MQ=255)
                                  cTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGt                                                                                                                     <  1:5557577/150‑1 (MQ=255)
                                  cTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGt                                                                                                                     >  1:4731055/1‑150 (MQ=255)
                                  cTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGt                                                                                                                     <  2:327165/150‑1 (MQ=255)
                                  cTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGt                                                                                                                     <  2:4808592/150‑1 (MQ=255)
                                   tGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGtt                                                                                                                    <  2:2145244/150‑1 (MQ=255)
                                     gCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTAc                                                                                                                  <  2:1512355/150‑1 (MQ=255)
                                       ccATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActt                                                                                                                >  2:4192219/1‑150 (MQ=255)
                                        cATGCTCGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActtc                                                                                                               <  2:3775856/150‑1 (MQ=255)
                                        cATGCACGCAAGATGGTGGTGAGCGGCTAGGCCTGAGCCTTGCTATTCCACGACGATTCGCGGGCACAGACAACACAACTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActtc                                                                                                               <  1:3682832/150‑1 (MQ=255)
                                         aTGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActtct                                                                                                              >  2:4436920/1‑150 (MQ=255)
                                           gCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActtcttc                                                                                                            <  1:410523/150‑1 (MQ=255)
                                             aCGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActtcttctt                                                                                                          >  1:2617100/1‑150 (MQ=255)
                                             aCGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActtcttctt                                                                                                          <  1:3580030/150‑1 (MQ=255)
                                             aCGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTActtcttctt                                                                                                          <  2:882127/150‑1 (MQ=255)
                                                 aaGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGc                                                                                                      <  2:3103039/150‑1 (MQ=255)
                                                      ggtggtGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTa                                                                                                 >  2:3622877/1‑150 (MQ=255)
                                                       gtggtgAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTGCGCACGTAc                                                                                                <  2:2000483/150‑1 (MQ=255)
                                                         ggtgAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAg                                                                                              <  2:798930/150‑1 (MQ=255)
                                                          gtgAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAgg                                                                                             >  2:4138982/1‑150 (MQ=255)
                                                          gtgAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAgg                                                                                             >  2:1370523/1‑150 (MQ=255)
                                                           tgAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGa                                                                                            <  1:4693974/150‑1 (MQ=255)
                                                           tgAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGa                                                                                            >  1:3040206/1‑150 (MQ=255)
                                                            gAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTGACTTCTTCTTACGCACGTACAGGAc                                                                                           <  1:2022067/150‑1 (MQ=255)
                                                                 gCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGAt                                                                                      >  2:224709/1‑150 (MQ=255)
                                                                  cAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGAtt                                                                                     <  1:5153250/150‑1 (MQ=255)
                                                                  cAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGAtt                                                                                     >  1:2194057/1‑150 (MQ=255)
                                                                          gAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAg                                                                                              <  1:1416558/133‑1 (MQ=255)
                                                                           aGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTccac                                                                            <  1:3883569/150‑1 (MQ=255)
                                                                            gCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCACAATCTGGTcctgt                                                                           <  1:1015828/150‑4 (MQ=255)
                                                                              cTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAg                                                                         >  2:3437284/1‑150 (MQ=255)
                                                                                tGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCt                                                                       <  2:75345/150‑1 (MQ=255)
                                                                                 gCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTc                                                                      <  1:5416382/150‑1 (MQ=255)
                                                                                      tCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGc                                                                 <  1:2716325/150‑1 (MQ=255)
                                                                                        cAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCa                                                               >  2:4846046/1‑150 (MQ=255)
                                                                                         aGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCAt                                                              >  2:514145/1‑150 (MQ=255)
                                                                                          ggACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATg                                                             <  2:1676707/150‑1 (MQ=255)
                                                                                              gATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACGTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCg                                                         >  2:4795739/1‑150 (MQ=255)
                                                                                                   gcgGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGAGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCAACAGCTCGAAGCCATGCTCGGCTAc                                                    >  2:3857459/1‑150 (MQ=255)
                                                                                                    cgGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACg                                                   >  2:5264473/1‑150 (MQ=255)
                                                                                                    cgGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGAACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGATACg                                                   >  2:4384192/1‑150 (MQ=255)
                                                                                                     gggCACAGACAACACATCTTTGAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGa                                                  <  2:3991541/150‑1 (MQ=255)
                                                                                                       gCACAGACAACACATGTTTTAGCCCCATAAAGCACAAAGGAGACCAGTAGGTCGGGTTTGCTATCCATAACGAAAAACCGTTTCTTCTTCTTAAGCACGTAAATGACGAGTTTGTGGTGCACCAGCTCGTAGCCATGCTCGGCTACGAtt                                                >  2:3318655/1‑150 (MQ=255)
                                                                                                       gCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGAtt                                                >  1:4610868/1‑150 (MQ=255)
                                                                                                       gCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGAtt                                                <  1:1535499/150‑1 (MQ=255)
                                                                                                       gCACAGACAACACATCTTTAAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGAtt                                                >  2:1129652/1‑150 (MQ=255)
                                                                                                         acaGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTc                                              >  1:4277782/1‑150 (MQ=255)
                                                                                                                 cacaTCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTCCAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGgc                                      >  2:1268754/1‑150 (MQ=255)
                                                                                                                  acaTCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGgcg                                     >  2:1944879/1‑150 (MQ=255)
                                                                                                                   caTCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGgcgc                                    <  2:4513655/150‑1 (MQ=255)
                                                                                                                    aTCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCt                                   >  1:890987/1‑150 (MQ=255)
                                                                                                                      catttAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTtc                                 <  2:4004638/148‑1 (MQ=255)
                                                                                                                       ttttAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATCACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTtct                                <  1:2189105/150‑1 (MQ=255)
                                                                                                                           aGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGAt                            >  2:1190856/1‑150 (MQ=255)
                                                                                                                             ccccAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCAACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGAttt                          >  1:597597/1‑150 (MQ=255)
                                                                                                                               ccAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCg                        >  2:3904138/1‑150 (MQ=255)
                                                                                                                               ccAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCg                        >  1:3122992/1‑150 (MQ=255)
                                                                                                                               ccAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCg                        >  2:1553188/1‑150 (MQ=255)
                                                                                                                                cAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCgg                       >  2:4997909/1‑150 (MQ=255)
                                                                                                                                 aGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTtct                                >  1:5000970/1‑140 (MQ=255)
                                                                                                                                       cacaAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCAt                <  1:5296569/150‑1 (MQ=255)
                                                                                                                                       cacaAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCAt                >  2:3501770/1‑150 (MQ=255)
                                                                                                                                        acaAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATg               >  2:2949818/1‑150 (MQ=255)
                                                                                                                                        acaAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATg               >  1:1971443/1‑150 (MQ=255)
                                                                                                                                          aaaGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGaa             >  1:1873355/1‑150 (MQ=255)
                                                                                                                                          aaaGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGaa             >  1:3991313/1‑150 (MQ=255)
                                                                                                                                           aaGGCGACCTTCGGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAAt            >  1:4351319/1‑150 (MQ=255)
                                                                                                                                           aaGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAAt            >  1:4679576/1‑150 (MQ=255)
                                                                                                                                            aGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAAtt           >  1:4807437/1‑150 (MQ=255)
                                                                                                                                              gCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCg         <  1:3453313/150‑1 (MQ=255)
                                                                                                                                               cGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGa        <  2:5115641/150‑1 (MQ=255)
                                                                                                                                               cGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGa        <  1:1781471/150‑1 (MQ=255)
                                                                                                                                               cGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGa        <  2:1334136/150‑1 (MQ=255)
                                                                                                                                                gACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGAt       >  1:5201844/1‑150 (MQ=255)
                                                                                                                                                 aCCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATg      <  2:1905768/150‑1 (MQ=255)
                                                                                                                                                  cccGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGa     >  2:5553079/1‑150 (MQ=255)
                                                                                                                                                  cccGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGa     <  1:5447782/150‑1 (MQ=255)
                                                                                                                                                     gAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGACtt  >  1:953331/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATGCTCGGCGGGATCGACCTGACCAAGGAATCCCTGGCCCATGCACGCAAGATGGTGGTGAGCGGCAAGGCCTGAGCCTTGCTTTTCCAGGACGATTCGCGGGCACAGACAACACATCTTTTAGCCCCAGAAAGCACAAAGGCGACCCGAAGGTCGCCTTTGCTATCCATAACGAAAAACCGTTACTTCTTCTTACGCACGTACAGGACCAGATTGTGGTCCACCAGCTCGAAGCCATGCTCGGCTACGATTTCGCGCTGGCGCTTCTCGATTTCGGCATCCATGAATTCGATGACTT  >  NC_002947/5385011‑5385309

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