Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,452,339 A→G 38.9% intergenic (‑3/+17) yfcS ← / ← yfcU putative periplasmic pilin chaperone/pseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,452,3390AG38.9% 46.7 / 91.3 85intergenic (‑3/+17)yfcS/yfcUputative periplasmic pilin chaperone/pseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base G (16/17);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

TGGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATG  >  NC_000913/2452241‑2452434
                                                                                                  |                                                                                               
tGGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGa                                                                                                 <  2:2163506/99‑1 (MQ=255)
tGGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGAt                                                                                                >  1:2284176/1‑100 (MQ=255)
   cGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCAt                                                                                             >  1:1968676/1‑100 (MQ=255)
       tgtaaCGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCt                                                                                         <  2:1188027/96‑1 (MQ=255)
       aCAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCt                                                                                         <  2:925354/100‑1 (MQ=255)
        cAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCtt                                                                                        <  2:1938389/100‑1 (MQ=255)
        cAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCt                                                                                         >  2:1921811/1‑99 (MQ=255)
        cAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCtt                                                                                        <  2:1595822/100‑1 (MQ=255)
        cAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCtt                                                                                        <  2:1757223/100‑1 (MQ=255)
          gACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTaa                                                                                      <  1:312791/100‑1 (MQ=255)
           aCGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAAt                                                                                     <  2:1190104/100‑1 (MQ=255)
             gCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCa                                                                                   <  1:1306517/100‑1 (MQ=255)
                aGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACg                                                                                 >  1:541398/1‑99 (MQ=255)
                 gggCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGtt                                                                               >  1:319501/1‑100 (MQ=255)
                  ggCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTa                                                                              >  1:2189316/1‑100 (MQ=255)
                  ggCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTa                                                                              >  1:1040606/1‑100 (MQ=255)
                  ggCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTa                                                                              <  1:1854135/100‑1 (MQ=255)
                  ggCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTa                                                                              >  1:1090226/1‑100 (MQ=255)
                   gCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTAt                                                                             >  1:1363408/1‑100 (MQ=255)
                      gACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTgg                                                                          >  2:1216128/1‑100 (MQ=255)
                       aCAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGt                                                                         <  1:1245284/100‑1 (MQ=255)
                       aCAAACTGGTAGCACTAACCAGTAACGCCAGCGTCATAGCACCCAATTTAGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGt                                                                         >  2:1477951/1‑100 (MQ=255)
                        cAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGtt                                                                        >  2:1415594/1‑100 (MQ=255)
                         aaaCTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTa                                                                       <  2:171232/100‑1 (MQ=255)
                         aaaCTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAAGCAGACATGGGTCATTCCTTAATCACGTTATTGGTTa                                                                       >  1:1682631/1‑100 (MQ=255)
                           aCTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGc                                                                     >  1:2185691/1‑100 (MQ=255)
                            cTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCt                                                                    >  1:1398298/1‑100 (MQ=255)
                            cTGGTAGCGCTAAGCAGTAACGCAAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGc                                                                     <  2:511793/99‑1 (MQ=255)
                              ggTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGCTGTTGGTTAGCTga                                                                  >  1:1534430/1‑100 (MQ=255)
                               gTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTgaa                                                                 <  1:1921833/100‑1 (MQ=255)
                                tAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTgaaa                                                                >  1:983851/1‑100 (MQ=255)
                                 aGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTgaaac                                                               >  1:1490547/1‑100 (MQ=255)
                                    gcTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTgaaacaga                                                            >  1:829454/1‑100 (MQ=255)
                                      tAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCATTAATCACGTTATTGGTTAGATgaaacagaaa                                                          >  2:1487209/1‑100 (MQ=255)
                                       aaGCAGTAACGCCACCGTCATAGCGCCCATTTTGGAAGAGCAAAGTAAATCTGACATGGATCATTCCTTAATCTCGTTATTGGTTAGCTgaaacagaaac                                                         <  2:404575/100‑1 (MQ=255)
                                            gTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAgggg                                                    <  1:770866/100‑1 (MQ=255)
                                            gTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAgggg                                                    >  2:2120657/1‑100 (MQ=255)
                                             tAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAggggg                                                   <  1:1690911/100‑1 (MQ=255)
                                               acaccAGCGGCATAGAACCAAATTTGGCAGAGAAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGAt                                                 <  1:2150664/97‑1 (MQ=255)
                                               aCGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGAt                                                 >  2:1956015/1‑100 (MQ=255)
                                                cGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATc                                                <  2:311871/100‑1 (MQ=255)
                                                cGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATc                                                >  1:1248250/1‑100 (MQ=255)
                                                 gCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCg                                               >  1:1098790/1‑100 (MQ=255)
                                                  ccAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCgg                                              >  1:151205/1‑100 (MQ=255)
                                                   cAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggt                                             <  2:2114117/100‑1 (MQ=255)
                                                    aGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtg                                            >  2:1295564/1‑100 (MQ=255)
                                                      cGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtggg                                          <  1:326878/100‑1 (MQ=255)
                                                      cGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtggg                                          <  1:2193033/100‑1 (MQ=255)
                                                      cGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtggn                                          <  2:217183/100‑2 (MQ=255)
                                                      cGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtggg                                          <  2:2181669/100‑1 (MQ=255)
                                                        tCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtgggtg                                        <  1:566977/100‑1 (MQ=255)
                                                           tgggagaaaaTGTGGCATAGCAAATGAAATCGGAAGTGGATCAGTCCTGAATCCCGTTAGTGGTTGGCTGAAACAGAAACAGGGGTAGCggtgggtg                                        <  2:2000493/89‑1 (MQ=255)
                                                           tAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCa                                     >  1:2315044/1‑100 (MQ=255)
                                                            aGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCAc                                    <  1:1839793/100‑1 (MQ=255)
                                                                cccAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGt                                >  2:550942/1‑100 (MQ=255)
                                                                  cAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTg                              <  2:1090216/100‑1 (MQ=255)
                                                                      tttgCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTc                          <  2:2032893/97‑1 (MQ=255)
                                                                      ttGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTc                          <  1:2012828/100‑1 (MQ=255)
                                                                      ttGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTc                          <  2:1261422/100‑1 (MQ=255)
                                                                      ttGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTc                          <  2:1414575/100‑1 (MQ=255)
                                                                       tGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCt                         <  2:1038896/100‑1 (MQ=255)
                                                                        ggCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTc                          <  1:550945/98‑1 (MQ=255)
                                                                        ggCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTg                        >  1:263689/1‑100 (MQ=255)
                                                                        ggCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTTACCCGTTTCGTCTg                        >  1:813818/1‑100 (MQ=255)
                                                                        ggCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTg                        >  2:1863359/1‑100 (MQ=255)
                                                                         gCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGc                       <  1:2404199/100‑1 (MQ=255)
                                                                         gCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGc                       >  2:1646223/1‑100 (MQ=255)
                                                                         gCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGCTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTAACCTGTTGCGTCTGc                       >  2:527245/1‑100 (MQ=255)
                                                                           agagCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTc                     >  2:2096194/1‑100 (MQ=255)
                                                                            gagCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTcc                    <  2:1264567/100‑1 (MQ=255)
                                                                                aaaGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGaa                >  1:2081675/1‑100 (MQ=255)
                                                                                aaaGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGaa                >  2:2098071/1‑100 (MQ=255)
                                                                                aaaGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGaa                <  1:924460/100‑1 (MQ=255)
                                                                                  aGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATg              >  2:1167126/1‑100 (MQ=255)
                                                                                  aGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATg              <  2:2178595/100‑1 (MQ=255)
                                                                                  aGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATg              <  2:1248237/100‑1 (MQ=255)
                                                                                   gTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGa             <  2:1696137/100‑1 (MQ=255)
                                                                                    tAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGAc            <  1:1440391/100‑1 (MQ=255)
                                                                                     aaaTCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACa           <  2:2081653/100‑1 (MQ=255)
                                                                                      aaTCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAg          >  1:1245937/1‑100 (MQ=255)
                                                                                           gACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCttt     >  2:546132/1‑100 (MQ=255)
                                                                                            aCATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTa    <  2:2227517/100‑1 (MQ=255)
                                                                                             cATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTAt   >  2:1540170/1‑100 (MQ=255)
                                                                                             cATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTAt   <  2:358888/100‑1 (MQ=255)
                                                                                              aTGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATg  >  1:1044672/1‑100 (MQ=255)
                                                                                              aTGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTAt   <  2:587777/99‑1 (MQ=255)
                                                                                                  |                                                                                               
TGGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATG  >  NC_000913/2452241‑2452434

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: