Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,191 C→A 11.8% intergenic (+200/+207) 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,1910CA11.8% 215.4 / 4.2 96intergenic (+200/+207)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (27/54);  new base A (7/6);  total (36/60)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCA  >  NC_000913/1552042‑1552339
                                                                                                                                                     |                                                                                                                                                    
atatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                                                                                                                                      >  1:141172/1‑150 (MQ=31)
atatTGTAGACAACGTAGGCTTTGTGCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATCAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                                                                                                                                      <  2:35163/150‑1 (MQ=14)
    tGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAAt                                                                                                                                                  <  1:1461646/150‑1 (MQ=34)
     gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATg                                                                                                                                                 <  2:92762/150‑1 (MQ=34)
      tAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATgg                                                                                                                                                <  1:1736741/150‑1 (MQ=34)
        gACCACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcg                                                                                                                                              <  2:880044/150‑1 (MQ=21)
        gACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcg                                                                                                                                              <  1:1144018/150‑1 (MQ=34)
         aCAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcgc                                                                                                                                             <  1:1619610/150‑1 (MQ=34)
           aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGATCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCtc                                                                                                                                           <  1:611795/150‑2 (MQ=21)
            aCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                                                                                                                          <  1:31955/150‑1 (MQ=35)
              gTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGAGTTTGTCATCAGTTAAAAATGGCTCTGTaa                                                                                                                                        >  2:330564/1‑150 (MQ=11)
                aGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCTGCATGAAAACCCTTCAAATCCAATAGATTGCAGTTAACGTGTAGGCCTGATAATAGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTATaaaa                                                                                                                                      >  2:648928/1‑150 (MQ=11)
                 ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTTAACGTGTAGGACTTATAAGCGTAGCGAATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAg                                                                                                                                     <  1:91378/150‑1 (MQ=255)
                 ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCCGCAGGCAATGTTGCGTTTGTCATCAGTTTGAAATGGCGCTGTAAAAg                                                                                                                                     >  1:1388870/1‑150 (MQ=255)
                 ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAg                                                                                                                                     <  2:1079044/150‑1 (MQ=40)
                 ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAg                                                                                                                                     <  1:935796/150‑1 (MQ=40)
                 ggCTTTGTTCATGCCGGATGAGGCGTGAACGCCTTATCCGGCATGAAAACACTTCAAATCCAATAGATTGCAGTTAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGAAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAg                                                                                                                                     <  2:292202/150‑1 (MQ=14)
                  gatttGTTCATGCCGGATGCGTCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAgg                                                                                                                                    <  1:1391432/148‑1 (MQ=255)
                  gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAgg                                                                                                                                    <  1:852715/150‑1 (MQ=40)
                      tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTc                                                                                                                                <  2:425927/150‑1 (MQ=255)
                       gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAACCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTTAAATGGCGCTGTAGAAggaggca                                                                                                                               >  1:1322306/1‑145 (MQ=255)
                       gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCa                                                                                                                               <  1:675331/150‑1 (MQ=255)
                       gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCa                                                                                                                               <  1:666556/150‑1 (MQ=255)
                           aTGCCGGATGCGGCTTGAACGCAGAATCCGTCATGAAAACCCTTCAAAACCTATATATAGCATTGAACGTGTAGCCCAGATAAGCGTAGCTCAGCAGGCAATGTTGCGTTTGTCATCAGCTTCAAATGGCGCTGTAAAAGGCGTCAtttt                                                                                                                           <  1:323293/150‑1 (MQ=14)
                           aTGCAGTATTCGGCGTGAACGCCTTATCCGTCATGAAAACCCTTCAAAACCAATAGATTGCAGTTAACGTGTAGGCCTGATAAGCTTAGCGCTTCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtttt                                                                                                                           <  2:951029/150‑1 (MQ=14)
                            tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                                                                                                          <  2:727025/150‑1 (MQ=255)
                            tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                                                                                                          <  2:141172/150‑1 (MQ=255)
                            tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                                                                                                          <  1:693870/150‑1 (MQ=255)
                             gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCa                                                                                                                         <  1:1381630/150‑1 (MQ=255)
                              cctgATGCGGCGTTAACGCCTAATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGTACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                                                                                                        <  1:1072921/147‑1 (MQ=17)
                              ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                                                                                                        <  2:576882/150‑1 (MQ=255)
                               cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCata                                                                                                                       <  2:1583689/150‑1 (MQ=255)
                               cGGATGCCGCGTGAACTCCTTATCCTGCATGAAAATCCATCAACTCCAATAGATATCAGTGAACGTGTAGTCTTGATAAGATTAGCGCATCAGGCAATTTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAATAGTCGTCATTTTCata                                                                                                                       <  2:1775695/150‑1 (MQ=25)
                                  aTGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTg                                                                                                                    <  1:1725819/150‑1 (MQ=255)
                                   tGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                                                                                                                   >  1:453122/1‑150 (MQ=255)
                                     cggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAATCGTAGCGCATCAGGCAATGTTGCGTTTGTAATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAg                                                                                                                 <  2:1180952/150‑1 (MQ=255)
                                       gcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGAc                                                                                                               <  1:1844233/150‑1 (MQ=255)
                                       gcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGAc                                                                                                               <  1:470929/150‑1 (MQ=255)
                                            aaCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGt                                                                                                          >  2:1559439/1‑150 (MQ=255)
                                              cGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTTTAGGCCTGATAATCGAAGCGCATCAGGCAAAGTTGCGTTTGTCAACCGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAg                                                                                                        <  2:1410611/150‑1 (MQ=18)
                                                  ttATCCGTAATGAAAACCCTTAAAATCCAATATATTGCAGTGAACGTGTAGGCCTGATAATCGTAGCGCTTCAGTCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCtt                                                                                                    <  2:160776/150‑1 (MQ=18)
                                                     tccgtgaTGAAAACCCTTCAAATCCAATAGCTTGCACTGAACTTGTAGGCCTGAGATGCGTATCTATTCAGTCAATGTTTATTTTCTTATTAGTTGAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGt                                                                                                 <  2:282412/144‑1 (MQ=16)
                                                      ccGGCATGAAAACCCTTCAAATCCAATAGTTTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGtt                                                                                                <  2:1756929/150‑1 (MQ=255)
                                                        ggCATGAAAACCCTTCATATCCAATAGATTGCAGTGAACGTGTAGTCCTGATAAGCGTATCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCa                                                                                              <  2:420100/150‑1 (MQ=255)
                                                          cATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATg                                                                                            >  2:638099/1‑150 (MQ=255)
                                                           aTGAAAACCCTTGAAATCCTATAGATTGATGTGTACGTTTAGTAATTATATTCGGAGCGCCTCATTTAATGTTGCTTTTGTTATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGc                                                                                           <  2:164184/150‑1 (MQ=12)
                                                              aataCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCTTAGCGATTCAGTCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCgg                                                                                        <  2:937286/147‑1 (MQ=21)
                                                              aaaaCGCTTCATTTCCAATAGATTGCCGTGAACGTGTAGGCCTGATATGGGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCgg                                                                                        <  2:1312320/150‑1 (MQ=21)
                                                              aaaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTATACAACGTAGGCTTTCTTCATGCCgg                                                                                        >  2:1615951/1‑150 (MQ=255)
                                                                aaCCCTTCAAATCCAATATATTGCAGTGTACGTGCAGTCAAGATAAGCGGATCTCATCATGAAATGTTTCTTTTTTTATAAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAt                                                                                      <  2:1299171/150‑1 (MQ=18)
                                                                aaCCCTTCAAATCCAATAGATTGAAGTGAACGTGTAGGCCTGATAATCGTAGCGCATCAGGAAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAt                                                                                      <  1:1818486/150‑1 (MQ=255)
                                                                  cccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgc                                                                                    >  1:1560977/1‑150 (MQ=255)
                                                                   ccTTCAAATCCAATAGATTGCAGTGAACGTGTAGTCCTGATAAGCGTATCTCATCAGGCAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcg                                                                                   <  1:1070095/150‑1 (MQ=255)
                                                                        gaaTCCAATAGATTGCACTGAACGTGTAGGCCTGATAAGCGTATCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTg                                                                              <  2:1520663/149‑1 (MQ=255)
                                                                              aaTAGATTGCAGTGAACGTGTAGGCCTTATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                        <  2:431285/150‑1 (MQ=255)
                                                                              aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCTTCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                        <  2:1727236/150‑1 (MQ=255)
                                                                              aaTAGATTGAAGTGAACGTGTAGGCCTGATAAGGGTAGCGAATGAGTCAATGTTGCTTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCAAGCCGGATGCGGCGTGAACGcc                                                                        <  2:1927585/150‑1 (MQ=21)
                                                                                     tGCAGTGAACGTGTAGGCCTGATAAGCGTATCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                 <  2:1643709/150‑1 (MQ=255)
                                                                                       ccgcgaccGTGTAGGCCTGATACGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                               <  1:1389955/143‑1 (MQ=21)
                                                                                         gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                             <  1:533923/150‑1 (MQ=255)
                                                                                          tGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCATg                                                            >  1:1370929/1‑150 (MQ=255)
                                                                                          tGAAAGTTTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTGCTTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATg                                                            <  2:1452713/150‑1 (MQ=255)
                                                                                               gtgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAc                                                       >  1:894471/1‑150 (MQ=255)
                                                                                               gtgtAGGCCTCATACGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAc                                                       <  1:706360/150‑1 (MQ=255)
                                                                                                     gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCa                                                 >  1:1736751/1‑150 (MQ=255)
                                                                                                      ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaa                                                <  1:782147/150‑1 (MQ=255)
                                                                                                      ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaa                                                <  1:1503361/150‑1 (MQ=255)
                                                                                                       cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGGCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaaa                                               <  2:1573055/150‑1 (MQ=255)
                                                                                                          aTAAGCGTAGCGCATCAGGCAATGTTGCGTTTTTCATAATTTTCAAATGGCTCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATcc                                            <  1:552870/150‑1 (MQ=255)
                                                                                                          aTAAGCGTAGCGCATCAGGAAATGTTGCTTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATcc                                            <  2:562041/150‑1 (MQ=255)
                                                                                                             aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAAt                                         <  2:453122/150‑1 (MQ=255)
                                                                                                             aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAATGAGTCATTTTAATATTGTATAAAACGTAGGCTTTGTACATGCCCGAACCGTCGTGATAGACATAAGCCGAAACAAAACCCTTcaaatcaaat                                         >  2:241190/1‑150 (MQ=12)
                                                                                                              gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTAAAAACCAATa                                        >  2:1932223/1‑150 (MQ=255)
                                                                                                                gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGGGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGa                                      >  1:1892675/1‑150 (MQ=255)
                                                                                                                gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACCACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATACGGCATGAAAAACCTTCAAATCAAATAGa                                      >  1:1742849/1‑150 (MQ=18)
                                                                                                                gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGa                                      >  1:458591/1‑150 (MQ=255)
                                                                                                                gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAGCCGTCAAATCCAATAGa                                      >  1:1711124/1‑150 (MQ=255)
                                                                                                                  aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATCCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGAtt                                    >  1:1765863/1‑150 (MQ=255)
                                                                                                                     tcATGAGGAAATGTTGTTTTTTTCTTCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCa                                 <  2:1933494/149‑1 (MQ=14)
                                                                                                                          aGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGaa                            >  1:1314885/1‑150 (MQ=255)
                                                                                                                           ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATTCAGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATAGCAGTGaaa                           >  1:1663218/1‑149 (MQ=255)
                                                                                                                            gCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACCTAGGCTTTGTTCATGCCGGATGAGTCGTGAAAGCCTTATACGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACg                          >  1:1160298/1‑150 (MQ=17)
                                                                                                                             cAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgt                         >  2:1975471/1‑150 (MQ=255)
                                                                                                                              aaTGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgtg                        >  2:1807591/1‑150 (MQ=255)
                                                                                                                                 gTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCATGAAAAACCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                     >  1:1246394/1‑150 (MQ=255)
                                                                                                                                   tGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGc                   >  2:837438/1‑150 (MQ=255)
                                                                                                                                    tcGTTTTTAATAAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                  <  2:1605163/149‑1 (MQ=255)
                                                                                                                                         tgaaatcaGATTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATa             <  2:627909/146‑1 (MQ=255)
                                                                                                                                         tGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATa             >  2:76609/1‑150 (MQ=255)
                                                                                                                                          gTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCTTGAACGCCATATCCGGCAGGAAAACCCTTAAAATCCAAAAGATATCAGTGAAAGTGTAGGCCTGATaa            >  1:304548/1‑150 (MQ=14)
                                                                                                                                             atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg         >  1:612399/1‑150 (MQ=39)
                                                                                                                                             atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg         >  1:1854249/1‑150 (MQ=39)
                                                                                                                                             atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg         >  2:667103/1‑150 (MQ=39)
                                                                                                                                             atcaGTTTCAAATGGCGCTGTAAAAGGCGTCAGTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAGGCGGCGTGAACGCCTGATCCGGCATGAAAACCCTTCAAATCCAATAGCTTGCAGTGAACGTGTAGGCCTGATAAgag         >  1:1783806/1‑148 (MQ=255)
                                                                                                                                              tcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTAATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt        >  2:238888/1‑150 (MQ=38)
                                                                                                                                                 gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATACGCGTAgc     >  2:1491720/1‑150 (MQ=25)
                                                                                                                                                 gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc     >  1:1047706/1‑150 (MQ=25)
                                                                                                                                                    tCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAGCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCa  >  1:1647853/1‑150 (MQ=18)
                                                                                                                                                     |                                                                                                                                                    
ATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCA  >  NC_000913/1552042‑1552339

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