Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,191,432 C→T 100% G72S (GGC→AGC)  yehC ← putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,191,4320CT96.1% 483.1 / 2.9 156G72S (GGC→AGC) yehCputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base C (4/3);  new base T (72/77);  total (76/80)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

TTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGT  >  NC_000913/2191336‑2191528
                                                                                                |                                                                                                
tttttACTTGTTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               <  1:3830946/100‑1 (MQ=255)
tttttACTTGCTGCCCGTAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               <  2:934185/100‑1 (MQ=255)
tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               >  1:1030328/1‑100 (MQ=255)
tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               <  2:3504653/100‑1 (MQ=255)
 ttttACTTTCTTCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCcatc                                                                                               <  1:2051840/99‑1 (MQ=255)
 ttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatca                                                                                              <  2:2468094/100‑1 (MQ=255)
 ttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatca                                                                                              >  2:914578/1‑100 (MQ=255)
 ttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatca                                                                                              <  1:955335/100‑1 (MQ=255)
   ttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTatcatc                                                                                            <  2:3038969/100‑1 (MQ=255)
     aCTTGCTTCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCaa                                                                                          <  2:639865/100‑1 (MQ=255)
     aCTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCaa                                                                                          <  1:3890935/100‑1 (MQ=255)
        tgctgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATcc                                                                                       <  1:2961120/100‑1 (MQ=255)
           tgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACg                                                                                    <  1:1752663/100‑1 (MQ=255)
           tgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACg                                                                                    <  1:3563235/100‑1 (MQ=255)
            gcCCGGAATTTGCCCCTATTTTTTTTTTTTTTTTCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                   >  2:1553955/1‑100 (MQ=255)
            gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                   >  2:2705990/1‑100 (MQ=255)
            gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                   >  2:681759/1‑100 (MQ=255)
            gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                   >  1:2403168/1‑100 (MQ=255)
             cccGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGcc                                                                                  >  1:1089492/1‑100 (MQ=255)
              ccGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCt                                                                                 >  1:2207653/1‑100 (MQ=255)
                   aTTTGCCGTTATTTTTGATACTGGTAGCGTTAACATGAAAGGACCGTGAATTTTTTCTGGTGGAAATGACGTATCGCTATCATCAATCCACGcc                                                                                  <  2:30692/94‑1 (MQ=255)
                   aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                            >  2:117000/1‑100 (MQ=255)
                   aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                            >  1:3294729/1‑100 (MQ=255)
                      tGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCagcagc                                                                         <  1:2299694/100‑1 (MQ=255)
                        ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCga                                                                       >  2:2181996/1‑100 (MQ=255)
                        ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCga                                                                       >  1:1462171/1‑100 (MQ=255)
                        ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCga                                                                       >  1:1459377/1‑100 (MQ=255)
                        ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCg                                                                        >  2:1762115/1‑99 (MQ=255)
                         cccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaa                                                                      <  2:3358671/100‑1 (MQ=255)
                          ccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCCATCCACGCCTGCAGCAGCgaag                                                                     >  2:3278886/1‑100 (MQ=255)
                              tttttGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACg                                                                 <  2:2105655/100‑1 (MQ=255)
                               ttttGCCACTGGTGGTGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACgg                                                                <  1:4111623/100‑1 (MQ=255)
                                tttgttacgggttgCGTTAAGTTGAAAGGAACTTGAATTTTTTCTGGGGGTAATGACGTTTTGCTATCTTCAATCCTCGCCTGCAGCAGCGAAGAACgg                                                                <  2:1218145/90‑1 (MQ=255)
                                  tGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGttt                                                             >  1:2666194/1‑100 (MQ=255)
                                  tGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGttt                                                             >  1:549332/1‑100 (MQ=255)
                                  tGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGttt                                                             >  2:2118407/1‑100 (MQ=255)
                                   gCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTc                                                            >  1:3505501/1‑100 (MQ=255)
                                     cACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTcc                                                           >  2:1887231/1‑99 (MQ=255)
                                     cACTGGTGGCATTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTccc                                                          >  2:1728281/1‑100 (MQ=255)
                                      aCTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGTTATCATCAATCCACGCCTGCAGCAGCGACGAACGGTTTCCCt                                                         <  1:2823067/100‑1 (MQ=255)
                                        tggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTgg                                                       >  2:2211/1‑100 (MQ=255)
                                         ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                      >  2:1820468/1‑100 (MQ=255)
                                         ggtggCGTTAACATGAAAGGAACCTGAATTTTTTATGGTGGTAATGACGTATAGCTACCATCAATCAACGCCTGCAGCAGAGAAGACCGGTTTCCCTGGt                                                      >  2:2480942/1‑100 (MQ=255)
                                          gtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGtt                                                     <  1:3909066/100‑1 (MQ=255)
                                           tggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtc                                                    >  1:4169244/1‑100 (MQ=255)
                                           tggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtc                                                    <  1:3431845/100‑1 (MQ=255)
                                            ggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtca                                                   >  2:3349455/1‑100 (MQ=255)
                                             gCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcat                                                  >  1:876313/1‑100 (MQ=255)
                                              cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                                 <  1:2132715/100‑1 (MQ=255)
                                               gTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatca                                                >  2:2849400/1‑100 (MQ=255)
                                               gTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatca                                                >  1:1297955/1‑100 (MQ=255)
                                               gTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatca                                                <  2:2404421/100‑1 (MQ=255)
                                               gTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatca                                                >  1:2705407/1‑100 (MQ=255)
                                                 tAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGTAGCAGCGAAGAACGGTTTCCCTGGTTCATCAAc                                              >  2:553221/1‑100 (MQ=255)
                                                  aaCATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACt                                             <  1:2659881/100‑1 (MQ=255)
                                                   aCATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTg                                            <  1:747539/100‑1 (MQ=255)
                                                    cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                           >  1:783009/1‑100 (MQ=255)
                                                    cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                           >  1:3123514/1‑100 (MQ=255)
                                                    cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                           >  2:420657/1‑100 (MQ=255)
                                                    cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                           >  2:3887560/1‑100 (MQ=255)
                                                     aTGAAAGGAACCTGACTTTTTTCTGGTGGTAATTACGTATCGCTATCATCAATCCACTCCTGCAGCAGGGAAGAACGGTTTCCCTGGTTCATCAACTGc                                           >  2:2919567/1‑99 (MQ=255)
                                                     aTGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCa                                          <  2:1171441/100‑1 (MQ=255)
                                                     aTGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCa                                          <  1:2159291/100‑1 (MQ=255)
                                                     aTGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCa                                          <  1:1630930/100‑1 (MQ=255)
                                                     aTGAAAGGAACCTGAATGTTTTCTGGTGCTAATGACGTATAGCTATCATCAATCCACGCGTGCAGGAGCGAAGAACGGTTTCCCTGGTTCCTCAACTGCa                                          >  2:938781/1‑100 (MQ=255)
                                                      ttaaaGGAACCTGAATTTTTTCTGGTGGTAATGCCGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                           <  1:2959015/96‑1 (MQ=255)
                                                       gAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCAcc                                        <  1:1127641/100‑1 (MQ=255)
                                                       gAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCAcc                                        >  2:44339/1‑100 (MQ=255)
                                                       gAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCAcc                                        <  1:3040483/100‑1 (MQ=255)
                                                         aaGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCAt                                      >  1:3195703/1‑100 (MQ=255)
                                                          aGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCAt                                      >  2:3062520/1‑99 (MQ=255)
                                                           ggAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCa                                    <  1:2233665/100‑1 (MQ=255)
                                                            gAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCAc                                   >  1:2678865/1‑100 (MQ=255)
                                                            gAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCAGCAACTGCACCATCAc                                   >  1:4113831/1‑100 (MQ=255)
                                                             aaCCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCCTCAATCCACGCCTGCCGCAGCGAAGAACGGTTTCCCTGGTTCATCACCTGCACCacccc                                   <  2:3323464/99‑5 (MQ=255)
                                                             aaCCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCAct                                  <  1:3462444/100‑1 (MQ=255)
                                                             aaCCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCAct                                  >  2:807854/1‑100 (MQ=255)
                                                               ccTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtc                                >  1:229621/1‑100 (MQ=255)
                                                               ccTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtc                                <  1:3049647/100‑1 (MQ=255)
                                                                cTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtct                               >  1:365791/1‑100 (MQ=255)
                                                                 tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGGAGGGAAGAACGGTTTCCATGGGTCATCAACTGCACCATCActtctt                              >  2:4046336/1‑100 (MQ=255)
                                                                  gAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCttt                             <  2:2226551/100‑1 (MQ=255)
                                                                    aTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTCTAt                           <  1:3062538/100‑1 (MQ=255)
                                                                     ttttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAtt                          >  2:258801/1‑100 (MQ=255)
                                                                     ttttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAtt                          >  1:2753596/1‑100 (MQ=255)
                                                                      tttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAttt                         <  2:2695896/100‑1 (MQ=255)
                                                                        tttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTc                       <  2:2787856/100‑1 (MQ=255)
                                                                         ttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCt                      <  2:2951313/100‑1 (MQ=255)
                                                                           cTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGTACCATCACTTCTTTATTTTCTGc                    <  1:2847947/100‑1 (MQ=255)
                                                                           cTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTg                     >  2:197667/1‑99 (MQ=255)
                                                                            tggtggTAATTACGTATCGCTATCATCAATCCACGTCTGCAGCAGCCAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGcc                   <  2:1374831/100‑1 (MQ=255)
                                                                              gtggtAACGACGTATCGCTATCCTCAATCCACGCCTGCAGCAGCTAAGACCGGTTTACCTTCTTCATCTACTGCACAATCACTTCTTTATTTTCtcac                   >  2:3946930/1‑95 (MQ=255)
                                                                                ggtAATGCCTTATGGCTATCATCAATCGACGCCTGCAGCAGCGAAGAAGGGTTTCTCTGATTCATCAACAGGACTAGCACTTCTTTATTTTCt                      >  2:3841039/1‑93 (MQ=255)
                                                                                  tAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTaa             <  1:2706006/100‑1 (MQ=255)
                                                                                    aTGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaat           >  1:557068/1‑100 (MQ=255)
                                                                                        cGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATg       >  2:3917006/1‑100 (MQ=255)
                                                                                         gTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATgn      <  2:3784470/100‑2 (MQ=255)
                                                                                          tATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATgcg     >  2:1429471/1‑100 (MQ=255)
                                                                                            tCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGcgcg   <  1:3955392/100‑1 (MQ=255)
                                                                                            tCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGcgcg   <  1:1073961/100‑1 (MQ=255)
                                                                                                |                                                                                                
TTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGT  >  NC_000913/2191336‑2191528

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: