Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,085 G→T 7.0% intergenic (+94/+313) fdnI → / ← yddM formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,0850GT 7.0% 1070.8 / 46.3 323intergenic (+94/+313)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (128/172);  new base T (5/18);  total (133/190)
Fisher's exact test for biased strand distribution p-value = 7.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 325 total aligned reads displayed.

TATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAC  >  NC_000913/1551938‑1552233
                                                                                                                                                   |                                                                                                                                                      
tATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                                                                                                      <  2:6364779/150‑1 (MQ=255)
tATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                                                                                                      <  1:2454517/150‑1 (MQ=255)
    gTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                                  <  1:6250795/150‑1 (MQ=255)
    gTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                                  <  2:3566428/150‑1 (MQ=255)
      gAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                <  2:3371251/150‑1 (MQ=255)
      gAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGccgc                                                                                                                                                <  1:5524917/150‑3 (MQ=255)
       aaaTCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                               >  1:166304/1‑150 (MQ=255)
         aTCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGccgcatc                                                                                                                                             <  1:4072870/150‑6 (MQ=37)
         aTCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGccgcatc                                                                                                                                             >  1:498562/1‑145 (MQ=37)
          tCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                            <  1:6130763/150‑1 (MQ=37)
           cGAGAAGGCAGAAGCGTATAAAGAGAGTGAAGAAGGGATTTAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGGTTTGTTCATGCGGGTTGCGGCGTGTACGCCTTATTCg                                                                                                                                           <  2:5145404/150‑1 (MQ=25)
           cGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATAGAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                           <  1:4844918/150‑1 (MQ=25)
             agaAGGCATACGCGAAAAATGATCGTGAAGACGGGATATAATCGCTTTTGAACTTTAAGCTGAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                        <  2:448643/150‑1 (MQ=18)
             agaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                         <  1:2038800/150‑1 (MQ=255)
              gaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                        <  2:819215/150‑1 (MQ=255)
              gaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                        <  2:668168/150‑1 (MQ=255)
                aGGCATAAGCGAAAAAAGTGAGTGACGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCATTTTTCTTTTTTTAGACAAGGTAGGCTTTGTTCTTGCCGGATAAGGCGTTCACGCCGCATCGGGCATg                                                                                                                                      <  1:4457932/150‑1 (MQ=255)
                 ggCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                          <  2:5194544/145‑1 (MQ=25)
                  gCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTGa‑‑a                                                                                                                                  >  1:660359/1‑150 (MQ=255)
                      aaGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGGGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAAcc                                                                                                                              <  2:3711676/150‑1 (MQ=255)
                      aaGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAAcc                                                                                                                              >  1:3730977/1‑150 (MQ=255)
                         cGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCtt                                                                                                                           <  1:2272958/150‑1 (MQ=255)
                          gAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTc                                                                                                                          >  1:2201276/1‑150 (MQ=255)
                              aaaGAGAGTGAAGAAGGGCTATCATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAAt                                                                                                                      <  1:3298131/150‑1 (MQ=255)
                              aaaGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAAt                                                                                                                      >  1:2020197/1‑150 (MQ=255)
                                 gagagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCa                                                                                                                   <  1:1405104/150‑1 (MQ=255)
                                   gagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAAt                                                                                                                 <  1:1901692/150‑1 (MQ=255)
                                    tgTGAAGATTGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAATGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCTGCATGA‑‑AAACCCTTCAAATCCAATa                                                                                                                <  1:5568297/149‑1 (MQ=255)
                                    agTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATa                                                                                                                >  1:1071540/1‑150 (MQ=255)
                                         agaagGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTg                                                                                                           <  2:6228649/150‑1 (MQ=255)
                                             gggATATACTCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGt                                                                                                       <  1:3627775/150‑1 (MQ=255)
                                                atatAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGaa                                                                                                    <  1:2886045/150‑1 (MQ=255)
                                                 tataATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAATAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGTACGCCTTATCCGGCATGT‑‑AAACCCTTCAAATCCAATAGTTTGGAGTgtag                                                                                                   >  2:1491259/1‑147 (MQ=17)
                                                 tataATCTATTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAAc                                                                                                   <  1:6458250/150‑1 (MQ=255)
                                                   taATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAAGCGATTGCAGTGAACgt                                                                                                 >  1:2772189/1‑150 (MQ=255)
                                                    aaTCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAATGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACgtg                                                                                                >  1:1251584/1‑150 (MQ=255)
                                                     aTCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACgtgt                                                                                               >  1:4167105/1‑150 (MQ=255)
                                                      tctcTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTa                                                                                              >  1:272811/1‑150 (MQ=255)
                                                       ctctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                                                                                             <  1:1811009/150‑1 (MQ=255)
                                                       ctctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                                                                                             >  1:6778699/1‑150 (MQ=255)
                                                       ctctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                                                                                             >  1:3353021/1‑150 (MQ=255)
                                                         ctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGc                                                                                           >  1:5729394/1‑150 (MQ=255)
                                                          ttttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                                                                                          >  2:1624811/1‑150 (MQ=255)
                                                          ttttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                                                                                          >  1:4566679/1‑150 (MQ=255)
                                                           tttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCt                                                                                         <  1:5618182/150‑1 (MQ=255)
                                                           tttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCt                                                                                         <  1:1624811/150‑1 (MQ=255)
                                                            ttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTg                                                                                        >  2:3366421/1‑150 (MQ=255)
                                                              gAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAt                                                                                      <  2:1155708/150‑1 (MQ=255)
                                                               aaCTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATa                                                                                     >  2:1715127/1‑150 (MQ=255)
                                                                  tttAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                                  <  1:5161144/150‑1 (MQ=40)
                                                                       gCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc                                                                             <  2:2185077/150‑1 (MQ=35)
                                                                       gCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc                                                                             >  2:5556199/1‑150 (MQ=35)
                                                                       gCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc                                                                             >  2:228561/1‑150 (MQ=35)
                                                                          gAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTAATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCa                                                                          >  1:3188049/1‑150 (MQ=18)
                                                                           aaaaTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGACCGTGTAGGCCTGATAAGCGTAGCGCAt                                                                         >  2:2078193/1‑150 (MQ=17)
                                                                           aaaaTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt                                                                         >  2:6032272/1‑150 (MQ=32)
                                                                              aTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                      <  1:872337/150‑1 (MQ=32)
                                                                                ggCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                    >  1:5317624/1‑150 (MQ=32)
                                                                                   gcTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaag                                                                 >  1:3956300/1‑149 (MQ=17)
                                                                                   gcTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGCATCCGGCATg                                                                                                                                                        <  2:6619031/65‑1 (MQ=11)
                                                                                     tGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGt                                                               >  2:3011033/1‑150 (MQ=32)
                                                                                      gTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACgt                                                                                                 <  1:3884087/115‑1 (MQ=32)
                                                                                      gTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGtt                                                              <  2:4167105/150‑1 (MQ=32)
                                                                                      gTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGtt                                                              <  2:4982156/150‑1 (MQ=32)
                                                                                      gTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGtt                                                              <  2:5772906/150‑1 (MQ=32)
                                                                                        aaaaGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGc                                                            <  2:6184990/150‑1 (MQ=32)
                                                                                           aGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGtt                                                         <  1:442441/150‑1 (MQ=32)
                                                                                            ggCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATACGCGTAGCGCATCAGGCAATGTTGCGttt                                                        <  2:4298395/150‑1 (MQ=17)
                                                                                              cgcCATTTTCATATTGTAGACACCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGt                                                      <  1:6652933/150‑1 (MQ=17)
                                                                                              cgcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                        >  1:425112/1‑70 (MQ=35)
                                                                                              cgcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                        <  2:425112/70‑1 (MQ=35)
                                                                                              cgcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGt                                                      <  1:4545236/150‑1 (MQ=32)
                                                                                               gcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                      <  1:3142153/133‑1 (MQ=31)
                                                                                               gcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                      >  2:3142153/1‑133 (MQ=31)
                                                                                               gcCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGc                                                            >  2:2487482/1‑143 (MQ=31)
                                                                                               gcCATTTTCATAGTGTAGACAACGTAGGCTGTGTTCATGCCGGATGCGGCGAGACCGCCTTCTCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGtc                                                     <  2:6716020/150‑1 (MQ=11)
                                                                                                ccATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGtca                                                    >  1:1131295/1‑150 (MQ=11)
                                                                                                       cATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCCTCAGGCAATGTTGCGTTTGTCATCAGttt                                             <  2:4642110/150‑1 (MQ=11)
                                                                                                         tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCa                                           >  2:1945929/1‑150 (MQ=32)
                                                                                                           ttGGAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACGCTTCAAATCCAATGGATTGCAGTGAACGGGTAGGCCTGATAAGCGTAGCGCATCGGGCAGTGTTGCGTTTGTCAGCAGTGTCaat                                         >  1:6053505/1‑149 (MQ=11)
                                                                                                              tAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATgg                                      <  1:3011033/150‑1 (MQ=34)
                                                                                                               aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgc                                     <  2:4034088/150‑1 (MQ=34)
                                                                                                                   aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                 <  2:1251584/150‑1 (MQ=35)
                                                                                                                     cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTa                               <  2:1201632/150‑1 (MQ=37)
                                                                                                                        aGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaaaa                            <  2:4541674/150‑1 (MQ=39)
                                                                                                                           cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGc                         <  1:980357/150‑1 (MQ=255)
                                                                                                                           cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGc                         <  1:2142586/150‑1 (MQ=255)
                                                                                                                           cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGc                         <  2:5490078/150‑1 (MQ=255)
                                                                                                                            tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                        >  1:6568036/1‑150 (MQ=255)
                                                                                                                               gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCa                     >  2:610453/1‑150 (MQ=255)
                                                                                                                                   aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCGGTTTCAAATGGCGCTGTAAAAGGCGTCAtt                   >  2:5457878/1‑150 (MQ=14)
                                                                                                                                    tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                <  1:1748972/150‑1 (MQ=255)
                                                                                                                                    tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCATCAAATCCAAAAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATAAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                >  1:1004049/1‑150 (MQ=255)
                                                                                                                                     gCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCa               <  1:5932285/150‑1 (MQ=255)
                                                                                                                                       cGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCata             >  1:4192636/1‑150 (MQ=255)
                                                                                                                                        ggATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCCGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat            >  2:2884815/1‑150 (MQ=255)
                                                                                                                                        ggATGCGGCATGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat            >  1:5441051/1‑150 (MQ=255)
                                                                                                                                            gcggcgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa        <  1:4717355/150‑1 (MQ=255)
                                                                                                                                                cgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGTAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACa    <  2:646558/150‑1 (MQ=255)
                                                                                                                                                  tGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAc  <  2:1645629/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAC  >  NC_000913/1551938‑1552233

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: