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

CTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGA  >  NC_000913/1004054‑1004247
                                                                                                     |                                                                                               
cTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGc‑c                                                                                                  >  1:1216962/1‑100 (MQ=255)
  ttCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑gg                                                                                              <  1:312063/100‑1 (MQ=255)
  ttCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑gg                                                                                              >  1:2134496/1‑100 (MQ=255)
        cagacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATGAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑gggggg                                                                                          <  1:784211/98‑1 (MQ=255)
        cagacagaGAATAGATATCTATTTTACTTATTAATACCAGTCAGGAAGCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑GGGGGGGa                                                                                        <  1:227177/100‑1 (MQ=255)
          gacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAACTCAACCGCCCAACGCGTTGGGC‑C‑‑GGGGGGGAAg                                                                                      >  1:413432/1‑100 (MQ=255)
            cagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑ggggggggaagc                                                                                    >  2:1190268/1‑95 (MQ=255)
              gagaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTACCAGAAAAATAAACCGCCCAACCCGTTTGGC‑C‑‑GGGGGGGAAGCGTc                                                                                  >  1:1252114/1‑100 (MQ=255)
               agaATCGATATCTATTTTACTGAATAATACCGGTCAGGAATCCTGTCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCg                                                                                    <  1:749033/98‑1 (MQ=255)
                gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCg                                                                                 <  2:878660/100‑1 (MQ=255)
                gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCg                                                                                 <  1:1866643/100‑1 (MQ=255)
                gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCg                                                                                 <  2:2339134/100‑1 (MQ=255)
                  aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGAc                                                                               <  2:2575067/100‑1 (MQ=255)
                    cGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑GGGGGGGAAGCGTCGACAGt                                                                            >  1:505060/1‑100 (MQ=255)
                        atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGcc                                                                         <  2:155533/100‑1 (MQ=255)
                          cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAg                                                                       >  2:2177237/1‑100 (MQ=255)
                          cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGACAg                                                                       >  1:564802/1‑100 (MQ=255)
                           tATTTTACTGACTAATACCAGTCAGGAAGCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑GGGGGGGAAGCGTCGACAGTGcc                                                                         <  1:1736002/96‑1 (MQ=255)
                              tttACTGACTAATACCTGTCAGGAATCCTGGCTTATTAACAGTAACATCAACCGCCAAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCAAGCAc                                                                    <  1:216745/99‑1 (MQ=255)
                              tttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCAcc                                                                   <  2:148110/100‑1 (MQ=255)
                               ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAAACAACCGCCCAACGCGTTGGGC‑C‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCAt                                                                 <  2:1809924/100‑1 (MQ=255)
                                tACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCAt                                                                 <  2:1705091/100‑1 (MQ=255)
                                tACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCAt                                                                 >  1:2260796/1‑100 (MQ=255)
                                  cTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTa                                                               >  1:881611/1‑100 (MQ=255)
                                   tGACTAATACCAGTCAGGACTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CGGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTa                                                               <  2:289759/100‑1 (MQ=255)
                                   tGACTAATACCAGTCAGGAATCCTGGCGTAGTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTAc                                                              <  1:591830/100‑1 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTAc                                                              >  1:1029954/1‑99 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACt                                                             >  2:1790378/1‑100 (MQ=255)
                                     aCTAATTCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                                                            <  2:1064999/100‑1 (MQ=255)
                                      cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg                                                           >  2:2020140/1‑100 (MQ=255)
                                      cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg                                                           >  1:751625/1‑100 (MQ=255)
                                       tAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑GGGGGCGAAGCGTCGACAGTGCCAGCACCATTACTCGc                                                          >  1:2448420/1‑99 (MQ=255)
                                        gatctcAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CGGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGc                                                          <  1:861591/95‑1 (MQ=255)
                                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                         <  2:1586082/100‑1 (MQ=255)
                                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                         >  2:1086412/1‑100 (MQ=255)
                                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                         <  2:909680/100‑1 (MQ=255)
                                            ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTc                                                     >  1:1240438/1‑100 (MQ=255)
                                             cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CGGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcagatc                                                     >  2:861594/1‑95 (MQ=255)
                                             cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                    <  2:2196792/100‑1 (MQ=255)
                                             cAGTCAGGAATCCTGGCTTATTAACAGTAAAATAAACAGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTAGACAGTGCCAGCACCATTACTGGCCGCTcc                                                    >  1:1096698/1‑100 (MQ=255)
                                              aGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CGGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                    >  2:969009/1‑100 (MQ=255)
                                              aGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg                                                   >  2:1479295/1‑100 (MQ=255)
                                                tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                                 <  2:465359/100‑1 (MQ=255)
                                                 cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa                                                >  2:1635242/1‑100 (MQ=255)
                                                   ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt                                             >  1:120221/1‑100 (MQ=255)
                                                       tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtc                                         <  2:1033110/100‑1 (MQ=255)
                                                         cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtct                                        >  2:439777/1‑100 (MQ=255)
                                                         cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtctc                                       >  2:2389517/1‑100 (MQ=255)
                                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTcc                                      <  2:752698/100‑1 (MQ=255)
                                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTcc                                      >  2:673101/1‑100 (MQ=255)
                                                            gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                     >  1:1981110/1‑100 (MQ=255)
                                                            gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                     <  1:2091756/100‑1 (MQ=255)
                                                            gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                    <  1:531030/100‑1 (MQ=255)
                                                             cTTATTAACAGTAAAATCAACCGCCCAACGGGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                    <  2:1094286/100‑1 (MQ=255)
                                                              ntattaACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGa                                   >  1:332254/2‑100 (MQ=255)
                                                               tattaACAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGaa                                  >  1:1889368/1‑100 (MQ=255)
                                                                attaACAGTAAAATAAACCGCCCCACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCTAGCACCATAACTTGCCGCTCCGCCACTTATTCTCTTGAAg                                 >  1:395726/1‑100 (MQ=255)
                                                                   aaCAGTAAAATCAACCGCCCAACGTGTTGGGT‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                               <  2:1256915/99‑1 (MQ=255)
                                                                    aCAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCg                             >  2:5748/1‑100 (MQ=255)
                                                                     cAGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCgg                            <  1:1431526/100‑1 (MQ=255)
                                                                      aGTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGt                           >  2:681980/1‑100 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                          >  2:589577/1‑100 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                          <  2:1971280/100‑1 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                          >  2:307936/1‑100 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTATCCTGAAGCCCGGTa                          >  2:2558524/1‑100 (MQ=255)
                                                                          aaaTCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                       >  1:314668/1‑100 (MQ=255)
                                                                           cccccccggCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGAc                      >  1:676740/9‑100 (MQ=255)
                                                                            aTCAACCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACa                     <  2:332256/100‑1 (MQ=255)
                                                                              ccccggCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                   >  1:211139/6‑100 (MQ=255)
                                                                              ccccggCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                   >  1:2443223/6‑100 (MQ=255)
                                                                               cccggCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                  <  2:676743/96‑1 (MQ=255)
                                                                               aaCCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACGATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                  >  1:1613725/1‑100 (MQ=255)
                                                                               aaCCGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                    <  1:542081/98‑1 (MQ=255)
                                                                                ccggCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa                 >  2:1160349/4‑100 (MQ=255)
                                                                                ccggCCCAACGCGTTGGGC‑C‑‑GGGGGGGAAGCGTCGACCGTGCCAGCACCAGTACGGGCCGCCCCGCCACTTATCCTCCTGAAGCCCGGTACGACAGGc                  >  1:2088966/4‑98 (MQ=255)
                                                                                ccggCCCAACGCGTTGGGC‑C‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCAc                >  1:1294855/4‑100 (MQ=255)
                                                                                 cggCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCAc                >  2:1108614/3‑100 (MQ=255)
                                                                                  ggCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACg               >  2:215780/2‑100 (MQ=255)
                                                                                  cGCCCAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACg               >  2:446064/1‑100 (MQ=255)
                                                                                     ccAACGCGTTGGGCGG‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg            <  2:504608/100‑1 (MQ=255)
                                                                                     ccAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg            >  2:67752/1‑100 (MQ=255)
                                                                                     ccAACGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg            >  2:734647/1‑100 (MQ=255)
                                                                                       aaCGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcg          >  2:1807206/1‑100 (MQ=255)
                                                                                       aaCGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcg          >  2:1578687/1‑100 (MQ=255)
                                                                                       aaCGCGTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgc           >  1:720489/1‑99 (MQ=255)
                                                                                          gcgTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCtt       <  2:211141/100‑1 (MQ=255)
                                                                                           cgTTGGGC‑T‑‑GGGGGGGAAGCGTCGTCCGTGCCAGCACCATTACGGGCCGCTACGCCACTTCTTGTCCGGAAGCCTGGTACGACAGGCACGTTGCGCt        <  1:2020343/97‑1 (MQ=255)
                                                                                            gTTGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc     <  2:1473201/100‑1 (MQ=255)
                                                                                             ttGGGC‑CG‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCt    <  2:915700/100‑1 (MQ=255)
                                                                                              tGGGC‑C‑‑GGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTga  <  2:539861/100‑1 (MQ=255)
                                                                                                     |                                                                                               
CTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGC‑C‑‑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: ‑