Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,106,518 (C)6→7 coding (709/1167 nt) wbbH ← putative O‑antigen polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,106,5121.C99.3% 1215.5 / ‑0.4 287coding (715/1167 nt)wbbHputative O‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (145/140);  total (146/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 329 total aligned reads displayed.

CAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACAC  >  NC_000913/2106412‑2106607
                                                                                                      |                                                                                               
cAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCAt                                                                                                  <  1:3066353/101‑1 (MQ=255)
 aaaaGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATc                                                                                                <  1:3066224/101‑1 (MQ=255)
 aaaaGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATc                                                                                                >  2:3066865/1‑101 (MQ=255)
   aaGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATccc                                                                                              <  1:3066520/101‑1 (MQ=255)
    aGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATcccc                                                                                             >  1:3066869/1‑101 (MQ=255)
    aGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATcccc                                                                                             <  2:3066466/101‑1 (MQ=255)
    aGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATcccc                                                                                             <  2:3066545/101‑1 (MQ=255)
     gACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATccccc                                                                                            >  2:3066871/1‑101 (MQ=255)
       cATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATccccccc                                                                                          <  2:3066194/101‑2 (MQ=255)
        aTGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATccccccca                                                                                         <  1:3066312/101‑3 (MQ=255)
        aTGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATccccccca                                                                                         >  1:3066877/1‑99 (MQ=255)
         tGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATcccccccag                                                                                        >  2:3066878/1‑98 (MQ=255)
           agaACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATcccccccaggc                                                                                      <  1:3066409/101‑6 (MQ=255)
           agaACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATcccccccaggc                                                                                      >  1:3066883/1‑96 (MQ=255)
             aaCTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAg                                                                                     >  1:3066886/1‑101 (MQ=255)
             aaCTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAg                                                                                     <  2:3066285/101‑1 (MQ=255)
             aaCTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAg                                                                                     <  1:3066365/101‑1 (MQ=255)
               cTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTc                                                                                   <  1:3066581/101‑1 (MQ=255)
                tGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTcc                                                                                  >  2:3066890/1‑101 (MQ=255)
                tGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTcc                                                                                  <  2:3066335/101‑1 (MQ=255)
                  gCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCAc                                                                                <  1:3066457/101‑1 (MQ=255)
                  gCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCAc                                                                                <  1:3066384/101‑1 (MQ=255)
                   cAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACg                                                                               <  2:3065939/101‑1 (MQ=255)
                    agagTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGt                                                                              >  2:3066893/1‑101 (MQ=255)
                    agagTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGt                                                                              >  1:3066894/1‑101 (MQ=255)
                      agTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTaa                                                                            >  2:3066900/1‑101 (MQ=255)
                      agTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTaa                                                                            <  1:3066571/101‑1 (MQ=255)
                      agTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTaa                                                                            <  1:3066494/101‑1 (MQ=255)
                       gTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTaaa                                                                           <  2:3066321/101‑1 (MQ=255)
                        ttAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGt                                                                                    >  1:3066902/1‑91 (MQ=255)
                           gATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAaga                                                                       <  1:3066401/101‑1 (MQ=255)
                           gATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAaga                                                                       >  1:3066906/1‑101 (MQ=255)
                            aTACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAagag                                                                      <  1:3066477/101‑1 (MQ=255)
                            aTACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATcc                                                                                               <  2:3066486/75‑1 (MQ=255)
                             tACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGc                                                                     >  1:3066909/1‑101 (MQ=255)
                               cTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCat                                                                   <  1:3066288/101‑1 (MQ=255)
                                ttGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCata                                                                  >  2:3066912/1‑101 (MQ=255)
                                 tGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCatat                                                                 >  1:3066915/1‑101 (MQ=255)
                                    tGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACa                                                              <  2:3066671/101‑1 (MQ=255)
                                    tGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACa                                                              >  1:3066920/1‑101 (MQ=255)
                                     gAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACaa                                                             >  2:3066923/1‑101 (MQ=255)
                                       aaaTAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAgg                                                           >  1:3066925/1‑101 (MQ=255)
                                        aaTAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGga                                                          <  1:3066532/101‑1 (MQ=255)
                                         aTAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGgag                                                         <  2:3066150/101‑1 (MQ=255)
                                         aTAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGgag                                                         <  2:3066326/101‑1 (MQ=255)
                                         aTAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGgag                                                         <  1:3066240/101‑1 (MQ=255)
                                         aTAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATccccc                                                                                            >  1:3066929/1‑65 (MQ=255)
                                          tAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGgaga                                                        >  2:3066931/1‑101 (MQ=255)
                                           aaaaCTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTaaa                                                                           >  1:3066934/1‑81 (MQ=255)
                                            aaaCTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGTTAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAAt                                                      <  1:3066179/101‑1 (MQ=255)
                                             aaCTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATa                                                     >  2:3066939/1‑101 (MQ=255)
                                             aaCTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATa                                                     <  2:3066370/101‑1 (MQ=255)
                                             aaCTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATa                                                     <  2:3066388/101‑1 (MQ=255)
                                              aCTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAg                                                    >  2:3066943/1‑101 (MQ=255)
                                               cTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGa                                                   <  2:3066434/101‑1 (MQ=255)
                                               cTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGa                                                   >  1:3066944/1‑101 (MQ=255)
                                                tCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGaa                                                  <  2:3066453/101‑1 (MQ=255)
                                                tCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGaa                                                  <  2:3066231/101‑1 (MQ=255)
                                                 ccTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAac                                                 >  2:3066947/1‑101 (MQ=255)
                                                     aaaCGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGCATAGAACACCt                                             <  1:3066397/101‑1 (MQ=255)
                                                     aaaCGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCt                                             >  2:3066951/1‑101 (MQ=255)
                                                      aaCGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTa                                            <  1:3066857/101‑1 (MQ=255)
                                                        cgcgATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACa                                          >  1:3066953/1‑101 (MQ=255)
                                                          cgATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGc                                        >  2:3066959/1‑101 (MQ=255)
                                                           gATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCt                                       >  1:3066961/1‑101 (MQ=255)
                                                            aTTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTg                                      >  2:3066965/1‑101 (MQ=255)
                                                              tataGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGta                                    <  1:3066651/101‑1 (MQ=255)
                                                              tataGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGta                                    >  1:3066967/1‑101 (MQ=255)
                                                               ataGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGtaa                                   >  2:3066970/1‑101 (MQ=255)
                                                                 aGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGtaata                                 >  1:3066976/1‑101 (MQ=255)
                                                                 aGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGtaata                                 <  1:3066429/101‑1 (MQ=255)
                                                                 aGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGtaata                                 <  2:3066414/101‑1 (MQ=255)
                                                                  gggCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGtaataa                                <  2:3066510/101‑1 (MQ=255)
                                                                   ggCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAg                               <  2:3066560/101‑1 (MQ=255)
                                                                   ggCTGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAg                               <  1:3066621/101‑1 (MQ=255)
                                                                      tGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATa                            >  1:3066982/1‑101 (MQ=255)
                                                                      tGACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATa                            >  2:3066981/1‑101 (MQ=255)
                                                                       gACCAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATaa                           <  1:3066700/101‑1 (MQ=255)
                                                                         ccAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAACAGCGCATATACAAGGAGAATAGAacac                                               >  1:3066988/1‑79 (MQ=255)
                                                                         ccAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAAc                         <  1:3066637/101‑1 (MQ=255)
                                                                         ccAAATACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAAc                         >  2:3066986/1‑101 (MQ=255)
                                                                           aaaTACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACAt                       >  2:3066991/1‑101 (MQ=255)
                                                                            aaTACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATa                      >  2:3066994/1‑101 (MQ=255)
                                                                              tACATTGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAAt                    >  1:3067001/1‑101 (MQ=255)
                                                                              tACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGga                                                          <  1:3066106/63‑1 (MQ=255)
                                                                              tACATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAAt                    >  2:3066999/1‑101 (MQ=255)
                                                                                cATGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGt                  >  2:3067004/1‑101 (MQ=255)
                                                                                 aTGGATAGATAATATGCCAACCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAac                                                 <  2:3066458/69‑1 (MQ=255)
                                                                                  tGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGttt                <  2:3066496/101‑1 (MQ=255)
                                                                                  tGGATAGATAATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGttt                >  2:3067008/1‑101 (MQ=255)
                                                                                     atagataATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTAc             <  2:3066629/101‑1 (MQ=255)
                                                                                     atagataATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTAc             >  1:3067012/1‑101 (MQ=255)
                                                                                      tagataATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTAct            <  2:3066831/101‑1 (MQ=255)
                                                                                       agataATATGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTActc           >  2:3067016/1‑101 (MQ=255)
                                                                                          taataTGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTAt        >  2:3067020/1‑101 (MQ=255)
                                                                                           aataTGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTAtt       >  2:3067025/1‑101 (MQ=255)
                                                                                           aataTGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTAtt       >  1:3067021/1‑101 (MQ=255)
                                                                                            atatGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTa      <  1:3066514/101‑1 (MQ=255)
                                                                                             tatGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTaa     <  1:3066544/101‑1 (MQ=255)
                                                                                              atGCCATCCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAac    >  1:3067027/1‑101 (MQ=255)
                                                                                                      |                                                                                               
CAAAAGACATGAGAACTGGCAGAGTTAGATACTTGCTGAAAATAAAACTCCTGAAACGCGATTATAGGGCTGACCAAATACATGGATAGATAATATGCCATCCCCCCAGGCAGTCCACGTAAAAAGAGCATATACAAGGAGAATAGAACACCTACAGCTGTAATAAGATAAACATAATGTTTTACTCTATTAACAC  >  NC_000913/2106412‑2106607

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