Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,026,541 T→C 9.2% intergenic (+16/+16) 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,5410TC 9.2% 235.1 / 11.8 86intergenic (+16/+16)helD/mgsADNA helicase IV/methylglyoxal synthase
Reads supporting (aligned to +/- strand):  ref base T (38/40);  new base C (4/4);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

AAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCG  >  NC_000913/1026392‑1026690
                                                                                                                                                     |                                                                                                                                                       
aagaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAATGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGttt                                                                                                                                                         >  1:819399/1‑150 (MQ=255)
  gaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGc                                                                                                                                                       <  2:1119851/150‑1 (MQ=255)
  gaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGc                                                                                                                                                       <  2:287149/150‑1 (MQ=255)
   aaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGcc                                                                                                                                                      <  2:518228/150‑1 (MQ=255)
   aaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGcc                                                                                                                                                      <  2:1429277/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:889042/150‑1 (MQ=255)
                        cgcTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAg                                                                                                                                 >  1:365092/1‑150 (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:207735/150‑1 (MQ=255)
                                                                                                    aaaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATg                                                   <  2:172908/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                           <  1:1709802/150‑1 (MQ=255)
                                                                                                                            aaaaCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtggtg                           <  2:64304/150‑1 (MQ=255)
                                                                                                                            aaaaCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtggtg                           <  1:862216/150‑1 (MQ=255)
                                                                                                                               aCCGTAAGAAACAGGTGGAGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCg                        <  2:968405/150‑1 (MQ=255)
                                                                                                                                  gTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGAcc                     >  1:1139541/1‑150 (MQ=255)
                                                                                                                                  gTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGAcc                     <  2:770107/150‑1 (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              >  1:1326324/1‑150 (MQ=255)
                                                                                                                                         aCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGt              >  2:381995/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)
                                                                                                                                                     tGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCTGAATATCGACCGCGTCGTTGAAATGCGGCGAGTGGATTATGATGTCTTCCGTTGCCACGTTGGTGGCGTCCGGTATTTTACATACCGTCg  >  2:121638/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATA‑ATC‑GGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCG  >  NC_000913/1026392‑1026690

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