Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,191,432 C→T 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,4320CT98.7% 532.4 / ‑2.9 159G72S (GGC→AGC) yehCputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (79/78);  total (79/80)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 172 total aligned reads displayed.

TTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGT  >  NC_000913/2191336‑2191528
                                                                                                |                                                                                                
tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAGAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtat                                                                                                <  1:2009266/99‑1 (MQ=255)
tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               >  2:1588079/1‑100 (MQ=255)
tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               >  1:2976414/1‑100 (MQ=255)
tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               >  1:2606017/1‑100 (MQ=255)
tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               >  1:197915/1‑100 (MQ=255)
 ttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCTCTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatca                                                                                              <  1:82014/100‑1 (MQ=255)
 ttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatca                                                                                              >  1:3227591/1‑100 (MQ=255)
 ttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatca                                                                                              <  2:1997439/100‑1 (MQ=255)
 ttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatca                                                                                              <  1:2180153/100‑1 (MQ=255)
  tttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCTCTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatcat                                                                                             <  1:1171270/100‑1 (MQ=255)
   ttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTatcatc                                                                                            <  1:1022596/100‑1 (MQ=255)
   ttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTatcatc                                                                                            <  1:101850/100‑1 (MQ=255)
    tACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTAtcatca                                                                                           <  1:2005579/100‑1 (MQ=255)
     aCTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCaa                                                                                          <  1:2914476/100‑1 (MQ=255)
     aCTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCaa                                                                                          >  2:1265487/1‑100 (MQ=255)
       ttGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATc                                                                                        <  1:89685/100‑1 (MQ=255)
        tgctgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCAGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATcc                                                                                       >  1:2438889/1‑100 (MQ=255)
          ctgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCAc                                                                                     >  2:238908/1‑100 (MQ=255)
          ctgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCAc                                                                                     <  2:696746/100‑1 (MQ=255)
            gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                   >  1:1581044/1‑100 (MQ=255)
             cccGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGcc                                                                                  <  2:39844/100‑1 (MQ=255)
              ccGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCt                                                                                 >  2:1774241/1‑100 (MQ=255)
               cGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTg                                                                                <  1:658521/100‑1 (MQ=255)
                ggAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgc                                                                               <  2:2548337/100‑1 (MQ=255)
                 gAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgca                                                                              <  2:1299365/100‑1 (MQ=255)
                      tGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCagcagc                                                                         <  1:1216437/100‑1 (MQ=255)
                      tGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCagcagc                                                                         >  1:399372/1‑100 (MQ=255)
                       gCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCg                                                                        <  1:2660427/100‑1 (MQ=255)
                        ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCga                                                                       >  2:2767042/1‑100 (MQ=255)
                        ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCga                                                                       <  1:2414650/100‑1 (MQ=255)
                          ccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaan                                                                     <  1:3090212/100‑2 (MQ=255)
                             aTTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAAc                                                                  <  2:2319279/100‑1 (MQ=255)
                              tttttGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACg                                                                 <  1:760074/100‑1 (MQ=255)
                               ttttGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACgg                                                                <  2:2655591/100‑1 (MQ=255)
                                tttGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGt                                                               >  2:2036952/1‑100 (MQ=255)
                                    ccACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGt                                                               <  2:2685661/96‑1 (MQ=255)
                                    ccACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTcc                                                           <  2:1867069/100‑1 (MQ=255)
                                    ccACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTcc                                                           >  1:141456/1‑100 (MQ=255)
                                     cACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCGGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTccc                                                          <  1:1718853/100‑1 (MQ=255)
                                      ttccgatctgTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCt                                                         <  1:3106452/91‑1 (MQ=255)
                                       cTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTg                                                        <  1:3146527/100‑1 (MQ=255)
                                         ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                      <  2:123595/100‑1 (MQ=255)
                                          gtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTg                                                        <  2:3325121/97‑1 (MQ=255)
                                           tggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtc                                                    <  1:422623/100‑1 (MQ=255)
                                           tggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtc                                                    <  1:3098768/100‑1 (MQ=255)
                                           tggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtc                                                    >  2:58687/1‑100 (MQ=255)
                                              cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                                 >  2:920704/1‑100 (MQ=255)
                                              cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                                 <  1:611206/100‑1 (MQ=255)
                                              cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                                 <  1:2574766/100‑1 (MQ=255)
                                              cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGCATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                                 >  1:1096165/1‑100 (MQ=255)
                                               gTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCtagatcggaa                                                >  2:3106502/1‑91 (MQ=255)
                                               gTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatca                                                <  1:2247969/100‑1 (MQ=255)
                                               gTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatca                                                >  1:3039993/1‑100 (MQ=255)
                                                ttAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCaa                                               <  2:454058/100‑1 (MQ=255)
                                                 tAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAAc                                              <  2:1974647/100‑1 (MQ=255)
                                                  aaCATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACt                                             >  1:417813/1‑100 (MQ=255)
                                                  aaCATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACt                                             <  2:2203940/100‑1 (MQ=255)
                                                   aCATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTg                                            <  2:2751362/100‑1 (MQ=255)
                                                    cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatcagatcg                                           >  2:77408/1‑95 (MQ=255)
                                                    cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                           >  2:35792/1‑100 (MQ=255)
                                                       gAAAGGAACCTGAATTTTTTCCGGTGGTAATGACGTATAGCTATCATCACTCCACGCCTGCAGCAGCGCAGAACGGTTTCCCTCGTTGATCAACTGCAcc                                        >  1:1530128/1‑100 (MQ=255)
                                                         aaGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCAt                                      <  1:1699622/100‑1 (MQ=255)
                                                         aaGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCAt                                      >  1:2012835/1‑100 (MQ=255)
                                                              aCCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAATAACGGTTTCCCTGGTTCATCAACTGCACCATCActt                                 >  2:3310538/1‑100 (MQ=255)
                                                              aCCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActt                                 >  2:1890991/1‑100 (MQ=255)
                                                               ccTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtc                                >  2:3250869/1‑100 (MQ=255)
                                                               ccTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtc                                <  2:807262/100‑1 (MQ=255)
                                                                cTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtc                                >  1:534606/1‑99 (MQ=255)
                                                                cTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACACCTGAAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtct                               >  2:1260953/1‑100 (MQ=255)
                                                                 tGAATTTTTTCTGGTGTTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                              <  1:2334416/100‑1 (MQ=255)
                                                                 tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                              <  2:1511722/100‑1 (MQ=255)
                                                                 tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                              <  2:2933582/100‑1 (MQ=255)
                                                                 tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                              >  2:3403421/1‑100 (MQ=255)
                                                                 tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                              <  1:1806540/100‑1 (MQ=255)
                                                                   aaTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTa                            <  1:146414/100‑1 (MQ=255)
                                                                     ttttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAt                           <  1:220885/99‑1 (MQ=255)
                                                                      tttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAttt                         <  1:2939/100‑1 (MQ=255)
                                                                       ttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAtttt                        >  2:1580321/1‑100 (MQ=255)
                                                                           cTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGc                    >  2:3285549/1‑100 (MQ=255)
                                                                            tggtggTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGcc                   <  2:1673516/100‑1 (MQ=255)
                                                                            tggtggTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGcc                   <  2:2821149/100‑1 (MQ=255)
                                                                             ggtggtAATGACGTATCGCTATCATCAATCGACGCCTGCAGCAGCGAAGAATGGTTTCCCTGGTTCATCAACTGCACGATCACTTCTTGATTTTCTGCCg                  <  1:2025328/100‑1 (MQ=255)
                                                                             ggtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCg                  >  2:1436122/1‑100 (MQ=255)
                                                                             ggtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCg                  <  2:776079/100‑1 (MQ=255)
                                                                             ggtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCg                  >  1:3389830/1‑100 (MQ=255)
                                                                             ggtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCg                  >  2:1776527/1‑100 (MQ=255)
                                                                              gtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCgg                 >  1:1146585/1‑100 (MQ=255)
                                                                              gtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCgg                 >  2:2200675/1‑100 (MQ=255)
                                                                                 gtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTa              <  2:2392473/100‑1 (MQ=255)
                                                                                  tAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTaa             <  2:1274087/100‑1 (MQ=255)
                                                                                   aaTGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTaaa            >  1:468141/1‑100 (MQ=255)
                                                                                    aTGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaat           <  2:2438922/100‑1 (MQ=255)
                                                                                    aTGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaat           >  2:2979811/1‑100 (MQ=255)
                                                                                     tGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaata          >  2:1293768/1‑100 (MQ=255)
                                                                                     tGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaata          >  2:1148817/1‑100 (MQ=255)
                                                                                        cGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATg       >  2:2670565/1‑100 (MQ=255)
                                                                                         gTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATgc      <  2:3125765/100‑1 (MQ=255)
                                                                                           aTCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATgcgc    >  1:127141/1‑100 (MQ=255)
                                                                                            tCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGcgcg   <  2:3389882/100‑1 (MQ=255)
                                                                                             cGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGt  >  2:2286887/1‑100 (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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: