Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,708,429 T→C 16.5% intergenic (‑724/+187) dppA ← / ← proK dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,708,4290TC16.5% 208.9 / 30.0 96intergenic (‑724/+187)dppA/proKdipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base T (37/43);  new base C (10/6);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTC  >  NC_000913/3708331‑3708528
                                                                                                  |                                                                                                   
gACGCTGACGCGTCTTATCAGGCCTACGTCGCGTCGAACCTCCCTGCCCCCGGAGGTTCTCCTCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTg                                                                                                    <  1:2559689/100‑1 (MQ=255)
gACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTg                                                                                                    <  1:462982/100‑1 (MQ=255)
gACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTg                                                                                                    >  1:2070174/1‑100 (MQ=255)
 aCGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAAAAACAAAAAAGCCTGCTCGTTGa                                                                                                   <  2:74487/100‑1 (MQ=255)
  cGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAg                                                                                                  >  1:81275/1‑100 (MQ=255)
  cGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAg                                                                                                  >  1:868465/1‑100 (MQ=255)
  cGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAg                                                                                                  <  1:1190089/100‑1 (MQ=255)
   gCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGc                                                                                                 <  1:1569377/100‑1 (MQ=255)
    cTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCa                                                                                                >  2:160523/1‑100 (MQ=255)
     gtcagcgTCGCATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCCTGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAg                                                                                               <  2:1364274/96‑1 (MQ=14)
        cgcgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCt                                                                                            <  2:996942/100‑1 (MQ=255)
        cgcgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCAACGAGCAGGCt                                                                                            <  1:1080719/100‑1 (MQ=17)
          cgccTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCttt                                                                                          <  2:1361779/97‑1 (MQ=255)
           gTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCtttt                                                                                         >  2:62015/1‑100 (MQ=255)
            tattATCAGGCCTACGTCGAGTCGAACTTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCAACGAGCAGGCt                                                                                            <  2:1056902/94‑1 (MQ=14)
            tCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTc                                                                                        <  1:121582/100‑1 (MQ=255)
             cTTATCAGGCCTACGTCGAGTTGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCg                                                                                       <  1:1934901/100‑1 (MQ=255)
              ttATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGa                                                                                      <  1:1681026/100‑1 (MQ=255)
               tATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGaa                                                                                     >  1:640373/1‑100 (MQ=255)
               tATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCAACGAGCAGGCttttttgc                                                                                     <  1:2341285/100‑5 (MQ=14)
               tATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCAACGAGCAGGCttttttg                                                                                      <  2:893524/99‑4 (MQ=14)
                aTCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAt                                                                                    >  2:2019917/1‑100 (MQ=255)
                aTCAGGCCTACGTCGACTCGAACCTCCCTGCCCCCAGAGGTTCTCATCCTCATGAGCGGCAGAAGCAAAAAAGCCTGCTCAACGAGCAGGCttttttgct                                                                                    <  2:1172537/100‑6 (MQ=14)
                 tCAGGCCTACGTCGATTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAtt                                                                                   <  2:1762932/100‑1 (MQ=255)
                   aGGCGTATGTCGAGTCGAACCTCCCTGCCCATGTAGGTTCTCATCCTCATGAGCTGTAGAAGAAAAATAGCATTCTCGTTGCGCAGGCTTTTCGAAtt                                                                                   <  2:1901493/98‑1 (MQ=11)
                   aGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTg                                                                                 >  1:1063932/1‑100 (MQ=255)
                   aGGCCTACGTCGAGTCGAACCTCCCTGCACCCGAAGTTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCAACGAGCAGGCttttttgcttct                                                                                 <  2:914396/100‑9 (MQ=14)
                    ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTgg                                                                                <  2:19444/100‑1 (MQ=255)
                     gCCTAGGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTgg                                                                                <  2:2304911/99‑1 (MQ=18)
                     ccttccGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCAACCTCATCAGCTGCATATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCCAATTTg                                                                                 <  2:511260/93‑1 (MQ=11)
                      ccTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTc                                                                              >  2:1950737/1‑100 (MQ=25)
                      ccTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTc                                                                              <  1:936358/100‑1 (MQ=25)
                       cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCg                                                                             <  1:2385339/100‑1 (MQ=25)
                       cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCg                                                                             >  1:1125235/1‑100 (MQ=25)
                          cGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTg                                                                          <  2:1383585/100‑1 (MQ=14)
                           gtcgagtcgaACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGa                                                                         <  2:2538713/100‑1 (MQ=11)
                           gtcgagtcgaACCTCACTCCCCCCGGAGGTTCTCATCCCCATGCGCTGCAGAAGCACAAACGTCTTCTCGTTGAGCACGCTTTTCGAATTTGGTctg                                                                            >  2:2423296/1‑95 (MQ=11)
                                        ncccTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTc                                                            >  1:2538571/2‑100 (MQ=18)
                                          cccgCCCCCGGAGGTTCTCATCCTCATGAGCTACAAATGCAAAAAAGCCTGCTCGCTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGc                                                          >  1:155171/4‑100 (MQ=11)
                                           cTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGcc                                                         <  1:1060463/100‑1 (MQ=255)
                                           cTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGcc                                                         <  1:1842037/100‑1 (MQ=255)
                                             gCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCAc                                                       <  1:1396575/100‑1 (MQ=255)
                                                 ccGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCg                                                   >  1:598482/1‑100 (MQ=255)
                                                  cGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGt                                                  <  2:307262/100‑1 (MQ=255)
                                                  cGGAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGt                                                  <  1:589404/100‑1 (MQ=255)
                                                    gAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTgg                                                <  1:1333546/100‑1 (MQ=255)
                                                    gAGGTTCTCATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTgg                                                <  1:2091413/100‑1 (MQ=255)
                                                         tctcATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGc                                           <  1:1839935/100‑1 (MQ=255)
                                                           tcATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTAGTCGCTACTTCGTGTCTCGc                                           >  2:89809/1‑98 (MQ=14)
                                                           tcATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGccc                                         <  2:1763265/100‑1 (MQ=255)
                                                           tcATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGccc                                         <  1:362392/100‑1 (MQ=255)
                                                           tcATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGcc                                          <  1:362923/99‑1 (MQ=255)
                                                            cATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCt                                        >  2:374930/1‑100 (MQ=255)
                                                            cATCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCt                                        >  2:248385/1‑100 (MQ=255)
                                                              tCCTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGc                                      >  2:2358467/1‑100 (MQ=255)
                                                                cTCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCgg                                    >  1:131136/1‑100 (MQ=255)
                                                                 tCATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCggg                                   >  1:1531560/1‑100 (MQ=255)
                                                                  cATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCggtc                                  <  1:2375652/100‑3 (MQ=255)
                                                                  cATGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGc                                  >  2:1910578/1‑100 (MQ=255)
                                                                   aTGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGcc                                 <  1:825820/100‑1 (MQ=255)
                                                                   aTGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGcc                                 <  1:1310525/100‑1 (MQ=255)
                                                                    tGAGCTGGAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCg                                <  2:1753346/100‑1 (MQ=255)
                                                                    tGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTg                                       >  2:282142/1‑93 (MQ=255)
                                                                    tGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCg                                <  1:1661348/100‑1 (MQ=255)
                                                                    tGAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCg                                <  1:1616381/100‑1 (MQ=255)
                                                                     gAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGt                               <  1:2358482/100‑1 (MQ=255)
                                                                     gAGCTGCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGt                               <  1:1142846/100‑1 (MQ=255)
                                                                          gCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTg                          >  1:101557/1‑100 (MQ=255)
                                                                          gCAGAAGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTg                          >  1:1568069/1‑100 (MQ=255)
                                                                             gaagCAAGAAACCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCGCTGCGGGCCGTTGCTg                          >  2:1340291/4‑97 (MQ=21)
                                                                             gaagCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGgcg                       <  2:1979025/97‑1 (MQ=255)
                                                                              aattcgaaaaGCCTGCTCAACGAGCAGGCTTGTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCTTTGCTGgc                        >  1:1271710/7‑98 (MQ=25)
                                                                              aagCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGAGAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGGGGGCGGTTGcggtcg                       >  2:1209198/3‑94 (MQ=21)
                                                                                ttcgaaaaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCaa                    >  1:342584/5‑100 (MQ=25)
                                                                                 tcgaaaaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGGGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAAc                   >  2:1707345/4‑100 (MQ=37)
                                                                                 cAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCTACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGgc                        >  2:1960491/1‑95 (MQ=255)
                                                                                 cAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCCAc                   >  1:2513643/1‑100 (MQ=255)
                                                                                  cgaaaaGCCTGCTCAACGAGCAGGCTTTTTGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGTCGCAACg                  >  2:41853/3‑100 (MQ=25)
                                                                                  cgaaaaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACg                  >  1:2544699/3‑100 (MQ=37)
                                                                                  cgaaaaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTAGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACg                  >  1:2087993/3‑100 (MQ=25)
                                                                                  cgaaaaGCCTGCTCAACGAGAAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACg                  >  2:2432097/3‑100 (MQ=37)
                                                                                  aaaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACg                  >  1:1531819/1‑100 (MQ=255)
                                                                                  aaaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTGGCTGGCGCAACg                  >  1:607300/1‑100 (MQ=255)
                                                                                   gaaaaGCCTGCTCAACGAGCAGGCTTTTCGACTTTGGTCGGTGATAGCGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGc                            >  2:1957546/2‑89 (MQ=25)
                                                                                   gaaaaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGt                 >  2:1992617/2‑100 (MQ=37)
                                                                                   aaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGt                 >  1:96953/1‑100 (MQ=255)
                                                                                     aaaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTtc               >  2:910456/1‑100 (MQ=37)
                                                                                        gCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGGGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTtctct            <  1:155234/100‑1 (MQ=255)
                                                                                        gCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTtctct            >  1:670230/1‑100 (MQ=255)
                                                                                         ccTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTAGTCGCTACTACGTGGGTGGCCCTGCGGGCGGTGGCTGGCGCACGGTtatctc           >  2:2467069/1‑100 (MQ=18)
                                                                                           tGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCtcgc         >  1:1910245/1‑100 (MQ=255)
                                                                                           tGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCtcgc         <  1:374934/100‑1 (MQ=255)
                                                                                           tGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCtcgc         >  2:858673/1‑100 (MQ=255)
                                                                                           tGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCtcgc         >  1:1957771/1‑100 (MQ=255)
                                                                                            gCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCtcgct        >  1:1435221/1‑100 (MQ=255)
                                                                                               cGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCtcgcttcg     <  2:963199/100‑1 (MQ=255)
                                                                                                 ttGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCtcgcttcgct   <  1:1218770/100‑1 (MQ=255)
                                                                                                  tGAGCAGGCGTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTc  <  1:884127/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTC  >  NC_000913/3708331‑3708528

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: