Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,191 C→A 24.4% 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,1910CA24.4% 135.4 / 15.8 74intergenic (+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 (20/34);  new base A (18/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 4.95e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAT  >  NC_000913/1552042‑1552340
                                                                                                                                                     |                                                                                                                                                     
atatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                                                                                                                                       <  2:689488/150‑1 (MQ=31)
 tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGACAATGTTGCGTTTGTCATCAGTTTCa                                                                                                                                                      <  2:1761332/150‑1 (MQ=17)
 tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGCGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCa                                                                                                                                                      <  2:1514650/150‑1 (MQ=17)
          cAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCt                                                                                                                                             <  2:1528903/150‑1 (MQ=34)
           aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                            <  2:1560007/150‑1 (MQ=35)
                aGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaaag                                                                                                                                       <  2:820078/150‑2 (MQ=38)
                aGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaaaa                                                                                                                                       >  2:955079/1‑150 (MQ=39)
                 ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTAATCAGTTTAAAATGGCGCTGTAAAAg                                                                                                                                      >  1:1591487/1‑150 (MQ=255)
                 ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGAATGAAAACCCTTCAAATCCTATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCCTTTGTCTTCAGTTTCAAATGGCGCTGTAAAAg                                                                                                                                      <  2:1166836/150‑1 (MQ=255)
                  gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAgg                                                                                                                                     <  2:1257319/150‑1 (MQ=40)
                   cTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGc                                                                                                                                    <  1:1121174/150‑1 (MQ=255)
                    tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                                                                                                                   <  1:634827/150‑1 (MQ=255)
                     ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGt                                                                                                                                  >  1:1458162/1‑150 (MQ=255)
                         tCATGCCGGATGCGGCGTTAACGCCTTATCCGGCATGAAAACCCTTCAAATCAAATAGATTGCAGTGAACGTGTATGACTGATAATCGTAGCTCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAATGCGTCAtt                                                                                                                              >  2:1630287/1‑150 (MQ=14)
                          cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAttt                                                                                                                             <  2:1411089/150‑1 (MQ=255)
                             gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCa                                                                                                                          <  2:869445/150‑1 (MQ=255)
                              ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAGGCGTCATTTTAat                                                                                                                         >  2:1480778/1‑150 (MQ=255)
                               cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTATCGCATAAGGCAATGTTGCGTTAGTCATCAGTTTAAAATGGCGCTGTATAATGAGTCAttttatga                                                                                                                        >  2:308958/1‑146 (MQ=11)
                                  ctGCGGCGTGAACGCCTTATCTTTCATGAAAACCATTCAAATCCAATAGATTCCAGTGATCTTTTACGCCTGATCAGCGAAGCGCATCATGCTATGTTTTGTTATTCTTCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTg                                                                                                                     <  2:83782/149‑1 (MQ=255)
                                       gcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCTATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGAc                                                                                                                <  2:701225/150‑1 (MQ=255)
                                              cGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCTCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAg                                                                                                         <  1:1742172/150‑1 (MQ=255)
                                                ccTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGc                                                                                                       <  1:520326/150‑1 (MQ=255)
                                                 cTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAGGAGTCATTTTATTAATATATACACCGTCGGCt                                                                                                      >  2:1690726/1‑150 (MQ=11)
                                                    aTCCTGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTg                                                                                                   <  2:769791/150‑1 (MQ=255)
                                                     tCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGt                                                                                                  >  1:902586/1‑150 (MQ=255)
                                                      cctgCATGAAAACCCTTCAAATCCTATAGATTGAATTGAACGTGTAGTCATGATAAGCGTATCTATTCATTAAATGTTGTGTTTGTCATAAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGtt                                                                                                 <  2:486384/147‑1 (MQ=21)
                                                         gCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTATACAACGAAGGCTTTGTTCAt                                                                                              >  2:892618/1‑150 (MQ=18)
                                                         gCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTATAAGGCGACATTTTAATATTGTATACACCGACGGCTTTGTTCAt                                                                                              >  2:1146038/1‑150 (MQ=14)
                                                            tGAAAACCCTTCACATCCAATAGATTGCATTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGcc                                                                                           <  2:1491641/150‑1 (MQ=255)
                                                             gAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTATACAACGTAGGCTTAGTTCATGccc                                                                                          >  2:534816/1‑149 (MQ=18)
                                                              aaaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCgg                                                                                         <  1:955079/150‑1 (MQ=255)
                                                                  cccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGACATTTTAATATTGTAGACAAAGAAGGCTTTGTACATGCCGGAtcc                                                                                     >  2:1008329/1‑148 (MQ=18)
                                                                  cccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTATAAAGCGACATTTTCATATTGAAGACAAAGTAGGCTTTGATCATGcacgaaca                                                                                     >  1:1662922/1‑143 (MQ=14)
                                                                      tCAAATCCAATAGATTGCTGTGAACGTGTAGGCCTGATAAGCGTATGGCATGAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggcg                                                                                 <  1:1254224/150‑1 (MQ=255)
                                                                         aaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGa                                                                              >  1:1306098/1‑150 (MQ=255)
                                                                            ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGATTCAAATGGCGCGGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                           <  2:824563/150‑1 (MQ=255)
                                                                            acattataTTGCTTTGAACGTTTAGGACAGATAAGCGTAGCTCATCAGGCAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                           <  2:113975/143‑1 (MQ=21)
                                                                              aaTAGATTGCAGTGCACGTGTAGGACTGATAAGCGTAGCTCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                         <  2:474235/150‑1 (MQ=255)
                                                                              aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                         <  2:1748779/150‑1 (MQ=255)
                                                                               aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                        <  2:1458162/150‑1 (MQ=255)
                                                                                   aTTGCAGTGAACGTGTAGTCCTGATAATCGGAGCGCATCAGGCAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                    <  2:1004480/150‑1 (MQ=255)
                                                                                       cAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                <  1:459429/150‑1 (MQ=255)
                                                                                        aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                               <  2:1614995/150‑1 (MQ=255)
                                                                                          tGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATg                                                             >  1:801603/1‑150 (MQ=255)
                                                                                           gAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAAGAGTCAATTTAATATTTTATAAAACGTATGCTTTTTTAATGCAGTATGCGGCGTGAACGCCTTATACgtaatta                                                            >  2:383076/1‑144 (MQ=17)
                                                                                              agtgtAGGCATGATAAGGGTATCGCATGAGTCAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGTATGCGGCGTGAACGCCTTATCCGGCATGaaaa                                                         <  2:1389320/149‑1 (MQ=18)
                                                                                               gtgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAc                                                        >  1:1506719/1‑150 (MQ=255)
                                                                                               gtgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACACCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTAATCCGGCATGAAAAc                                                        >  2:25590/1‑150 (MQ=255)
                                                                                                 gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAACGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAATAccc                                                      >  2:1290133/1‑150 (MQ=255)
                                                                                                 gtACGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAACGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAccc                                                      <  1:1529458/150‑1 (MQ=255)
                                                                                                  tAGTCCTGATAAGCGTAGCGCATCAGGAAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCt                                                     <  2:732757/150‑1 (MQ=255)
                                                                                                    ggCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAAGCCTTATCCGGCATGAAAAcaatgc                                                   >  2:991677/1‑145 (MQ=255)
                                                                                                      ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTATACAAAGAAGGCTTTTTACATTCCGGATGCGTCGTGAAACACATAAACCGCAAGAAAACACTTCaa                                                 >  2:651380/1‑150 (MQ=16)
                                                                                                           tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCa                                            <  2:1306098/150‑1 (MQ=255)
                                                                                                              gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGAGTTATCCGGCATGAAAACCCTTCAAATCCAATa                                         >  1:213792/1‑150 (MQ=255)
                                                                                                              gCGTAGCGCATCAGGCAATCTTTCGTGTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATa                                         <  2:1102412/150‑1 (MQ=255)
                                                                                                                gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGGCATTTTAATATTGAATACAACGAAGGCTTTGTTAATGCCGGAAGACGCGTGAACGCATTATACGGCATGAAAACCGAGCACATCAAATAGa                                       >  1:1099088/1‑150 (MQ=21)
                                                                                                                    cgcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                   <  2:801603/150‑1 (MQ=255)
                                                                                                                     gcATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGAAGACAACGAAGGCTTATTACATTCCGAAAGACGCGTGAAAGACATATCAGGAAAGAAAAACATTAAAATACAAAAGATTGCa                                  >  1:460501/1‑150 (MQ=12)
                                                                                                                           ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGACATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATACGGAATGAAAACCCTTCAAATCCAATAGATTGCAGTTAAc                            >  1:1701640/1‑150 (MQ=17)
                                                                                                                            gCAATGTTGCGTTTGTCATCAGTTTCAGATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATTCCGGATGCGGCGAGAACGACTAATCCGGCATGAAAACCCTTCAAATCAAATAGAATGCAGTTAACg                           >  2:810891/1‑150 (MQ=17)
                                                                                                                            gCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTATAAAAAGTAGGCTTTGTTCATTACAGATGAGGAGTGAAAGCATTATCCTGAATGAATACACTTCAAATAAAATAGATTTCAGAGaaag                           >  1:269964/1‑148 (MQ=255)
                                                                                                                                tgtttgtttttTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTa                       <  2:1327046/143‑1 (MQ=255)
                                                                                                                                    gCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAAGCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                   >  2:719314/1‑150 (MQ=255)
                                                                                                                                     cGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTTAACGTGTAGGCCt                  >  2:1373555/1‑150 (MQ=255)
                                                                                                                                     cGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCAGGATGCGGCGTGATCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCTGTGAAGGTGTAGGCCt                  >  2:890540/1‑150 (MQ=255)
                                                                                                                                          gTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa             >  1:1319089/1‑150 (MQ=255)
                                                                                                                                            catcaGATTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc           >  1:5094/1‑150 (MQ=38)
                                                                                                                                             atcaGTTTCCAATGGCGCTGTAAAAGGCGTCATTTTCAGATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGTGGCGTGCACGCCTTATCCGGCATGAAAACCCTTCATATCCAATAGTTAGCAGTGAACTTTTAGGCCTGATAAGct          >  1:270732/1‑149 (MQ=11)
                                                                                                                                             atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg          >  1:716992/1‑150 (MQ=39)
                                                                                                                                               caGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa        >  1:448010/1‑150 (MQ=38)
                                                                                                                                                 gcttCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc      <  1:1055064/148‑1 (MQ=21)
                                                                                                                                                 gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTACAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc      >  1:1506699/1‑150 (MQ=25)
                                                                                                                                                    tCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCa   >  1:293005/1‑150 (MQ=33)
                                                                                                                                                    tCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCa   >  1:401261/1‑150 (MQ=33)
                                                                                                                                                     cAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACACTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt  >  1:1129699/1‑150 (MQ=17)
                                                                                                                                                     |                                                                                                                                                     
ATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAT  >  NC_000913/1552042‑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: