Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,945,985 C→T 5.2% intergenic (‑143/+96) tcdA ← / ← mltA tRNA threonylcarbamoyladenosine dehydratase/membrane‑bound lytic murein transglycosylase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,945,9850CT 5.2% 307.1 / 3.9 96intergenic (‑143/+96)tcdA/mltAtRNA threonylcarbamoyladenosine dehydratase/membrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base C (42/48);  new base T (3/2);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCACCCAGA  >  NC_000913/2945840‑2946130
                                                                                                                                                 |                                                                                                                                                 
cATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                               <  1:1416188/150‑1 (MQ=37)
  tAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                             >  1:357658/1‑150 (MQ=255)
    gCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                           >  1:404628/1‑150 (MQ=255)
    gCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                           <  1:678275/150‑1 (MQ=255)
     cACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                                          <  1:405305/150‑3 (MQ=37)
      aCCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                                                                                                         <  2:361889/150‑1 (MQ=255)
           aGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaa                                                                                                                                    <  2:44419/150‑1 (MQ=255)
            gAATTAAAAAGGGCGAGCTTAGCCCTCACCCTAAACGCATAAATGCCTGATCCGCTACGCTTATCAGGCCTACATAACACCTCAAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaagaac                                                                                                                                   >  2:1490343/1‑150 (MQ=21)
             aaTTAAAAAGGGCGAGGTTATCACTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTAAATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcg                                                                                                                                  >  2:1446142/1‑150 (MQ=255)
                tAAAAAGGGCGAGGTTATCCCTCACCCAAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACgcgc                                                                                                                               <  2:819403/150‑1 (MQ=255)
                 aaaaaGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCa                                                                                                                              <  2:804800/150‑1 (MQ=255)
                 aaaaaGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCa                                                                                                                              >  2:799093/1‑150 (MQ=255)
                       ggCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGt                                                                                                                        <  1:1569670/150‑1 (MQ=255)
                       ggCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGt                                                                                                                        >  2:435716/1‑150 (MQ=255)
                       ggCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGt                                                                                                                        >  1:79827/1‑150 (MQ=255)
                           aGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCaat                                                                                                                    >  2:252937/1‑150 (MQ=255)
                            ggTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCaata                                                                                                                   <  2:1135884/150‑1 (MQ=255)
                                  cccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAgg                                                                                                             <  2:1583250/150‑1 (MQ=255)
                                  cccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAgg                                                                                                             >  2:956567/1‑150 (MQ=255)
                                  cccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAgg                                                                                                             <  2:856711/150‑1 (MQ=255)
                                   ccTCCCCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGa                                                                                                            >  2:1119515/1‑150 (MQ=255)
                                   ccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCa                                                                                                               <  1:203257/147‑1 (MQ=255)
                                   ccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCa                                                                                                               >  2:203257/1‑147 (MQ=255)
                                   ccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGa                                                                                                            <  1:1554474/150‑1 (MQ=255)
                                   ccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGa                                                                                                            <  2:293323/150‑1 (MQ=255)
                                   ccTCACCCTAAACGCATAAATGCCTGACGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGa                                                                                                            <  2:1296478/150‑1 (MQ=255)
                                        cccTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGa                                                                                                       <  1:546186/150‑1 (MQ=255)
                                        cccTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGa                                                                                                       >  1:1463767/1‑150 (MQ=255)
                                        cccTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGa                                                                                                       >  1:225709/1‑150 (MQ=255)
                                          cTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAgg                                                                                                     <  1:162358/150‑1 (MQ=255)
                                               cGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAcc                                                                                                <  2:660758/150‑1 (MQ=255)
                                                 cATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAcccc                                                                                              <  1:1249499/150‑1 (MQ=255)
                                                    aaaTGCCTGATGCGCTACGCTTATCAGGCCTACATATCTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCa                                                                                           >  2:211264/1‑150 (MQ=255)
                                                     aaTGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAc                                                                                          <  1:956567/150‑1 (MQ=255)
                                                           tGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTc                                                                                    >  1:909828/1‑150 (MQ=255)
                                                            gATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCa                                                                                   <  2:623249/150‑1 (MQ=255)
                                                               gcgcTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCatat                                                                                >  1:962455/1‑150 (MQ=255)
                                                                cgcTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATc                                                                               <  1:252937/150‑1 (MQ=255)
                                                                cgcTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATc                                                                               >  2:622504/1‑150 (MQ=255)
                                                                cgcTACGCTTATCAGGCCTACAAAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATc                                                                               <  1:1247668/150‑1 (MQ=255)
                                                                 gcTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATcc                                                                              >  2:1464340/1‑150 (MQ=255)
                                                                 gcTACGATTAACAGGCCTACATAACTACTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGAAAAGGCCTTCACGCCGCATCCGTCATGAACAACGCGAACTAAGACAATAATCACGATGTGAGGTTACACCTCACCCTGTCATATcc                                                                              >  2:1391029/1‑150 (MQ=255)
                                                                   tACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGt                                                                            <  2:1411223/150‑1 (MQ=255)
                                                                         tATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAAc                                                                      >  1:1513577/1‑150 (MQ=255)
                                                                          aTCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACg                                                                     <  2:404628/150‑1 (MQ=255)
                                                                            cAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGc                                                                   >  2:1272584/1‑150 (MQ=255)
                                                                            cAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTACCCCTCACCCTGTCATATCCGTAAAAACGGCa                                                                  <  2:357658/150‑1 (MQ=255)
                                                                             aGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCa                                                                  <  1:22327/150‑1 (MQ=255)
                                                                             aGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCa                                                                  <  1:1272584/150‑1 (MQ=255)
                                                                              ggCCTACATAACTCCTGCAATCTTTTGGATTTGCATGAGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACACCGCGCACGTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCAt                                                                 <  2:1427848/150‑1 (MQ=255)
                                                                              ggCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCAt                                                                 >  2:1617049/1‑150 (MQ=255)
                                                                              ggCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCAt                                                                 >  2:877133/1‑150 (MQ=255)
                                                                                  tACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAg                                                             <  2:1124783/150‑1 (MQ=255)
                                                                                    cATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGaa                                                           >  1:1105755/1‑150 (MQ=255)
                                                                                      tAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAata                                                         <  1:1464340/150‑1 (MQ=255)
                                                                                         cTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATca                                                      >  2:11451/1‑150 (MQ=255)
                                                                                         cTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATca                                                      <  2:553409/150‑1 (MQ=255)
                                                                                            cTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACAt                                                   >  1:739245/1‑150 (MQ=255)
                                                                                                 atatTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGcc                                              <  2:709641/150‑1 (MQ=255)
                                                                                                 atatTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGcc                                              <  2:1105755/150‑1 (MQ=255)
                                                                                                    ttttGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCt                                           >  1:383191/1‑150 (MQ=255)
                                                                                                     tttGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTa                                          <  1:706632/150‑1 (MQ=255)
                                                                                                     tttGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTa                                          >  1:341675/1‑150 (MQ=255)
                                                                                                        ggATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAg                                       <  2:348643/150‑1 (MQ=255)
                                                                                                        ggATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAg                                       <  2:342674/150‑1 (MQ=255)
                                                                                                         gATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACACTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGa                                      >  2:1149949/1‑150 (MQ=255)
                                                                                                           tttGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACg                                    <  1:1296257/150‑1 (MQ=255)
                                                                                                            ttGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGt                                   <  2:1463767/150‑1 (MQ=255)
                                                                                                             tGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGtt                                  >  1:321452/1‑150 (MQ=255)
                                                                                                               cATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTAc                                >  2:78150/1‑150 (MQ=255)
                                                                                                                aTGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTAcc                               <  1:11451/150‑1 (MQ=255)
                                                                                                                 tGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGATAAAGACggcgtca                              >  1:1525928/1‑144 (MQ=255)
                                                                                                                   aTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTgc                            >  2:367514/1‑150 (MQ=255)
                                                                                                                    tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTgcg                           <  1:374437/150‑1 (MQ=255)
                                                                                                                     gttcatgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTgcgc                          >  1:173463/7‑150 (MQ=38)
                                                                                                                      tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGcc                         <  1:1446142/150‑1 (MQ=255)
                                                                                                                       gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGccc                        >  2:1428603/1‑150 (MQ=255)
                                                                                                                         atgCCGGATGCGGCGTGAACGCCTTATCTGGCATGAACACCGCGCACTTAGTCAATAATCAGGATGTGAGGTTCCCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCCCATCAGCCGCTAAAGACGTTCCCTGCGCCCgg                      >  1:1305037/3‑150 (MQ=255)
                                                                                                                         aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCgg                      >  1:800500/1‑150 (MQ=255)
                                                                                                                         aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCgg                      <  1:78150/150‑1 (MQ=255)
                                                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCggg                     <  1:622504/150‑1 (MQ=255)
                                                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCgggg                    >  2:89023/1‑150 (MQ=255)
                                                                                                                            ccGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGc                   >  1:1055198/1‑150 (MQ=255)
                                                                                                                            ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGc                   >  2:525979/1‑150 (MQ=255)
                                                                                                                            ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGc                   <  2:791087/150‑1 (MQ=255)
                                                                                                                              ggATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCgg                 <  2:739245/150‑1 (MQ=255)
                                                                                                                              ggATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCgg                 <  2:856844/150‑1 (MQ=255)
                                                                                                                                   aGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTc            <  1:1521219/150‑1 (MQ=255)
                                                                                                                                   aGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTc            >  1:1306517/1‑150 (MQ=255)
                                                                                                                                   aGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTc            <  1:1617049/150‑1 (MQ=255)
                                                                                                                                   aGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTc            >  1:277266/1‑150 (MQ=255)
                                                                                                                                    ggCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCa           >  1:431080/1‑150 (MQ=255)
                                                                                                                                      cGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGc         <  1:1605684/150‑1 (MQ=255)
                                                                                                                                       gTTCATGCCGGATGCGGCGTGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCa        >  2:1570649/1‑150 (MQ=255)
                                                                                                                                        ttCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCAc       >  1:630427/1‑150 (MQ=255)
                                                                                                                                         tCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCAcc      >  2:309176/1‑150 (MQ=255)
                                                                                                                                            cgccgcATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCACCCAg   <  1:1200453/150‑1 (MQ=255)
                                                                                                                                             gccgcATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCACCCAGa  >  1:239241/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCACCCAGA  >  NC_000913/2945840‑2946130

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: