Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,004,160:1 (G)7→8 100% coding (214/711 nt) ycbF → putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,004,1531.G81.2% 240.8 / 51.5 85coding (207/711 nt)ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base . (8/8);  new base G (40/29);  total (48/37)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGA  >  NC_000913/1004054‑1004247
                                                                                                     |                                                                                               
cTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGcc                                                                                                  >  1:1216962/1‑100 (MQ=255)
  ttCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑gg                                                                                              <  1:312063/100‑1 (MQ=255)
  ttCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑gg                                                                                              >  1:2134496/1‑100 (MQ=255)
        cagacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATGAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑gggggg                                                                                          <  1:784211/98‑1 (MQ=255)
        cagacagaGAATAGATATCTATTTTACTTATTAATACCAGTCAGGAAGCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGa                                                                                        <  1:227177/100‑1 (MQ=255)
          gacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAACTCAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGAAg                                                                                      >  1:413432/1‑100 (MQ=255)
            cagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑ggggggggaagc                                                                                    >  2:1190268/1‑95 (MQ=255)
              gagaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTACCAGAAAAATAAACCGCCCAACCCGTTTGGCC‑‑GGGGGGGAAGCGTc                                                                                  >  1:1252114/1‑100 (MQ=255)
               agaATCGATATCTATTTTACTGAATAATACCGGTCAGGAATCCTGTCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCg                                                                                    <  1:749033/98‑1 (MQ=255)
                gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCg                                                                                 <  2:878660/100‑1 (MQ=255)
                gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCg                                                                                 <  1:1866643/100‑1 (MQ=255)
                gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCg                                                                                 <  2:2339134/100‑1 (MQ=255)
                  aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGAc                                                                               <  2:2575067/100‑1 (MQ=255)
                    cGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACAGt                                                                            >  1:505060/1‑100 (MQ=255)
                        atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGcc                                                                         <  2:155533/100‑1 (MQ=255)
                          cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAg                                                                       >  2:2177237/1‑100 (MQ=255)
                          cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGACAg                                                                       >  1:564802/1‑100 (MQ=255)
                           tATTTTACTGACTAATACCAGTCAGGAAGCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACAGTGcc                                                                         <  1:1736002/96‑1 (MQ=255)
                              tttACTGACTAATACCTGTCAGGAATCCTGGCTTATTAACAGTAACATCAACCGCCAAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCAAGCAc                                                                    <  1:216745/99‑1 (MQ=255)
                              tttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCAcc                                                                   <  2:148110/100‑1 (MQ=255)
                               ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAAACAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCAt                                                                 <  2:1809924/100‑1 (MQ=255)
                                tACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                                                                 <  2:1705091/100‑1 (MQ=255)
                                tACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCAt                                                                 >  1:2260796/1‑100 (MQ=255)
                                  cTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTa                                                               >  1:881611/1‑100 (MQ=255)
                                   tGACTAATACCAGTCAGGACTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTa                                                               <  2:289759/100‑1 (MQ=255)
                                   tGACTAATACCAGTCAGGAATCCTGGCGTAGTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTAc                                                              <  1:591830/100‑1 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTAc                                                              >  1:1029954/1‑99 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACt                                                             >  2:1790378/1‑100 (MQ=255)
                                     aCTAATTCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                                                            <  2:1064999/100‑1 (MQ=255)
                                      cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg                                                           >  2:2020140/1‑100 (MQ=255)
                                      cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg                                                           >  1:751625/1‑100 (MQ=255)
                                       tAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑GGGGGCGAAGCGTCGACAGTGCCAGCACCATTACTCGc                                                          >  1:2448420/1‑99 (MQ=255)
                                        gatctcAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGc                                                          <  1:861591/95‑1 (MQ=255)
                                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                         <  2:1586082/100‑1 (MQ=255)
                                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                         >  2:1086412/1‑100 (MQ=255)
                                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                         <  2:909680/100‑1 (MQ=255)
                                            ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTc                                                     >  1:1240438/1‑100 (MQ=255)
                                             cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcagatc                                                     >  2:861594/1‑95 (MQ=255)
                                             cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                    <  2:2196792/100‑1 (MQ=255)
                                             cAGTCAGGAATCCTGGCTTATTAACAGTAAAATAAACAGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTAGACAGTGCCAGCACCATTACTGGCCGCTcc                                                    >  1:1096698/1‑100 (MQ=255)
                                              aGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                    >  2:969009/1‑100 (MQ=255)
                                              aGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg                                                   >  2:1479295/1‑100 (MQ=255)
                                                tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                                 <  2:465359/100‑1 (MQ=255)
                                                 cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa                                                >  2:1635242/1‑100 (MQ=255)
                                                   ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt                                             >  1:120221/1‑100 (MQ=255)
                                                       tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtc                                         <  2:1033110/100‑1 (MQ=255)
                                                         cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtct                                        >  2:439777/1‑100 (MQ=255)
                                                         cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtctc                                       >  2:2389517/1‑100 (MQ=255)
                                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTcc                                      <  2:752698/100‑1 (MQ=255)
                                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTcc                                      >  2:673101/1‑100 (MQ=255)
                                                            gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                     >  1:1981110/1‑100 (MQ=255)
                                                            gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                     <  1:2091756/100‑1 (MQ=255)
                                                            gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                    <  1:531030/100‑1 (MQ=255)
                                                             cTTATTAACAGTAAAATCAACCGCCCAACGGGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                    <  2:1094286/100‑1 (MQ=255)
                                                              ntattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGa                                   >  1:332254/2‑100 (MQ=255)
                                                               tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGaa                                  >  1:1889368/1‑100 (MQ=255)
                                                                attaACAGTAAAATAAACCGCCCCACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCTAGCACCATAACTTGCCGCTCCGCCACTTATTCTCTTGAAg                                 >  1:395726/1‑100 (MQ=255)
                                                                   aaCAGTAAAATCAACCGCCCAACGTGTTGGGTCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                               <  2:1256915/99‑1 (MQ=255)
                                                                    aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCg                             >  2:5748/1‑100 (MQ=255)
                                                                     cAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCgg                            <  1:1431526/100‑1 (MQ=255)
                                                                      aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGt                           >  2:681980/1‑100 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                          >  2:589577/1‑100 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                          <  2:1971280/100‑1 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                          >  2:307936/1‑100 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTATCCTGAAGCCCGGTa                          >  2:2558524/1‑100 (MQ=255)
                                                                          aaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                       >  1:314668/1‑100 (MQ=255)
                                                                           cccccccggCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGAc                      >  1:676740/9‑100 (MQ=255)
                                                                            aTCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACa                     <  2:332256/100‑1 (MQ=255)
                                                                              ccccggCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                   >  1:211139/6‑100 (MQ=255)
                                                                              ccccggCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                   >  1:2443223/6‑100 (MQ=255)
                                                                               cccggCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                  <  2:676743/96‑1 (MQ=255)
                                                                               aaCCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACGATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                  >  1:1613725/1‑100 (MQ=255)
                                                                               aaCCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                    <  1:542081/98‑1 (MQ=255)
                                                                                ccggCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa                 >  2:1160349/4‑100 (MQ=255)
                                                                                ccggCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACCGTGCCAGCACCAGTACGGGCCGCCCCGCCACTTATCCTCCTGAAGCCCGGTACGACAGGc                  >  1:2088966/4‑98 (MQ=255)
                                                                                ccggCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCAc                >  1:1294855/4‑100 (MQ=255)
                                                                                 cggCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCAc                >  2:1108614/3‑100 (MQ=255)
                                                                                  ggCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACg               >  2:215780/2‑100 (MQ=255)
                                                                                  cGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACg               >  2:446064/1‑100 (MQ=255)
                                                                                     ccAACGCGTTGGGCGG‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg            <  2:504608/100‑1 (MQ=255)
                                                                                     ccAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg            >  2:67752/1‑100 (MQ=255)
                                                                                     ccAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg            >  2:734647/1‑100 (MQ=255)
                                                                                       aaCGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcg          >  2:1807206/1‑100 (MQ=255)
                                                                                       aaCGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcg          >  2:1578687/1‑100 (MQ=255)
                                                                                       aaCGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgc           >  1:720489/1‑99 (MQ=255)
                                                                                          gcgTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCtt       <  2:211141/100‑1 (MQ=255)
                                                                                           cgTTGGGCT‑‑GGGGGGGAAGCGTCGTCCGTGCCAGCACCATTACGGGCCGCTACGCCACTTCTTGTCCGGAAGCCTGGTACGACAGGCACGTTGCGCt        <  1:2020343/97‑1 (MQ=255)
                                                                                            gTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc     <  2:1473201/100‑1 (MQ=255)
                                                                                             ttGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCt    <  2:915700/100‑1 (MQ=255)
                                                                                              tGGGCC‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTga  <  2:539861/100‑1 (MQ=255)
                                                                                                     |                                                                                               
CTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGA  >  NC_000913/1004054‑1004247

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