Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 934,749 A→G Q509R (CAA→CGA)  ftsK → cell division DNA translocase FtsK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913934,7490AG99.7% 1252.6 / ‑4.2 305Q509R (CAA→CGA) ftsKcell division DNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (166/138);  total (166/139)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

Only 100 of 316 total aligned reads displayed.

TGTGGAGCCTGAACCCGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGTACCTCCA  >  NC_000913/934651‑934846
                                                                                                  |                                                                                                 
tgtgGAGCCTGAACCCGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGaacgac                                                                                                 >  1:1735944/1‑101 (MQ=255)
tgtgGAGCCTGAACCCGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGaacgac                                                                                                 >  2:1735945/1‑101 (MQ=255)
     aGCCTGAACCCGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCg                                                                                            <  1:1735516/101‑1 (MQ=255)
      gCCTGAACCCGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCgg                                                                                           >  1:1735952/1‑101 (MQ=255)
        cTGAACCCGTTGTAGAAGAGACAAAACCCTCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGcc                                                                                         <  2:1735471/101‑1 (MQ=255)
        cTGAACCCGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGcc                                                                                         >  1:1735957/1‑101 (MQ=255)
        cTGAACCCGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGcc                                                                                         >  2:1735956/1‑101 (MQ=255)
         tGAACCCGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCt                                                                                        <  1:1735792/101‑1 (MQ=255)
         tGAACCCGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCt                                                                                        <  1:1735448/101‑1 (MQ=255)
           aaCCCGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTgg                                                                                      <  1:1735432/101‑1 (MQ=255)
           aaCCCGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTgg                                                                                      >  1:1735962/1‑101 (MQ=255)
           aaCCCGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTgg                                                                                      >  2:1735963/1‑101 (MQ=255)
               cGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTgg                                                                                      >  2:1735970/1‑97 (MQ=255)
                gTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCa                                                                                 >  1:1735971/1‑101 (MQ=255)
                gTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCa                                                                                 <  2:1735600/101‑1 (MQ=255)
                 ttGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCaa                                                                                >  1:1735976/1‑101 (MQ=255)
                     agaagaGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGa                                                                            >  2:1735978/1‑101 (MQ=255)
                       aagaGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGAtt                                                                          <  1:1735552/101‑1 (MQ=255)
                         gagaCAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTcc                                                                        <  2:1735597/101‑1 (MQ=255)
                         gagaCAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTcc                                                                        >  2:1735982/1‑101 (MQ=255)
                              aaaaCCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAAc                                                                   >  1:1735987/1‑101 (MQ=255)
                              aaaaCCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAAc                                                                   >  2:1735988/1‑101 (MQ=255)
                               aaccccGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATc                                                                                  <  2:1735589/82‑1 (MQ=255)
                               aaaCCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCaa                                                                                <  1:1735289/87‑1 (MQ=255)
                               aaaCCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAAcc                                                                  <  2:1735629/101‑1 (MQ=255)
                                aaCCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTcc                                                                        <  1:1735361/94‑1 (MQ=255)
                                 aCCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCgg                                                                >  2:1735993/1‑101 (MQ=255)
                                  cccGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAAcc                                                                  <  2:1735573/98‑1 (MQ=255)
                                   ccGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAAcc                                                                  >  1:1735997/1‑97 (MQ=255)
                                    cgcgcgTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTa                                                             <  2:1735638/101‑1 (MQ=255)
                                     gcgcgTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCg                                                                             >  2:1735999/1‑84 (MQ=255)
                                     gcgcgTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTcc                                                                        <  2:1735695/89‑1 (MQ=255)
                                      cgcgTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGc                                                                                          >  2:1736002/1‑70 (MQ=255)
                                      cgcgTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTaaa                                                           <  2:1735522/101‑1 (MQ=255)
                                      cgcgTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTaaa                                                           <  2:1735490/101‑1 (MQ=255)
                                       gcgTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAag                                                          >  1:1736004/1‑101 (MQ=255)
                                        cgTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCg                                                                 <  2:1735650/93‑1 (MQ=255)
                                         gTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAagaa                                                        >  1:1736010/1‑101 (MQ=255)
                                         gTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAagaa                                                        <  2:1735750/101‑1 (MQ=255)
                                         gTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGCGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAagaa                                                        >  2:1736008/1‑101 (MQ=255)
                                          tCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAagaac                                                       <  1:1735635/101‑1 (MQ=255)
                                            cgccgcTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAacaacca                                                     >  1:1736013/1‑101 (MQ=255)
                                             gccgcTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAagaaccag                                                    >  2:1736017/1‑101 (MQ=255)
                                              ccgcTTTACTACTTTGAAGAAGTTTAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAagaaccaga                                                   <  2:1735669/101‑1 (MQ=255)
                                              ccgcTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTAt                                                                                   <  1:1735441/69‑1 (MQ=255)
                                              ccgcTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAagaaccaga                                                   >  1:1736019/1‑101 (MQ=255)
                                               cgcTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAagaaccagaa                                                  <  1:1735462/101‑1 (MQ=255)
                                                 cTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAagaaccagaacc                                                >  1:1736025/1‑101 (MQ=255)
                                                 cTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAagaaccagaacc                                                >  2:1736024/1‑101 (MQ=255)
                                                  tttACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGAt                                                                           >  1:1736030/1‑73 (MQ=255)
                                                  tttACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCg                                                                 <  1:1735709/83‑1 (MQ=255)
                                                  tttACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCg                                               <  2:1735529/101‑1 (MQ=255)
                                                   ttACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGa                                              >  2:1736031/1‑101 (MQ=255)
                                                     actactTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGTCTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATc                                            <  2:1735612/101‑1 (MQ=255)
                                                     actactTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGtt                                                              >  2:1736037/1‑83 (MQ=255)
                                                       tactTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCCGGTATCAACCg                                                                             >  1:1736038/1‑66 (MQ=255)
                                                          tttGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATc                                       <  2:1735433/101‑1 (MQ=255)
                                                           ttGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCt                                      >  2:1736041/1‑101 (MQ=255)
                                                            tGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCgg                                                                                           >  2:1736047/1‑47 (MQ=255)
                                                            tGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCg                                                                 <  1:1735539/73‑1 (MQ=255)
                                                            tGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAAt                                        <  1:1735682/98‑1 (MQ=255)
                                                            tGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGGATCAACCGATTCCAGAACCGGTTAAaga                                                         >  1:1736046/1‑81 (MQ=255)
                                                             gaagaaGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATc                                            <  2:1735461/93‑1 (MQ=255)
                                                             gaagaaGTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTc                                    <  1:1735483/101‑1 (MQ=255)
                                                              aagaagTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCg                                   >  2:1736051/1‑101 (MQ=255)
                                                               agaagTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGc                                  >  1:1736054/1‑101 (MQ=255)
                                                                 aagTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTg                                >  2:1736058/1‑101 (MQ=255)
                                                                  agTTGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGa                               >  1:1736060/1‑101 (MQ=255)
                                                                    ttGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGaaa                             >  2:1736064/1‑101 (MQ=255)
                                                                     tGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAg                            <  2:1735568/101‑1 (MQ=255)
                                                                     tGAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAg                            >  1:1736067/1‑101 (MQ=255)
                                                                      gAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGc                           >  1:1736073/1‑101 (MQ=255)
                                                                      gAAGAGAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGc                           >  2:1736071/1‑101 (MQ=255)
                                                                        agagAAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCAc                         >  2:1736078/1‑101 (MQ=255)
                                                                         gataaGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCAcc                        >  1:1736079/4‑101 (MQ=255)
                                                                          agaAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCt                       <  1:1735602/101‑1 (MQ=255)
                                                                          agaAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGAATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTg                                >  1:1736083/1‑92 (MQ=255)
                                                                           gaAGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCtt                      <  1:1735645/101‑1 (MQ=255)
                                                                            aaGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTc                     >  2:1736085/1‑101 (MQ=255)
                                                                            aaGCGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTc                     <  1:1735627/101‑1 (MQ=255)
                                                                               cGAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGt                  <  1:1735494/101‑1 (MQ=255)
                                                                                gAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGtt                 >  1:1736091/1‑101 (MQ=255)
                                                                                gAGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGtt                 >  2:1736090/1‑101 (MQ=255)
                                                                                 aGCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTg                <  2:1735677/101‑1 (MQ=255)
                                                                                  gCCCGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTgc               >  1:1736097/1‑101 (MQ=255)
                                                                                   cccGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTgca              <  1:1735389/101‑1 (MQ=255)
                                                                                    ccGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCt                                 <  1:1735718/81‑1 (MQ=255)
                                                                                    ccGTGAACGTGAACGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTgcag             >  2:1736098/1‑101 (MQ=255)
                                                                                     cgtgaacgtgaaCGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTgcagc            >  1:1736101/1‑101 (MQ=255)
                                                                                       tgaacgtgaacgACTTGCGGCCTTGTATCAACCGATTCCAGaa                                                                    >  2:1736108/1‑43 (MQ=255)
                                                                                       tgaacgtgaacgACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGcagcag          <  1:1735817/101‑1 (MQ=255)
                                                                                       tgaacgtgaacgACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGcagcag          <  2:1735739/101‑1 (MQ=255)
                                                                                       tgaacgtgaacgACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGcagcag          >  1:1736106/1‑101 (MQ=255)
                                                                                       tgaacgtgaacgACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGcagcag          <  2:1735543/101‑1 (MQ=255)
                                                                                        gaacgtgaacGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGt         <  2:1735657/101‑1 (MQ=255)
                                                                                          acgtgaacGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGTAc       >  1:1736111/1‑101 (MQ=255)
                                                                                           cgtgaacGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGTAcc      >  1:1736114/1‑101 (MQ=255)
                                                                                            gtgaacGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGTACCt     >  1:1736117/1‑101 (MQ=255)
                                                                                             tgaacGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGTACCTc    <  2:1735795/101‑1 (MQ=255)
                                                                                             tgaacGACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGTACCTc    <  2:1735762/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
TGTGGAGCCTGAACCCGTTGTAGAAGAGACAAAACCCGCGCGTCCGCCGCTTTACTACTTTGAAGAAGTTGAAGAGAAGCGAGCCCGTGAACGTGAACAACTTGCGGCCTGGTATCAACCGATTCCAGAACCGGTTAAAGAACCAGAACCGATCAAATCTTCGCTGAAAGCACCTTCTGTTGCAGCAGTACCTCCA  >  NC_000913/934651‑934846

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