Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,191 C→A 25.7% 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,1910CA25.7% 129.7 / 24.1 81intergenic (+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 (31/28);  new base A (18/4);  total (49/32)
Fisher's exact test for biased strand distribution p-value = 2.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAT  >  NC_000913/1552043‑1552340
                                                                                                                                                    |                                                                                                                                                     
tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCa                                                                                                                                                      <  2:1297235/150‑1 (MQ=32)
     tAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAAGAAAACCCTTAAAATCAAATAGATTGCAGTGAACGATAAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAAttg                                                                                                                                                 >  2:607508/1‑148 (MQ=11)
        aCAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcgc                                                                                                                                              <  2:1705638/150‑1 (MQ=34)
          aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                            <  1:212050/150‑1 (MQ=35)
          aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGGAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                            <  2:1894432/150‑1 (MQ=25)
           aCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                                                                                                                           <  2:446599/150‑1 (MQ=35)
             gTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaa                                                                                                                                         <  2:1805298/150‑1 (MQ=38)
                 gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTATGCCTGATAAGCGTAGCGCATAAGGCAATGTTGCCTTTGTCATCAGTTAAAAATGGCGCTGTaaaatc                                                                                                                                     >  2:328994/1‑148 (MQ=11)
                           tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGAGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                                                                                                           <  2:429/150‑1 (MQ=255)
                             ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                                                                                                         <  1:892503/150‑1 (MQ=255)
                             ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                                                                                                         <  2:1558135/150‑1 (MQ=255)
                               gttttcggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCATTATATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                                                                                                       <  1:110501/145‑1 (MQ=255)
                                gATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATAtt                                                                                                                      <  1:1024426/150‑1 (MQ=255)
                                  tGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                                                                                                                    >  1:953529/1‑150 (MQ=255)
                                    cggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAGGCGTCATTTTCATATTgaat                                                                                                                  >  2:615342/1‑147 (MQ=255)
                                      gcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGTCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGCGACATTTTAATATTGAAGAc                                                                                                                >  1:757675/1‑150 (MQ=14)
                                         tGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAc                                                                                                             <  2:953529/150‑1 (MQ=255)
                                          gAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACg                                                                                                            <  1:901246/150‑1 (MQ=255)
                                          gAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGCCGCTGTAAAAAGAGACATTTTAATATTGTATAAAAct                                                                                                            >  2:1797359/1‑149 (MQ=11)
                                             cGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAATGAGTCATTTTAATATTGAATACACAGTAg                                                                                                         >  2:399003/1‑150 (MQ=11)
                                              gCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAgg                                                                                                        <  1:26089/150‑1 (MQ=255)
                                               ccTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGc                                                                                                       <  2:596644/150‑1 (MQ=255)
                                               ccTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGc                                                                                                       >  2:645042/1‑150 (MQ=255)
                                                     ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGtt                                                                                                 <  2:1059984/150‑1 (MQ=255)
                                                     ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTAGACAACGTAGGCtttcta                                                                                                 >  1:1556221/1‑147 (MQ=255)
                                                     ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAAGCGACATTTTCATATTGTAGACAACGTAGGCTTTGtt                                                                                                 >  2:1879339/1‑150 (MQ=255)
                                                      cgtaaTGAAAACCCTACTAATCCTATAGATTGCTTTGAACGTGTAGTCCTTATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTc                                                                                                <  1:46873/146‑1 (MQ=21)
                                                         cATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAGGCGACATTTTAATATTGTATACAACGTAGGCTTTGTTCAtc                                                                                             >  2:1756309/1‑149 (MQ=17)
                                                         cATGAAAACCCTTCAAATCAAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGTCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATg                                                                                             <  2:1553248/150‑1 (MQ=255)
                                                           tGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGGCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGcc                                                                                           <  1:509709/150‑1 (MQ=255)
                                                            gAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCg                                                                                          <  1:802659/150‑1 (MQ=255)
                                                             aaaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGACATTTTAATATTGAATACAACGTCGGCTTTATACATACCgc                                                                                         >  2:1636650/1‑149 (MQ=17)
                                                              aaaCCCTTCAAATCCAATATATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCATTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGa                                                                                        <  2:750472/150‑1 (MQ=255)
                                                              aaaCCATTCAAATCCTATAGATTGCATTGAACGTGTATGCCTTATAATCTTAGCGATTCAGGAAATGTTGAGTTTGTCATAAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGa                                                                                        <  2:597317/150‑1 (MQ=21)
                                                                     tCAAATCCAATAGATTGCAGTGAACGTGGAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggcg                                                                                 <  1:1843510/150‑1 (MQ=255)
                                                                        aaTCCAATAGATTGCATTGCAAGTTTAGTACTTATAATCGTAGCTCATCAGTAAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGa                                                                              <  2:1216692/150‑1 (MQ=21)
                                                                          tCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                            <  2:1451276/150‑1 (MQ=255)
                                                                           ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGAAAACGTAGTCATTGTTCATGCCGGATGCGGCGTGAAcc                                                                           >  2:1905984/1‑149 (MQ=255)
                                                                            cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                          >  1:464444/1‑150 (MQ=255)
                                                                             aaTATATTGCAGTGAACGTCTAGTACTGATAAGCGTAGCTCATCAGGAAATGTTGCTTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                         <  2:1323301/150‑1 (MQ=21)
                                                                               tAGATTGCAGTGAACGTGTAGTCCTCATAAGCGTAGCGCATCATGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                       <  1:1005014/150‑1 (MQ=255)
                                                                                        gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCTTGAACGCCTTATCCGGCAt                                                              >  2:1386279/1‑150 (MQ=255)
                                                                                        gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCAGGATGCGGCGTGAAAGCCTTATCCGGCAt                                                              >  1:582086/1‑150 (MQ=255)
                                                                                                  aGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTACAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTATACAAAGAAGGCTTTATTCATGCCGGATGCGGCTTTATAGAATTATACGGAATGAAAACCCtt                                                    >  1:357627/1‑150 (MQ=21)
                                                                                                    gCCTGATAAGCGTAGCTCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCTCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCGTGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCa                                                  <  2:1518857/150‑1 (MQ=255)
                                                                                                     ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaa                                                 >  1:1430309/1‑150 (MQ=255)
                                                                                                     ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCATGAAAACCCTTCaa                                                 >  1:1571011/1‑150 (MQ=255)
                                                                                                      cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaaa                                                >  2:190866/1‑150 (MQ=255)
                                                                                                            aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCGTTATCCGCCATGAAAACCCTTCTACTCCAAt                                          >  2:480174/1‑150 (MQ=255)
                                                                                                             ccGTAGCGCATCAGGCAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATa                                         <  2:912006/149‑1 (MQ=255)
                                                                                                              cGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGAAGACAAAGAAGGCTTTGTACATGCAGGATCCGGATTTAAAGGCATATACGGAATGAAATCCATACAAAACAAATAg                                        >  2:1203960/1‑150 (MQ=16)
                                                                                                                tAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATACGGCATGAAAACCCTTCAAATCCAATAGAt                                      >  1:1782551/1‑150 (MQ=255)
                                                                                                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATACGGCATCAAAACCCTTCAAATCCAATAGAtt                                     >  1:1794010/1‑150 (MQ=255)
                                                                                                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTAGACAACGAAGGCTTTGTTCATTCCGGATGCGGCTTGAACGCCTTATCCGGAATGAAAACCCTTCAAATCCAAAAGAtt                                     >  1:838761/1‑150 (MQ=18)
                                                                                                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTATACAAAGTAGGCTTTTTTAATGCCCGATCCGGCGTGAAAGACATATACGGCATGAAAACACATCAAATACAATAGAtt                                     >  2:495052/1‑150 (MQ=18)
                                                                                                                  gcgcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTg                                    >  1:1541702/1‑150 (MQ=255)
                                                                                                                  gcgcATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTTAACGACTTATCCGGCATGAACACCCTTCAAATACAATAGATTg                                    >  2:1495510/1‑150 (MQ=17)
                                                                                                                   cgcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCGTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                   >  1:1496468/1‑150 (MQ=255)
                                                                                                                   cgcATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAGGAGACATTTTAATATTGAATAAAAAGAAGGCTTAGTTCATGCAGAATGCGGCGTGAACGCAGTTTCAGGCATGAAAACAATTCAAATCAAATAGATTga                                   >  2:801094/1‑149 (MQ=16)
                                                                                                                   cgcATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAAGAGACATTTTAATATTGTTTACAACGAAGGCTTTGTACATGCAGTAAGCTCCGTAGACGCAATATGAGGCAAGAAAACCATTCAAATCAAATAGAATgg                                   >  2:1029283/1‑149 (MQ=12)
                                                                                                                      aTCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGt                                >  1:1435030/1‑150 (MQ=255)
                                                                                                                       tCAGGCAATGTTGGGTTTTGTATCAGTTGCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTg                               <  2:1773022/150‑1 (MQ=14)
                                                                                                                         aGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAAGCGGCGTGAACGACATATCCGGCATGAACACCCTTCAAATCCAATAGATTGCAGTGaa                             >  2:186073/1‑150 (MQ=255)
                                                                                                                         aGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTTTAAAACGTAGGCTTTGTTCATGCCGGATGCGGCGTTAAAGACATATCCGGCATGAAAGCAATTCAAATCCAATAGATTGCAGTGaa                             >  2:492092/1‑150 (MQ=17)
                                                                                                                         aGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAGGGGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTAATGCCGGATGCGGCGTGAAAGCCATATGCGGCAAGAAAACCATTAAAATAAAATAGATTGCAGTGaa                             >  2:758873/1‑150 (MQ=17)
                                                                                                                          ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCAAGAAAACCCTTCAAATCCATTACATTCCATTGAAc                            >  1:342998/1‑150 (MQ=17)
                                                                                                                          ggCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAATGAGGCATTTTAATATTATAGACAAAGAAGGAATTGTACATTAAGAATGATTAGTGAACGCCTTATCCGGGATGAACCCCATTAAAATACAAAAGATTGCAGTGAAc                            >  2:1249311/1‑150 (MQ=12)
                                                                                                                                  tGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGACATTTTATTATTGAAGACACCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCGGAAAACGGAATGTACACCCTTCAAATGCAATAGATTGCAGTGAACGTGTATGc                    >  2:284684/1‑150 (MQ=14)
                                                                                                                                   gCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                   >  1:1853721/1‑150 (MQ=255)
                                                                                                                                     ttgtaggatcaGCTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTg                 <  2:927301/143‑1 (MQ=255)
                                                                                                                                          tcatcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACTTGTAGGCCTGATAAg            >  1:1578188/1‑150 (MQ=255)
                                                                                                                                           catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc           >  2:198346/1‑150 (MQ=40)
                                                                                                                                           catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc           >  2:861323/1‑150 (MQ=40)
                                                                                                                                            atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg          >  2:496546/1‑150 (MQ=39)
                                                                                                                                            atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg          >  1:778738/1‑150 (MQ=39)
                                                                                                                                            atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAACCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg          >  2:397547/1‑150 (MQ=255)
                                                                                                                                             tcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTGGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCATGAAAACCCTTCAAATCAAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt         >  2:46507/1‑150 (MQ=255)
                                                                                                                                             tcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt         >  1:56303/1‑150 (MQ=38)
                                                                                                                                              caGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa        >  2:367513/1‑150 (MQ=38)
                                                                                                                                              aaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa        <  2:21579/149‑1 (MQ=37)
                                                                                                                                                gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCATGCTAAGCGTAgc      >  1:1213808/1‑150 (MQ=25)
                                                                                                                                                    cAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt  >  2:1864862/1‑150 (MQ=32)
                                                                                                                                                    |                                                                                                                                                     
TATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAT  >  NC_000913/1552043‑1552340

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: