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,4320CT99.4% 549.4 / ‑4.7 163G72S (GGC→AGC) yehCputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (74/88);  total (75/88)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 166 total aligned reads displayed.

TTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGT  >  NC_000913/2191336‑2191528
                                                                                                |                                                                                                
tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               >  2:1779842/1‑100 (MQ=255)
tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               <  2:370499/100‑1 (MQ=255)
  tttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatcat                                                                                             >  2:3491339/1‑100 (MQ=255)
   ttACTTGCTGCCTGTATTTGGGCCCTATTTTTGCCGCTGGTGGCGTTAACATGAAAGGAACCTGATTTTTTTCTGGTGGTAATGACGTATCGCTatcatc                                                                                            <  1:2003828/100‑1 (MQ=255)
      cTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGTTTTTTATTACGTATCGCTATCATCAAt                                                                                         >  1:2659935/1‑100 (MQ=255)
      cTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAAt                                                                                         >  2:535415/1‑100 (MQ=255)
       ttGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATc                                                                                        <  2:2468672/100‑1 (MQ=255)
       ttGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATc                                                                                        >  2:2489090/1‑100 (MQ=255)
           tgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATc                                                                                        <  1:474918/96‑1 (MQ=255)
           tgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACg                                                                                    <  2:3301427/100‑1 (MQ=255)
            gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                   >  2:1169902/1‑100 (MQ=255)
            gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                   >  1:4032487/1‑100 (MQ=255)
            gcCCGGAATTTGCCCCTATTGTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                   >  2:4022708/1‑100 (MQ=255)
             cccGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGcc                                                                                  <  2:465654/100‑1 (MQ=255)
             cccGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGcc                                                                                  >  2:1886240/1‑100 (MQ=255)
               cGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGcc                                                                                  <  1:1548830/98‑1 (MQ=255)
                ggAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgc                                                                               <  1:294611/100‑1 (MQ=255)
                ggAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgc                                                                               <  2:2704719/100‑1 (MQ=255)
                   aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                            <  2:4252935/100‑1 (MQ=255)
                   aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                            >  1:1579281/1‑100 (MQ=255)
                   aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                            <  2:370741/100‑1 (MQ=255)
                   aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                            <  1:1879940/100‑1 (MQ=255)
                   aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                            >  2:2059090/1‑100 (MQ=255)
                   aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                            <  2:2307531/100‑1 (MQ=255)
                   aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgca                                                                              >  1:3031114/1‑98 (MQ=255)
                   aTTTGCCCCTATTTTTGCCACTGGTGGCGATAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                            <  2:3200010/100‑1 (MQ=255)
                      tGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCagcagc                                                                         >  2:1117664/1‑100 (MQ=255)
                        ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTTTTTTTTTTATGACGTATCGCTATCATCAATCCACGCCTGCAGCATCga                                                                       >  1:689858/1‑100 (MQ=255)
                        ccccTATTTATGCCACTGGTGGCGTTATCATGAAAGGCACTTGAATTTTTTCTTGTGATAATGACCTATCGCTATCTTCAATCCACGCCTGCAGCAACg                                                                        <  1:1676450/99‑1 (MQ=255)
                         cccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaa                                                                      >  2:2913002/1‑100 (MQ=255)
                          ccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaag                                                                     >  2:2413463/1‑100 (MQ=255)
                          ccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaag                                                                     >  2:3526110/1‑100 (MQ=255)
                          ccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaag                                                                     <  2:1626880/100‑1 (MQ=255)
                             aTTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAAc                                                                  <  2:2406991/100‑1 (MQ=255)
                             aTTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAAc                                                                  <  2:1291312/100‑1 (MQ=255)
                               ttttGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACgg                                                                >  2:2163617/1‑100 (MQ=255)
                               ttttGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACgg                                                                <  1:1716622/100‑1 (MQ=255)
                                    ccACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTcc                                                           >  1:1030005/1‑100 (MQ=255)
                                     cACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTccc                                                          <  2:2584002/100‑1 (MQ=255)
                                      aCTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTcc                                                           <  1:3835730/98‑1 (MQ=255)
                                        tggtggCGTTAACCTGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTgg                                                       <  2:2933132/100‑1 (MQ=255)
                                        tggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTgg                                                       <  2:3087921/100‑1 (MQ=255)
                                         ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                      <  2:934306/100‑1 (MQ=255)
                                         ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                      <  1:1995619/100‑1 (MQ=255)
                                          gtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGtt                                                     >  2:3451946/1‑100 (MQ=255)
                                          gtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGtt                                                     >  2:2256918/1‑100 (MQ=255)
                                           tggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtc                                                    <  1:3219255/100‑1 (MQ=255)
                                            ggCGTTAACATGAAAGGAACCTGAATTTTTTCTTGTAGTCATGACGTATCGCTATCATCCATCCACGCCTGCAGCGGCGAAGAACGGTTTCCCTGGTtca                                                   <  1:1055272/100‑1 (MQ=255)
                                            ggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtca                                                   <  2:2197648/100‑1 (MQ=255)
                                            ggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtca                                                   >  2:1954638/1‑100 (MQ=255)
                                              cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                                 >  2:2670886/1‑100 (MQ=255)
                                               gTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACTTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatca                                                <  2:1494550/100‑1 (MQ=255)
                                                  aaCATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACt                                             <  1:511767/100‑1 (MQ=255)
                                                  aaCATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACt                                             <  2:702377/100‑1 (MQ=255)
                                                  aaCATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCACCGAAGAACGGTTTCCCTGGTTCATCAACt                                             <  1:4059716/100‑1 (MQ=255)
                                                   aCATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTg                                            <  1:4114660/100‑1 (MQ=255)
                                                    cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                           >  2:1803934/1‑100 (MQ=255)
                                                    cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                           >  1:1184159/1‑100 (MQ=255)
                                                      tGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCAc                                         <  1:1089991/100‑1 (MQ=255)
                                                       gAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCAcc                                        >  1:4016867/1‑100 (MQ=255)
                                                        aaaGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCa                                       <  2:3610848/100‑1 (MQ=255)
                                                         aatgAACCTGAATTTTTTCTGGTGGTAATCACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCa                                       <  1:3457732/96‑1 (MQ=255)
                                                           ggAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCa                                    <  2:4032507/100‑1 (MQ=255)
                                                            gAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCAc                                   >  2:672444/1‑100 (MQ=255)
                                                            gAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCAc                                   <  1:2277887/100‑1 (MQ=255)
                                                            gAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCAc                                   >  2:2829850/1‑100 (MQ=255)
                                                             aaCCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCatt                                  >  2:795230/1‑98 (MQ=255)
                                                                 tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCGTCAACTGCACCATCActtctt                              >  2:3667529/1‑100 (MQ=255)
                                                                 tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                              >  1:3820800/1‑100 (MQ=255)
                                                                 tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                              >  1:262844/1‑100 (MQ=255)
                                                                 tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                              >  2:168128/1‑100 (MQ=255)
                                                                  gAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCttt                             >  2:785144/1‑100 (MQ=255)
                                                                  gAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCttt                             <  1:4259472/100‑1 (MQ=255)
                                                                    aTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAt                           >  2:277731/1‑100 (MQ=255)
                                                                    aTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTACAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAt                           >  2:146617/1‑100 (MQ=255)
                                                                        tttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTc                       <  2:1452898/100‑1 (MQ=255)
                                                                         ttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGATCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCt                      <  1:369400/100‑1 (MQ=255)
                                                                         ttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCt                      <  1:371696/100‑1 (MQ=255)
                                                                         ttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCt                      <  1:2680169/100‑1 (MQ=255)
                                                                          tCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTg                     <  1:1102077/100‑1 (MQ=255)
                                                                           cTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTgt                    >  1:3214398/1‑99 (MQ=255)
                                                                             ggtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCg                  <  1:2840998/100‑1 (MQ=255)
                                                                             ggtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCg                  >  2:4011323/1‑100 (MQ=255)
                                                                              gtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCgg                 >  2:3624194/1‑100 (MQ=255)
                                                                              gtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCgg                 <  1:3893538/100‑1 (MQ=255)
                                                                               tggtAATGACGTATCGCTATAATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTg                     >  1:4224254/1‑95 (MQ=255)
                                                                                ggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGt               <  2:633285/100‑1 (MQ=255)
                                                                                ggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGt               >  1:3793168/1‑100 (MQ=255)
                                                                                ggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGt               <  2:812489/100‑1 (MQ=255)
                                                                                ggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCGTTATTTTCTGCCgg                 >  1:2493128/1‑98 (MQ=255)
                                                                                  tAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTaa             >  1:3927241/1‑100 (MQ=255)
                                                                                     tGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaata          >  1:798893/1‑100 (MQ=255)
                                                                                     tGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaata          >  1:2602706/1‑100 (MQ=255)
                                                                                     tGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGGAaata          >  2:194749/1‑100 (MQ=255)
                                                                                      gACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaataa         >  2:2157701/1‑100 (MQ=255)
                                                                                       aCGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaataat        <  2:3166777/100‑1 (MQ=255)
                                                                                        cGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaataatn       <  2:2524566/100‑2 (MQ=255)
                                                                                         gTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATgc      <  1:2413457/100‑1 (MQ=255)
                                                                                          tATCGCTATCATCAATTCACGCCAGCAGCAGCTAATAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTGATTTTCTGCCGGGTAAATAATgc      <  1:1615857/99‑1 (MQ=255)
                                                                                             cGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGt  >  1:3227528/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: