Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,965,089 A→G 15.0% V368A (GTG→GCG) ‡ flhB ← flagellin export apparatus, substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,965,0890AG15.0% 304.6 / 40.6 132V368A (GTG→GCG) flhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (63/48);  new base G (13/7);  total (77/55)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

Only 100 of 134 total aligned reads displayed.

TGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTAC  >  NC_000913/1964992‑1965187
                                                                                                 |                                                                                                  
tGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGGTAATAAAATCCAGGGCTTCCGGCAcc                                                                                                  >  2:3381181/1‑100 (MQ=255)
 gCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCAcc                                                                                                  >  2:1244470/1‑99 (MQ=255)
 gCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCg                                                                                                 >  2:471551/1‑100 (MQ=255)
 gCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCg                                                                                                 >  1:2351805/1‑100 (MQ=255)
   cATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGa                                                                                               >  1:3961865/1‑100 (MQ=255)
      tgtgTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGa                                                                                            <  1:2910593/100‑1 (MQ=255)
       ntgtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATACAATCCAGGGCTTCCGGCACCGGAAGAt                                                                                           >  2:2988479/2‑100 (MQ=255)
        tgtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATg                                                                                          <  1:1104598/100‑1 (MQ=255)
        tgtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATg                                                                                          <  1:3384818/100‑1 (MQ=255)
        tgtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATg                                                                                          <  1:4135283/100‑1 (MQ=255)
        tgtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATg                                                                                          <  1:1343881/100‑1 (MQ=255)
         gtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGTGCCGGAAgccct                                                                                         <  1:3931930/100‑5 (MQ=255)
         gtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGTGCCGGAAgccct                                                                                         <  1:1142347/100‑5 (MQ=255)
         gtCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGTGCCGGAAgccct                                                                                         <  1:477261/100‑5 (MQ=255)
          tCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGTGCCGGAAgccctg                                                                                        <  1:3730666/100‑6 (MQ=255)
          tCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGTGCCGGAAgccctg                                                                                        <  2:2621180/100‑6 (MQ=255)
           cGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGt                                                                                       >  2:1104303/1‑100 (MQ=255)
             aTTTCAGATTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAg                                                                                     >  2:2166519/1‑100 (MQ=255)
              tttCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAgg                                                                                    <  1:3578566/100‑1 (MQ=255)
               ttCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGt                                                                                   <  2:1108609/100‑1 (MQ=255)
                 cAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTg                                                                                 <  1:536682/100‑1 (MQ=255)
                 cAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTg                                                                                 <  2:1324880/100‑1 (MQ=255)
                  aGGTTTGCGGGCAGGCGCAGCAGCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTg                                                                                 <  2:786681/99‑1 (MQ=255)
                    gTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTAc                                                                              >  1:858392/1‑100 (MQ=255)
                      tggCGTGCAGGGGCAGCATCGCGGCCAGAGTACTCATGGGTCGGTTTATCGTTAATAAGATCCAGGGCTGCCGGCACCGGAAGATGAGTAGGTTGTACAg                                                                            <  2:655515/98‑1 (MQ=255)
                         cGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGgc                                                                         >  1:3705992/1‑100 (MQ=255)
                         cGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATACAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGgc                                                                         >  2:2163352/1‑100 (MQ=255)
                            gCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCCCCGGTAGATGAGTAGATTGTACAGGGCGCt                                                                      >  2:954133/1‑100 (MQ=255)
                             cAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTg                                                                     <  1:2334764/100‑1 (MQ=255)
                              aGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGt                                                                    <  2:111622/100‑1 (MQ=255)
                               ggCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTc                                                                   >  1:2749681/1‑100 (MQ=255)
                               ggCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTc                                                                   <  1:471555/100‑1 (MQ=255)
                               ggCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTc                                                                   >  1:423604/1‑100 (MQ=255)
                                gcgcAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTcc                                                                  <  2:545536/100‑1 (MQ=255)
                                gcgcAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTcc                                                                  >  2:4026936/1‑100 (MQ=255)
                                  gcagcaTCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTccac                                                                >  1:3572510/1‑100 (MQ=255)
                                  gcagcaTCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTccac                                                                <  1:734603/100‑1 (MQ=255)
                                   cagcaTCGCGGCCAGATTACTCATGGGCCGGTTTCTCGTTATTACAATTCAGGGCTTCCGGCACCGGCCGATGGGTGGGTTGTTCTGGGCGCCGtgcctc                                                               <  2:2406537/100‑6 (MQ=255)
                                        tCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATCAAATCACGTGCTTCCCGCACCGGAAGATGAGTCGGTTGTACAGGGCGCTCTCCAccag                                                             <  2:2463830/97‑1 (MQ=255)
                                        tCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCAccagcca                                                          >  2:210791/1‑100 (MQ=255)
                                        tCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCAccagcc                                                           <  2:3939271/99‑1 (MQ=255)
                                           cgGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGAc                                                       >  1:3770737/1‑100 (MQ=255)
                                             gCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGc                                                     >  1:3533631/1‑100 (MQ=255)
                                             gCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGc                                                     >  2:553762/1‑100 (MQ=255)
                                               cAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCa                                                   <  1:769448/100‑1 (MQ=255)
                                                 gATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGc                                                 >  1:3445090/1‑100 (MQ=255)
                                                    tACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGtt                                              <  1:1065852/100‑1 (MQ=255)
                                                      cTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTc                                            <  2:2821207/100‑1 (MQ=255)
                                                      cTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTc                                            >  2:3257590/1‑100 (MQ=255)
                                                        cATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAg                                          <  1:3577731/100‑1 (MQ=255)
                                                           gggTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTg                                       <  1:946258/100‑1 (MQ=255)
                                                            ggTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGc                                      <  2:2321517/100‑1 (MQ=255)
                                                            ggTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGc                                      >  1:1510904/1‑100 (MQ=255)
                                                            ggTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGc                                      <  2:2558151/100‑1 (MQ=255)
                                                             gTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGcc                                     >  1:809244/1‑100 (MQ=255)
                                                             gTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGcc                                     <  2:2123725/100‑1 (MQ=255)
                                                                ggTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGa                                  >  2:3035181/1‑100 (MQ=255)
                                                                ggTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGa                                  >  2:1842033/1‑100 (MQ=255)
                                                                ggTTTCTCGTTAATAAAATCCAGGGCGTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCACCGTCTCAGTTGCCAGa                                  >  2:1917426/1‑100 (MQ=255)
                                                                 gTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGAc                                 <  2:347949/100‑1 (MQ=255)
                                                                    tctcGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCa                              >  2:3370126/1‑100 (MQ=255)
                                                                     ctcGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCa                              >  2:861967/1‑99 (MQ=255)
                                                                     ctcGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAg                             >  1:2855165/1‑100 (MQ=255)
                                                                      tcGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAgg                            >  1:2637241/1‑100 (MQ=255)
                                                                       cGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGc                           >  2:2276973/1‑100 (MQ=255)
                                                                         ttAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGcc                          <  2:2809268/99‑1 (MQ=255)
                                                                         ttAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCa                         >  1:2665016/1‑100 (MQ=255)
                                                                          taataaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCan                        <  1:2988497/100‑2 (MQ=255)
                                                                          taataaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAg                        <  1:2431857/100‑1 (MQ=255)
                                                                          taataaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAg                        >  1:4016334/1‑100 (MQ=255)
                                                                            ataaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCa                      >  1:3979572/1‑100 (MQ=255)
                                                                            ataaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCa                      >  1:1790121/1‑100 (MQ=255)
                                                                            ataaAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCa                      >  1:362176/1‑100 (MQ=255)
                                                                                  tactcatcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCg                >  1:521950/8‑100 (MQ=255)
                                                                                  tCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCg                <  1:472235/100‑1 (MQ=255)
                                                                                  tCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCg                <  2:2007589/100‑1 (MQ=255)
                                                                                  tCCAGGGCTTCCGGCACCAGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCg                <  2:3246813/100‑1 (MQ=255)
                                                                                   actcatcTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGc               >  2:4030897/7‑100 (MQ=255)
                                                                                    cAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGcc              >  1:959962/1‑100 (MQ=255)
                                                                                    cAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGcc              >  1:959952/1‑100 (MQ=255)
                                                                                    cAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGcc              >  1:1903294/1‑100 (MQ=255)
                                                                                     tcatcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCa             >  2:2081633/5‑100 (MQ=255)
                                                                                     tcatcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCa             >  2:81556/5‑100 (MQ=255)
                                                                                     aGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGACCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCa             >  1:822523/1‑100 (MQ=255)
                                                                                      gggCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCAc            <  1:152024/100‑1 (MQ=255)
                                                                                      gggCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCAc            >  1:1234546/1‑100 (MQ=255)
                                                                                      gggCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCAc            >  1:3639950/1‑100 (MQ=255)
                                                                                      catcTTCCGGTGCCGGAAGATGCGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCAc            >  2:3805878/4‑100 (MQ=255)
                                                                                      catcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGTCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCAc            >  2:2084008/4‑100 (MQ=255)
                                                                                      catcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCAc            >  1:1955180/4‑100 (MQ=255)
                                                                                      catcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCAc            >  1:640167/4‑100 (MQ=255)
                                                                                      catcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCAc            >  1:1205798/4‑100 (MQ=255)
                                                                                      catcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCAc            >  1:1743477/4‑100 (MQ=255)
                                                                                      catcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCAc            >  1:2680584/4‑100 (MQ=255)
                                                                                       atcTTCCGGTGCCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCAcc           >  1:3660409/3‑100 (MQ=255)
                                                                                        gCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACcg          >  2:1277102/1‑100 (MQ=255)
                                                                                              gggaCCGGAAGTTGAGTAGGGTGGACAGGGCGCTGTCCGCCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCg                <  2:3800895/85‑1 (MQ=255)
                                                                                              ggCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGt    >  1:3021581/1‑100 (MQ=255)
                                                                                              ggCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGt    >  2:243464/1‑100 (MQ=255)
                                                                                                cACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTAc  <  1:3081077/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
TGCCATTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTAC  >  NC_000913/1964992‑1965187

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