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

CTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACC  >  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)
                                                       ggCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCg                                                                                               >  2:764084/1‑150 (MQ=255)
                                                        gCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCgg                                                                                              <  2:1252856/150‑1 (MQ=255)
                                                         cACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCggg                                                                                             >  1:1053358/1‑150 (MQ=255)
                                                          aCTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCgggg                                                                                            <  2:869614/150‑1 (MQ=255)
                                                             gTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATc                                                                                         >  1:1139461/1‑150 (MQ=255)
                                                               ttAAGAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAGATCggaagagc                                                                                        <  2:301713/150‑7 (MQ=255)
                                                               ttAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAg                                                                                       <  2:1629748/150‑1 (MQ=255)
                                                                tAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACTGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGa                                                                                      >  2:702745/1‑150 (MQ=255)
                                                                 aaCAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGACGGTGGCGAGAAAACCGTAAGAACCAGGTGGCGTTTGCCACATGTGCCATATTACTTCAGACGCTCCGCGAGAAAACGCTGATCATCGCGGAGCAGaa                                                                                     >  1:279924/1‑150 (MQ=255)
                                                                    aaaGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAatat                                                                                  >  1:426487/1‑150 (MQ=255)
                                                                     aaGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCAGCGAGATAACGCTGATAATCGTGGATCAGAatata                                                                                 >  1:1084248/1‑149 (MQ=255)
                                                                      agagAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCTCGATATAACGCTGATAATCGGGGCTCTGAATATCg                                                                                >  1:148503/1‑150 (MQ=255)
                                                                      agagAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCg                                                                                >  2:1677890/1‑150 (MQ=255)
                                                                       gagaATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGa                                                                               <  1:1426507/150‑1 (MQ=255)
                                                                       cagaATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGa                                                                               >  1:172908/2‑150 (MQ=255)
                                                                         gaATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGAcc                                                                             <  2:1478101/150‑1 (MQ=255)
                                                                         gaATCCCTCTACCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGAAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGAcc                                                                             <  2:1688547/150‑1 (MQ=255)
                                                                              ccTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGcgtc                                                                        >  2:808472/1‑150 (MQ=255)
                                                                               ctctCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGcgtcg                                                                       >  2:456280/1‑150 (MQ=255)
                                                                                 ctcCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGtt                                                                     <  2:1712658/150‑1 (MQ=255)
                                                                                   cccTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGCATATCGACCGCGTCGTTGa                                                                   >  1:1089035/1‑150 (MQ=255)
                                                                                   cccTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGa                                                                   <  1:772805/150‑1 (MQ=255)
                                                                                          tgGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATgcgg                                                            <  1:738385/150‑1 (MQ=255)
                                                                                           ggAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTGGAAATgcggc                                                           >  1:666834/1‑150 (MQ=255)
                                                                                              aaTACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGCTCAGAATATCGACCGCGTCGTTGAAATGCGGCGAc                                                        >  1:1159039/1‑150 (MQ=255)
                                                                                               aTACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTAAGACGGTCCGCGAGATTACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGATGAAAAGCGGCGACt                                                       >  2:635820/1‑150 (MQ=255)
                                                                                                 aCTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTgg                                                     <  1:759662/150‑1 (MQ=255)
                                                                                                   tGAAAAATCTGGATGTTCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAt                                                   <  2:1476686/150‑1 (MQ=255)
                                                                                                   tGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAt                                                   >  1:1192972/1‑150 (MQ=255)
                                                                                                       aaaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATg                                               <  2:172908/150‑1 (MQ=255)
                                                                                                       aaaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATg                                               <  2:207735/150‑1 (MQ=255)
                                                                                                       aaaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATg                                               >  2:691435/1‑150 (MQ=255)
                                                                                                        aaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGa                                              <  1:1430034/150‑1 (MQ=255)
                                                                                                         aTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGaa                                             >  2:855050/1‑150 (MQ=255)
                                                                                                                 gtgCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCg                                     >  1:592746/1‑150 (MQ=255)
                                                                                                                 gtgCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCg                                     <  2:1316324/150‑1 (MQ=255)
                                                                                                                  tgCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGt                                    >  1:1384972/1‑150 (MQ=255)
                                                                                                                         ggCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACg                             >  2:1541534/1‑150 (MQ=255)
                                                                                                                           cGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCAGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGt                            >  1:64102/1‑150 (MQ=255)
                                                                                                                           cGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACAGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGtt                           >  1:1228818/1‑150 (MQ=255)
                                                                                                                            gagaAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGAGCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtg                          >  1:333340/1‑150 (MQ=255)
                                                                                                                              gaAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGACCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGAGGAGACTGGATTATGAAGTCTGCAGTTGCCACCTtgct                        >  2:785513/1‑148 (MQ=255)
                                                                                                                               aaaaCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtggtg                       <  2:64304/150‑1 (MQ=255)
                                                                                                                               aaaaCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtggtg                       <  1:862216/150‑1 (MQ=255)
                                                                                                                               aaaaCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtggtg                       <  1:1709802/150‑1 (MQ=255)
                                                                                                                                  aCCGTAAGAAACAGGTGGAGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCg                    <  2:968405/150‑1 (MQ=255)
                                                                                                                                     gTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGAcc                 <  2:770107/150‑1 (MQ=255)
                                                                                                                                     gTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGAcc                 >  1:1139541/1‑150 (MQ=255)
                                                                                                                                     atattgcaCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGAcc                 >  1:399042/8‑150 (MQ=255)
                                                                                                                                         gAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGCCCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCAGTTGCCACGTTGGTGGCGACCGGaa             >  1:1603639/1‑150 (MQ=255)
                                                                                                                                         gAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGaa             <  2:1192972/150‑1 (MQ=255)
                                                                                                                                           caCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATg           >  1:1605825/2‑150 (MQ=255)
                                                                                                                                            aCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGt          >  2:381995/1‑150 (MQ=255)
                                                                                                                                            aCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGt          >  1:1326324/1‑150 (MQ=255)
                                                                                                                                             cAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGtt         >  2:565751/1‑150 (MQ=255)
                                                                                                                                              cggTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGGCCGCGAGATAACGCTGATAAGCGGGGAGCAGAATATCGCCCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGGCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTc        <  1:1054978/149‑1 (MQ=255)
                                                                                                                                               ggTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTcc       >  1:929887/1‑150 (MQ=255)
                                                                                                                                                  ggCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATa    <  1:721448/150‑1 (MQ=255)
                                                                                                                                                    cGTTTCCCACCTGTGCAATATTACTTCAGACCGTCCGCGAGATAACGCTGAAACTCGCGGACCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATAcc  <  1:1069076/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CTGAAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACC  >  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: