Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,026,538 G→A 9.0% intergenic (+13/+19) helD → / ← mgsA DNA helicase IV/methylglyoxal synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,026,5380GA 9.0% 261.6 / 11.0 89intergenic (+13/+19)helD/mgsADNA helicase IV/methylglyoxal synthase
Reads supporting (aligned to +/- strand):  ref base G (40/41);  new base A (4/4);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACC  >  NC_000913/1026389‑1026686
                                                                                                                                                     |                                                                                                                                                      
cTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGAGAGAAAACCGTAAGAAACAGGTGGCg                                                                                                                                                        >  1:1071547/1‑150 (MQ=255)
  gaagatCGGCGGGTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACACAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAATAAACAGGTGGCGtt                                                                                                                                                      <  1:975151/150‑1 (MQ=255)
  gaagaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGTTGGCGAGAAAACCGTAAGAAACAGGTGGCGtt                                                                                                                                                      >  2:862216/1‑150 (MQ=255)
   aagaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAATGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGttt                                                                                                                                                     >  1:819399/1‑150 (MQ=255)
     gaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGc                                                                                                                                                   <  2:1119851/150‑1 (MQ=255)
     gaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGc                                                                                                                                                   <  2:287149/150‑1 (MQ=255)
      aaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGcc                                                                                                                                                  <  2:1429277/150‑1 (MQ=255)
      aaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGcc                                                                                                                                                  <  2:518228/150‑1 (MQ=255)
      aaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGcc                                                                                                                                                  <  1:1233383/150‑1 (MQ=255)
           cggTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtg                                                                                                                                             <  1:1667944/150‑1 (MQ=255)
            ggTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAGTACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgt                                                                                                                                            <  2:867965/150‑1 (MQ=255)
                taTGTAGGTGGCGCTGAGCCGGGCAGGCCTTCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGTAAATTCTGGATGTGCCGGTGGCGAGATAACCGTAAGAAACAGGTGGCAAACGCCACCtgtttct                                                                                                                                        <  2:1660982/149‑5 (MQ=255)
                   gTACGTGGCGCTGACCCGGGCTCGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgtttcttac                                                                                                                                     <  1:380987/150‑8 (MQ=255)
                    tACGTGGCGCTGACCCGTGCACGCCAGCGCGTAGGGGCACTGTTTAACAAAGAGCATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAACCCGTAAGCAACAGGTGGCGTTTGCCACCTGTGCAATAtt                                                                                                                                    <  2:1147368/150‑1 (MQ=255)
                    tACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATAtt                                                                                                                                    <  2:816370/150‑1 (MQ=255)
                      cGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAGTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTAc                                                                                                                                  >  1:667703/1‑150 (MQ=255)
                      cGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTAc                                                                                                                                  >  2:721448/1‑150 (MQ=255)
                      cGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTAc                                                                                                                                  <  2:1239397/150‑1 (MQ=255)
                           ggcTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAg                                                                                                                             <  2:1334883/149‑1 (MQ=255)
                           cggtGACGCGGGCACGCCATCGTGTATGGGCACTGTTTAACAAAGAGAATCGCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAg                                                                                                                             <  2:116492/147‑1 (MQ=255)
                           cgcTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAg                                                                                                                             >  1:365092/1‑150 (MQ=255)
                           cgcTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAg                                                                                                                             <  1:889042/150‑1 (MQ=255)
                              tGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACg                                                                                                                          >  1:56090/1‑150 (MQ=255)
                                 cccGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAACCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTc                                                                                                                       >  1:1316324/1‑150 (MQ=255)
                                 cccGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTc                                                                                                                       >  1:661232/1‑150 (MQ=255)
                                  ccGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTcc                                                                                                                      <  1:134359/150‑1 (MQ=255)
                                   cGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAAAATTACTTCAGACGGTCcg                                                                                                                     >  2:200711/1‑150 (MQ=255)
                                     ggCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGAGGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCcgcg                                                                                                                   <  1:1668983/150‑1 (MQ=255)
                                            cATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCGGGATGTGCCGGTTGCGAGAAGACCGTAAGAAACAGGTGGCGATTGCCACCTGTGCAATATTACATCATACGGTCCGCGAGATAAc                                                                                                            >  2:313720/1‑150 (MQ=255)
                                               cGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCt                                                                                                         <  1:516851/150‑1 (MQ=255)
                                                gggTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTg                                                                                                        <  2:1378116/150‑1 (MQ=255)
                                                       ggCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑g                                                                                               >  2:764084/1‑150 (MQ=255)
                                                        gCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑gg                                                                                              <  2:1252856/150‑1 (MQ=255)
                                                         cACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑ggg                                                                                             >  1:1053358/1‑150 (MQ=255)
                                                          aCTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑gggg                                                                                            <  2:869614/150‑1 (MQ=255)
                                                             gTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATc                                                                                         >  1:1139461/1‑150 (MQ=255)
                                                               ttAAGAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAGATC‑ggaagagc                                                                                        <  2:301713/150‑7 (MQ=255)
                                                               ttAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAg                                                                                       <  2:1629748/150‑1 (MQ=255)
                                                                tAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACTGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGa                                                                                      >  2:702745/1‑150 (MQ=255)
                                                                 aaCAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGACGGTGGCGAGAAAACCGTAAGAACCAGGTGGCGTTTGCCACATGTGCCATATTACTTCAGACGCTCCGCGAGAAAACGCTGATC‑ATC‑GCGGAGCAGaa                                                                                     >  1:279924/1‑150 (MQ=255)
                                                                    aaaGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAatat                                                                                  >  1:426487/1‑150 (MQ=255)
                                                                     aaGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCAGCGAGATAACGCTGATA‑ATC‑GTGGATCAGAatata                                                                                 >  1:1084248/1‑149 (MQ=255)
                                                                      agagAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCTCGATATAACGCTGATA‑ATC‑GGGGCTCTGAATATCg                                                                                >  1:148503/1‑150 (MQ=255)
                                                                      agagAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCg                                                                                >  2:1677890/1‑150 (MQ=255)
                                                                       gagaATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGa                                                                               <  1:1426507/150‑1 (MQ=255)
                                                                       cagaATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGa                                                                               >  1:172908/2‑150 (MQ=255)
                                                                         gaATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGAcc                                                                             <  2:1478101/150‑1 (MQ=255)
                                                                         gaATCCCTCTACCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGAAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGAcc                                                                             <  2:1688547/150‑1 (MQ=255)
                                                                              ccTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGcgtc                                                                        >  2:808472/1‑150 (MQ=255)
                                                                               ctctCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGcgtcg                                                                       >  2:456280/1‑150 (MQ=255)
                                                                                 ctcCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGtt                                                                     <  2:1712658/150‑1 (MQ=255)
                                                                                   cccTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGCATATCGACCGCGTCGTTGa                                                                   >  1:1089035/1‑150 (MQ=255)
                                                                                   cccTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGa                                                                   <  1:772805/150‑1 (MQ=255)
                                                                                          tgGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATgcgg                                                            <  1:738385/150‑1 (MQ=255)
                                                                                           ggAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTGGAAATgcggc                                                           >  1:666834/1‑150 (MQ=255)
                                                                                              aaTACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGCTCAGAATATCGACCGCGTCGTTGAAATGCGGCGAc                                                        >  1:1159039/1‑150 (MQ=255)
                                                                                               aTACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTAAGACGGTCCGCGAGATTACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGATGAAAAGCGGCGACt                                                       >  2:635820/1‑150 (MQ=255)
                                                                                                 aCTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTgg                                                     <  1:759662/150‑1 (MQ=255)
                                                                                                   tGAAAAATCTGGATGTTCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAt                                                   <  2:1476686/150‑1 (MQ=255)
                                                                                                   tGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAt                                                   >  1:1192972/1‑150 (MQ=255)
                                                                                                       aaaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATg                                               <  2:172908/150‑1 (MQ=255)
                                                                                                       aaaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATg                                               <  2:207735/150‑1 (MQ=255)
                                                                                                       aaaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATg                                               >  2:691435/1‑150 (MQ=255)
                                                                                                        aaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGa                                              <  1:1430034/150‑1 (MQ=255)
                                                                                                         aTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGaa                                             >  2:855050/1‑150 (MQ=255)
                                                                                                                 gtgCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCg                                     >  1:592746/1‑150 (MQ=255)
                                                                                                                 gtgCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCg                                     <  2:1316324/150‑1 (MQ=255)
                                                                                                                  tgCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGt                                    >  1:1384972/1‑150 (MQ=255)
                                                                                                                         ggCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACg                             >  2:1541534/1‑150 (MQ=255)
                                                                                                                           cGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATCAGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGt                            >  1:64102/1‑150 (MQ=255)
                                                                                                                           cGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACAGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGtt                           >  1:1228818/1‑150 (MQ=255)
                                                                                                                            gagaAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGAGCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtg                          >  1:333340/1‑150 (MQ=255)
                                                                                                                              gaAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGACCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGAGGAGACTGGATTATGAAGTCTGCAGTTGCCACCTtgct                        >  2:785513/1‑148 (MQ=255)
                                                                                                                               aaaaCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtggtg                       <  2:64304/150‑1 (MQ=255)
                                                                                                                               aaaaCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtggtg                       <  1:862216/150‑1 (MQ=255)
                                                                                                                               aaaaCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtggtg                       <  1:1709802/150‑1 (MQ=255)
                                                                                                                                  aCCGTAAGAAACAGGTGGAGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCg                    <  2:968405/150‑1 (MQ=255)
                                                                                                                                     gTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGAcc                 <  2:770107/150‑1 (MQ=255)
                                                                                                                                     gTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGAcc                 >  1:1139541/1‑150 (MQ=255)
                                                                                                                                     atattgcaCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGAcc                 >  1:399042/8‑150 (MQ=255)
                                                                                                                                         gAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGCCCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCAGTTGCCACGTTGGTGGCGACCGGaa             >  1:1603639/1‑150 (MQ=255)
                                                                                                                                         gAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGaa             <  2:1192972/150‑1 (MQ=255)
                                                                                                                                           caCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATg           >  1:1605825/2‑150 (MQ=255)
                                                                                                                                            aCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGt          >  2:381995/1‑150 (MQ=255)
                                                                                                                                            aCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGt          >  1:1326324/1‑150 (MQ=255)
                                                                                                                                             cAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGtt         >  2:565751/1‑150 (MQ=255)
                                                                                                                                              cggTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGGCCGCGAGATAACGCTGATA‑AGC‑GGGGAGCAGAATATCGCCCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGGCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTc        <  1:1054978/149‑1 (MQ=255)
                                                                                                                                               ggTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTcc       >  1:929887/1‑150 (MQ=255)
                                                                                                                                                  ggCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATa    <  1:721448/150‑1 (MQ=255)
                                                                                                                                                    cGTTTCCCACCTGTGCAATATTACTTCAGACCGTCCGCGAGATAACGCTGAAA‑CTC‑GCGGACCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATAcc  <  1:1069076/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACC  >  NC_000913/1026389‑1026686

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: ‑