Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,945,985 C→T 5.3% 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.3% 236.2 / 2.0 75intergenic (‑143/+96)tcdA/mltAtRNA threonylcarbamoyladenosine dehydratase/membrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base C (32/39);  new base T (2/2);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

CAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCACCCAGACA  >  NC_000913/2945836‑2946132
                                                                                                                                                     |                                                                                                                                                   
caGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                                     <  2:1273383/150‑1 (MQ=255)
  gACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                                                   >  1:374071/1‑150 (MQ=255)
    cATAGCACCTCAGATTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                 <  1:1075532/150‑1 (MQ=37)
      tAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                               >  2:250906/1‑150 (MQ=255)
       aGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                              >  1:1344526/1‑150 (MQ=255)
       aGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                              <  1:48919/150‑1 (MQ=255)
        gCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                             <  1:464125/150‑1 (MQ=255)
         cACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                            <  2:360109/150‑1 (MQ=255)
         cACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                            <  1:822464/150‑1 (MQ=255)
             tCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctacaa                                                                                                                                        <  1:255213/150‑7 (MQ=25)
             tCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                                                                                                                                        <  1:704978/150‑1 (MQ=255)
                gAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                                                                                                                     >  1:954684/1‑150 (MQ=255)
                  aTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcgc                                                                                                                                   >  1:1388236/1‑150 (MQ=255)
                  aTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcgc                                                                                                                                   >  2:1407541/1‑150 (MQ=255)
                    tAAAAAGGGCGAGGTTATCCCTCACACTAAACGCATAAATCCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGACCGGATAAGGCGTTCACGCCGCATCAGGCATGAACAACgcgc                                                                                                                                 >  2:1132237/1‑150 (MQ=255)
                     aaaaaGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCa                                                                                                                                <  2:469154/150‑1 (MQ=255)
                     aaaaaGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCa                                                                                                                                <  2:1095660/150‑1 (MQ=255)
                      aaaaGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCAc                                                                                                                               >  2:451438/1‑150 (MQ=255)
                       aaaGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACt                                                                                                                              <  1:756361/150‑1 (MQ=255)
                        aaGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACtt                                                                                                                             >  1:309494/1‑150 (MQ=255)
                          gggCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAg                                                                                                                           >  2:459023/1‑150 (MQ=255)
                           ggCGAGGTTATCCCTCACACTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGt                                                                                                                          >  1:666439/1‑150 (MQ=255)
                             cGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCa                                                                                                                        <  2:253791/150‑1 (MQ=255)
                               aGGTTATCCCTCACCCTAAACGCATAAATGCCTGATTCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCaat                                                                                                                      <  1:95675/150‑1 (MQ=255)
                                  ttATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCaataat                                                                                                                   >  1:338210/1‑150 (MQ=255)
                                        cTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAt                                                                                                             >  1:1110135/1‑150 (MQ=255)
                                            cccTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGa                                                                                                         >  2:1269218/1‑150 (MQ=255)
                                             ccTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAg                                                                                                        <  2:82425/150‑1 (MQ=255)
                                              cTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAgg                                                                                                       >  2:722055/1‑150 (MQ=255)
                                              cTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAgg                                                                                                       <  2:666439/150‑1 (MQ=255)
                                              cTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAgg                                                                                                       <  1:250906/150‑1 (MQ=255)
                                                   cGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTTCATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAcc                                                                                                  >  1:1038003/1‑150 (MQ=255)
                                                   cGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAcc                                                                                                  <  2:1388236/150‑1 (MQ=255)
                                                         aaTGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAc                                                                                            >  2:953266/1‑150 (MQ=255)
                                                         aaTGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAc                                                                                            <  1:1132237/150‑1 (MQ=255)
                                                           tGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAccc                                                                                          >  2:1392418/1‑150 (MQ=255)
                                                           tGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAccc                                                                                          <  2:1308880/150‑1 (MQ=255)
                                                           tGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAccc                                                                                          <  2:1344526/150‑1 (MQ=255)
                                                              cTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAAGTGAGGTTACCCCTCACCCTGt                                                                                       >  2:1113081/1‑150 (MQ=255)
                                                                      cTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCg                                                                               <  2:240880/150‑1 (MQ=255)
                                                                      cTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCg                                                                               <  2:16500/150‑1 (MQ=255)
                                                                        aCGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTa                                                                             >  1:558348/1‑150 (MQ=255)
                                                                          gCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTaaa                                                                           >  2:542804/1‑150 (MQ=255)
                                                                             tATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAAc                                                                        <  1:953266/150‑1 (MQ=255)
                                                                             tATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAAc                                                                        >  1:1264592/1‑150 (MQ=255)
                                                                             tATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAAc                                                                        <  2:1251993/150‑1 (MQ=255)
                                                                                cAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGc                                                                     <  1:497270/150‑1 (MQ=255)
                                                                                     cTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATcc                                                                                <  1:169569/134‑1 (MQ=255)
                                                                                     cTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCATGTCATATcc                                                                                >  2:169569/1‑134 (MQ=255)
                                                                                            aCTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATc                                                         >  2:109733/1‑150 (MQ=255)
                                                                                                cTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACAt                                                     <  1:722055/150‑1 (MQ=255)
                                                                                                       atTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGc                                              >  1:314242/1‑150 (MQ=255)
                                                                                                         tttGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATTAGTTTTTTATGTTAACGGTGACCATGTATTATAAGTAAAATCGTCATACTTAATATCAAATAAGctgtta                                            <  2:380113/150‑6 (MQ=12)
                                                                                                         tttGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTa                                            <  1:1407541/150‑1 (MQ=255)
                                                                                                           tGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTaaa                                          >  1:127073/1‑150 (MQ=255)
                                                                                                           tGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTaaa                                          <  2:374071/150‑1 (MQ=255)
                                                                                                            ggATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAg                                         <  2:954684/150‑1 (MQ=255)
                                                                                                                ttGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGt                                     <  1:231960/150‑1 (MQ=255)
                                                                                                                 tGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGtt                                    <  2:558348/150‑1 (MQ=255)
                                                                                                                  gCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCTCTAAAGACGTTa                                   <  2:842263/150‑1 (MQ=255)
                                                                                                                  gCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTa                                   >  2:8416/1‑150 (MQ=255)
                                                                                                                   cATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTAc                                  >  1:897199/1‑150 (MQ=255)
                                                                                                                     tGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCt                                <  1:454095/150‑1 (MQ=255)
                                                                                                                        tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTgcg                             >  2:1368727/1‑150 (MQ=255)
                                                                                                                           gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGccc                          >  1:234790/1‑150 (MQ=255)
                                                                                                                            catgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCg                         >  2:553384/4‑150 (MQ=38)
                                                                                                                              ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCggg                       <  2:338210/150‑1 (MQ=255)
                                                                                                                                ccGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGCCGTTACCTGCGCCCGGGGc                     >  2:737452/1‑150 (MQ=255)
                                                                                                                                   gATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGt                  <  1:1166270/150‑1 (MQ=255)
                                                                                                                                       aGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAGCGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTc              <  1:451438/150‑1 (MQ=255)
                                                                                                                                        ggCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCa             >  1:205021/1‑150 (MQ=255)
                                                                                                                                        ggCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCa             <  2:259512/150‑1 (MQ=255)
                                                                                                                                         gCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAg            <  2:234790/150‑1 (MQ=255)
                                                                                                                                          cGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGc           >  2:275408/1‑150 (MQ=255)
                                                                                                                                            ttCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCAc         <  1:542804/150‑1 (MQ=255)
                                                                                                                                             tCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCAcc        <  1:1392418/150‑1 (MQ=255)
                                                                                                                                                  ccgcATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTg                               <  2:406513/122‑1 (MQ=255)
                                                                                                                                                  ccgcATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTg                               >  1:406513/1‑122 (MQ=255)
                                                                                                                                                   cgcATCCGGCATGAACAACGCGCACTTACTCAATAATCAGCATGCGAGGATACCCCTCACCCTGTCATATCCCTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCACCCAGaca  <  1:796478/150‑1 (MQ=255)
                                                                                                                                                   cgcATCCGGCATGAACAACGCCCACTTAGTCAATAATCAGGATGCGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCACCCAGaca  <  1:801430/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCACCCAGACA  >  NC_000913/2945836‑2946132

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: