Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,190,398 T→G 51.5% M173R (ATG→AGG)  yqiH → putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,190,3980TG51.5% 24.9 / 269.7 186M173R (ATG→AGG) yqiHputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (48/42);  new base G (45/51);  total (93/93)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

Only 100 of 188 total aligned reads displayed.

GGCCAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAAT  >  NC_000913/3190299‑3190495
                                                                                                   |                                                                                                 
ggcCAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAt                                                                                                   <  1:1633777/100‑1 (MQ=255)
   cAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAATGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGc                                                                                                 >  2:2290091/1‑99 (MQ=255)
   cAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATg                                                                                                  <  1:1806033/98‑1 (MQ=255)
   cAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGcc                                                                                                <  1:2751998/100‑1 (MQ=255)
   cAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGcc                                                                                                <  1:2416651/100‑1 (MQ=255)
   cAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGc                                                                                                 >  2:1925957/1‑99 (MQ=255)
    aaaaGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCt                                                                                               >  2:1805148/1‑100 (MQ=255)
    aaaaGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCt                                                                                               <  1:3366729/100‑1 (MQ=255)
     aaaGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCtt                                                                                              <  2:9804/100‑1 (MQ=255)
      aaGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTa                                                                                             <  1:4131687/100‑1 (MQ=255)
      aaGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTa                                                                                             <  2:3470363/100‑1 (MQ=255)
       aGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTAt                                                                                            <  2:1967311/100‑1 (MQ=255)
         ccATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCa                                                                                          <  1:3811913/100‑1 (MQ=255)
         ccATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCa                                                                                          <  1:2445004/100‑1 (MQ=255)
          cATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCAc                                                                                         >  1:3078849/1‑100 (MQ=255)
          cATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCAGCTTTAACAATAAGAGGCCTTATAAc                                                                                         >  2:2009581/1‑100 (MQ=255)
           aTTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACg                                                                                        >  1:775715/1‑100 (MQ=255)
              gAATTAGAAGGGGATGGGGTAATGACTTACGAAATAGTTGATGTCATTCGTAGAAATGACGGTTCGATCTGCTTTAACAATAAGATGCCTTATCACGTTa                                                                                     <  2:2406770/100‑1 (MQ=255)
              gAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTa                                                                                     >  1:3402111/1‑100 (MQ=255)
                aTTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACg                                                                                   >  1:2076421/1‑100 (MQ=255)
                 ttAGACGACGATGGGGTAATGACTTACGAAAAAGTTGAACTCATTCGTAGAAATGATGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGc                                                                                  >  2:3638822/1‑100 (MQ=255)
                   agacgaCGATAGGGTAATGACTTACGAAAAATTTGAAGTCATTAGTAGAAATGACGCTTCGATCAGCTTTAACAATAAGATACCTTATCACGGTagacg                                                                                 >  2:3833182/1‑95 (MQ=255)
                   agaagaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTg                                                                                >  2:4025247/1‑100 (MQ=255)
                   agaagaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTg                                                                                <  1:1529924/100‑1 (MQ=255)
                      agaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGt                                                                             <  1:1087853/100‑1 (MQ=255)
                      agaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGt                                                                             <  1:892334/100‑1 (MQ=255)
                      agaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGt                                                                             <  1:3514880/100‑1 (MQ=255)
                        aCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTta                                                                           <  1:2604297/100‑1 (MQ=255)
                          gATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTtata                                                                         >  2:3148965/1‑100 (MQ=255)
                            tGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATc                                                                       <  1:642186/100‑1 (MQ=255)
                              gggTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCgg                                                                     >  2:4117112/1‑100 (MQ=255)
                               ggTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGa                                                                    >  1:383801/1‑100 (MQ=255)
                               ggTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGa                                                                    >  1:184313/1‑100 (MQ=255)
                               ggTAAGGACTTACGAAAAAGTTGAAATCATTCGGAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCgg                                                                     <  2:3079243/99‑1 (MQ=255)
                                gTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGaa                                                                   >  2:3886993/1‑100 (MQ=255)
                                gTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGaa                                                                   >  2:3496765/1‑100 (MQ=255)
                                  aaTGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACa                                                                 >  2:4113681/1‑100 (MQ=255)
                                  aaTGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACa                                                                 >  2:2574666/1‑100 (MQ=255)
                                  aaTGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACa                                                                 <  2:744740/100‑1 (MQ=255)
                                  aaTGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACa                                                                 <  2:847033/100‑1 (MQ=255)
                                     gACTTACGAAAAAGTTGAAATCATTCGTAGCAAGGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATACCGGAACAAAt                                                              >  2:2726015/1‑100 (MQ=255)
                                       cTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATgg                                                            <  1:2860230/100‑1 (MQ=255)
                                       cTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATgg                                                            >  1:977501/1‑100 (MQ=255)
                                       cTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATgg                                                            >  1:3249546/1‑100 (MQ=255)
                                         tACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATgggg                                                          <  1:721859/100‑1 (MQ=255)
                                         tACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATgggg                                                          <  2:1202077/100‑1 (MQ=255)
                                          aCGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGt                                                         >  2:1693346/1‑100 (MQ=255)
                                           cGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTa                                                        >  1:473590/1‑100 (MQ=255)
                                           cGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTa                                                        >  2:1339276/1‑100 (MQ=255)
                                              aaaaGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACg                                                     <  1:3778397/100‑1 (MQ=255)
                                               aaaGTTGAAATCATTCGTAGAAATGACGGTTCGATCAGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGa                                                    <  2:906013/100‑1 (MQ=255)
                                                aaGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGAt                                                   >  1:1214076/1‑100 (MQ=255)
                                                 aGTTGAATTCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATg                                                  <  1:567743/100‑1 (MQ=255)
                                                 aGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATg                                                  >  1:1141235/1‑100 (MQ=255)
                                                 aGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATg                                                  <  2:3402091/100‑1 (MQ=255)
                                                     gAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTAc                                              >  1:2188382/1‑100 (MQ=255)
                                                     gAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTAc                                              <  1:1468694/100‑1 (MQ=255)
                                                     gAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAAcatcgttac                                              >  1:3551493/1‑92 (MQ=255)
                                                       aaTCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACca                                            >  1:3029998/1‑100 (MQ=255)
                                                       aaTCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACca                                            >  1:2731303/1‑100 (MQ=255)
                                                            ttCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAAc                                       <  2:109696/100‑1 (MQ=255)
                                                            ttCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAAc                                       <  1:738320/100‑1 (MQ=255)
                                                            ttCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAAc                                       <  1:2184368/100‑1 (MQ=255)
                                                              cGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGc                                     >  1:3895805/1‑100 (MQ=255)
                                                              cGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGCCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGc                                     <  1:3089960/100‑1 (MQ=255)
                                                               gTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCa                                    >  2:1311426/1‑100 (MQ=255)
                                                               gTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCa                                    >  2:3732603/1‑100 (MQ=255)
                                                                  gAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAg                                 >  1:2060046/1‑100 (MQ=255)
                                                                  gAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAg                                 <  1:3481524/100‑1 (MQ=255)
                                                                     aTGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCtt                              <  1:1673148/100‑1 (MQ=255)
                                                                        aCGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAAt                           <  1:1391528/100‑1 (MQ=255)
                                                                         cGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACACAAAGCTTAATg                          <  1:2523071/100‑1 (MQ=255)
                                                                         cGGTTCGATCCGCTTTAACAATAAGATGCCGTATAACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATg                          <  1:1975556/100‑1 (MQ=255)
                                                                           gTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGt                        <  2:3754397/100‑1 (MQ=255)
                                                                           gTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGt                        <  1:1334493/100‑1 (MQ=255)
                                                                           gTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGt                        >  1:4041227/1‑100 (MQ=255)
                                                                             tCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTGATGGTaa                      <  2:3537121/100‑1 (MQ=255)
                                                                              cGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAAc                     >  2:53527/1‑100 (MQ=255)
                                                                              cGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAAc                     <  1:2322066/100‑1 (MQ=255)
                                                                               gATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAAcc                    <  2:478199/100‑1 (MQ=255)
                                                                                aTCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAAccc                   >  1:850565/1‑100 (MQ=255)
                                                                                aTCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAAccc                   >  2:1747113/1‑100 (MQ=255)
                                                                                aTCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAAcc                    <  2:2465690/99‑1 (MQ=255)
                                                                                aTCCGCTTTAACAATAAGAGGCCTTATAACGTTACGCTGGGGTATACCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAAcc                    <  2:401403/99‑1 (MQ=255)
                                                                                   cGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCtt                >  1:3647653/1‑100 (MQ=255)
                                                                                    gCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCttt               >  1:1997379/1‑100 (MQ=255)
                                                                                      tttAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTa             <  2:147538/100‑1 (MQ=255)
                                                                                      tttAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTa             <  2:507317/100‑1 (MQ=255)
                                                                                       ttAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAg            <  1:4007797/100‑1 (MQ=255)
                                                                                         aaCAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCt          >  2:1308976/1‑100 (MQ=255)
                                                                                          aCAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTa         >  2:1609647/1‑100 (MQ=255)
                                                                                           cAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTAt        >  1:2659375/1‑100 (MQ=255)
                                                                                           cAATAAGAGGCCTTATCCCGTTACTCTGGGTTATATCGTAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTATCCCCTTTTAGCt          <  2:3670314/98‑1 (MQ=255)
                                                                                            aaTAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATg       <  1:1288902/100‑1 (MQ=255)
                                                                                             aTAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGc      <  2:2122205/100‑1 (MQ=255)
                                                                                             aTAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGc      >  2:2314170/1‑100 (MQ=255)
                                                                                             aTAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGc      >  1:2705203/1‑100 (MQ=255)
                                                                                             aTAAGATGCCTGATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGc      <  1:4043966/100‑1 (MQ=255)
                                                                                              tAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCa     >  1:3392162/1‑100 (MQ=255)
                                                                                                 gATGCCTTATCACGTTACGCTGGGATATATTGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGATATGCa     >  2:2819916/1‑97 (MQ=255)
                                                                                                   |                                                                                                 
GGCCAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAAT  >  NC_000913/3190299‑3190495

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: