Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,134,118 A→C Q106P (CAG→CCG)  flgF → flagellar component of cell‑proximal portion of basal‑body rod

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,134,1180AC98.0% 697.9 / 0.6 198Q106P (CAG→CCG) flgFflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base C (84/110);  total (85/113)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

ACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGC  >  NC_000913/1134021‑1134214
                                                                                                 |                                                                                                
acCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCgg                                                                                                <  2:3058044/100‑1 (MQ=255)
   tCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTg                                                                                             <  1:3890145/100‑1 (MQ=255)
    cgcgCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGa                                                                                            >  1:5486867/1‑100 (MQ=255)
    cgcgCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGa                                                                                            >  1:284199/1‑100 (MQ=255)
     gcgcCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGAt                                                                                           <  1:5370774/100‑1 (MQ=255)
     gcgcCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGAt                                                                                           >  1:694166/1‑100 (MQ=255)
     gcgcCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGAt                                                                                           <  1:2324484/100‑1 (MQ=255)
       gcCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATcc                                                                                         <  1:373938/100‑1 (MQ=255)
          cGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCAc                                                                                      >  2:1703447/1‑100 (MQ=255)
           gCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCAcc                                                                                     <  2:4612660/100‑1 (MQ=255)
              ggACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCggg                                                                                  <  2:163825/100‑1 (MQ=255)
               gACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGc                                                                                 <  2:1530292/100‑1 (MQ=255)
               gACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGc                                                                                 <  1:2851587/100‑1 (MQ=255)
                aCGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCa                                                                                <  1:2120979/100‑1 (MQ=255)
                aCGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCa                                                                                <  1:4126241/100‑1 (MQ=255)
                aCGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCa                                                                                >  1:2615575/1‑100 (MQ=255)
                  gTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAAc                                                                              >  1:857386/1‑100 (MQ=255)
                   tCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACt                                                                             >  2:3530477/1‑100 (MQ=255)
                    cgcgTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTg                                                                            >  2:4334342/1‑100 (MQ=255)
                      cgTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGAc                                                                          <  2:126285/100‑1 (MQ=255)
                      cgTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGAc                                                                          <  2:4242617/100‑1 (MQ=255)
                      cgTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGAc                                                                          >  2:1144220/1‑100 (MQ=255)
                       gTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACa                                                                         >  2:5364076/1‑100 (MQ=255)
                        ttGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACaa                                                                        <  2:2740496/100‑1 (MQ=255)
                         tGCAGCAGGATGGCTGACTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAAt                                                                       <  1:4579760/100‑1 (MQ=255)
                          gcagcaGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAAtt                                                                      <  2:5393397/100‑1 (MQ=255)
                           cagcagGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACGATTc                                                                     <  1:3158278/100‑1 (MQ=255)
                           cagcagGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTc                                                                     <  1:2938988/100‑1 (MQ=255)
                           cagcagGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTc                                                                     <  1:574840/100‑1 (MQ=255)
                            agcagGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCa                                                                    <  2:2498971/100‑1 (MQ=255)
                               agGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGTAACTGACGATTCAggg                                                                 <  2:5322703/100‑1 (MQ=255)
                               agGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAggg                                                                 >  1:3894100/1‑100 (MQ=255)
                                 gATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGc                                                               <  2:2842148/100‑1 (MQ=255)
                                     gctggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATcc                                                           >  1:2350264/1‑100 (MQ=255)
                                       tggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCgg                                                         <  1:3359512/100‑1 (MQ=255)
                                        ggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGt                                                        >  2:740215/1‑100 (MQ=255)
                                        ggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGt                                                        <  1:5313029/100‑1 (MQ=255)
                                        ggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGt                                                        <  2:5286458/100‑1 (MQ=255)
                                          ctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGa                                                      <  1:5042876/100‑1 (MQ=255)
                                          ctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGa                                                      >  1:515362/1‑100 (MQ=255)
                                          ctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGa                                                      >  1:516513/1‑100 (MQ=255)
                                           tggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCCGTGAt                                                     <  2:4841841/100‑1 (MQ=255)
                                             gcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAg                                                   <  1:5095792/100‑1 (MQ=255)
                                              ncGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAgg                                                  >  1:4165820/2‑100 (MQ=255)
                                                 tGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGa                                               <  1:1484156/100‑1 (MQ=255)
                                                  gCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAg                                              <  1:2153900/100‑1 (MQ=255)
                                                   cAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAgg                                             >  2:4716093/1‑100 (MQ=255)
                                                    aGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGc                                            <  1:749776/100‑1 (MQ=255)
                                                     gACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCt                                           >  2:927429/1‑100 (MQ=255)
                                                     gACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCt                                           >  1:3640618/1‑100 (MQ=255)
                                                       ccGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTgg                                         >  2:5118863/1‑100 (MQ=255)
                                                       ccGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTgg                                         >  1:4711379/1‑100 (MQ=255)
                                                       ccGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTgg                                         >  2:1491967/1‑100 (MQ=255)
                                                       ccGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTgg                                         >  2:1285416/1‑100 (MQ=255)
                                                        cGCTGACGGCAGTGAAGTGTAAAGGCGTAATGGCAGCATTCAGGTTGATCCCACCGGTCAACTGACACGTCAGGGGCATCCGGTGATAGGCGAGGCTggc                                        <  2:834715/100‑2 (MQ=255)
                                                        cGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTggg                                        <  1:1180733/100‑1 (MQ=255)
                                                         gCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGc                                       <  2:4563599/100‑1 (MQ=255)
                                                          cTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGcc                                      <  2:1134685/100‑1 (MQ=255)
                                                            gACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCCTCCGGTGATAGGCGAGGCTGGGCCaa                                    >  1:4060322/1‑100 (MQ=255)
                                                            gACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCag                                    >  2:2618115/1‑99 (MQ=255)
                                                             aCGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAAt                                   >  1:4347577/1‑100 (MQ=255)
                                                             aCGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAAt                                   <  1:2643779/100‑1 (MQ=255)
                                                              cGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAAtt                                  >  2:615014/1‑100 (MQ=255)
                                                              cGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAAtt                                  >  1:178430/1‑100 (MQ=255)
                                                              cGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAAtt                                  >  2:3136662/1‑100 (MQ=255)
                                                                gCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGc                                >  1:2672192/1‑100 (MQ=255)
                                                                gCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGc                                <  1:3291723/100‑1 (MQ=255)
                                                                gCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGc                                <  2:730344/100‑1 (MQ=255)
                                                                 cAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCt                               >  2:3844433/1‑100 (MQ=255)
                                                                  aGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCtg                              >  1:474727/1‑100 (MQ=255)
                                                                  aGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCtg                              >  1:3162695/1‑100 (MQ=255)
                                                                    cGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCtgtn                            <  1:4423673/100‑2 (MQ=255)
                                                                    cGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCtgtg                            >  1:5291667/1‑100 (MQ=255)
                                                                    cGAAGGGTACACGCGTGATGGCAGCGTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGTTAGGCGAGGCTGGGCCAATTGCtgtg                            <  2:2042855/100‑1 (MQ=255)
                                                                     gAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGc                           >  2:3419644/1‑100 (MQ=255)
                                                                      aaGGGTATACGCGACATGGCAGCATCCCGGTTGATCCAACCGGGAAACTGACCATTCAGGGGCATCCGGTGATAGGCGAGGCTGCGCCAATTGCTGTGcc                          <  2:1363957/100‑1 (MQ=255)
                                                                       aGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCg                         <  2:1453558/100‑1 (MQ=255)
                                                                        gggTATACGCGTAATGGTAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCGATTGCTGTGCCgg                        >  1:145749/1‑100 (MQ=255)
                                                                         ggTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGa                       >  2:4073662/1‑100 (MQ=255)
                                                                         ggTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGa                       >  1:2053969/1‑100 (MQ=255)
                                                                           tataCGCGTAATGGCAGCATTCCGGTTGATCCCACGGGGCGACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCGCTTGCTGTGCCGGAAg                     >  2:1042487/1‑100 (MQ=255)
                                                                           tataCGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAg                     <  1:1317775/100‑1 (MQ=255)
                                                                           tataCGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGCGCCAATTGCTGTGCCGGAAg                     <  2:3720267/100‑1 (MQ=255)
                                                                            ataCGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCGACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAgg                    >  1:4863136/1‑100 (MQ=255)
                                                                            ataCGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAgg                    >  1:3106385/1‑100 (MQ=255)
                                                                            ataCGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAgg                    <  2:509579/100‑1 (MQ=255)
                                                                              aCGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAgggg                  <  1:1977062/100‑1 (MQ=255)
                                                                               cgcgTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGc                 <  1:605141/100‑1 (MQ=255)
                                                                               cgcgTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGACTGGGCCAATTGCTGTGCCGGAAGGGGc                 >  1:1334516/1‑100 (MQ=255)
                                                                                       ggCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGGAGGGGCGGAAATCa         >  2:2265520/1‑100 (MQ=255)
                                                                                       ggCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCa         >  2:1838813/1‑100 (MQ=255)
                                                                                       ggCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCa         >  1:1030153/1‑100 (MQ=255)
                                                                                        gcagcaTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCAc        <  1:1026197/100‑1 (MQ=255)
                                                                                           gcaTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTAt     >  2:1958597/1‑100 (MQ=255)
                                                                                            caTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTAtt    <  1:2198074/100‑1 (MQ=255)
                                                                                            caTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTAtt    <  2:2590035/100‑1 (MQ=255)
                                                                                             aTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATtg   >  1:23678/1‑100 (MQ=255)
                                                                                             aTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATtg   <  1:3558460/100‑1 (MQ=255)
                                                                                              ttCCGGTTGATCNCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATtgc  <  1:4929736/100‑1 (MQ=255)
                                                                                              ttCCGGTTGATCCCACCGGGCAACTGACAATTCAGAGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATtgc  <  2:4632705/100‑1 (MQ=255)
                                                                                                 |                                                                                                
ACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGC  >  NC_000913/1134021‑1134214

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: