Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,798,235 A→T 15.5% intergenic (+28/+4) waaL → / ← waaU O‑antigen ligase/lipopolysaccharide core biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,798,2350AT15.5% 391.5 / 59.0 154intergenic (+28/+4)waaL/waaUO‑antigen ligase/lipopolysaccharide core biosynthesis
Reads supporting (aligned to +/- strand):  ref base A (72/58);  new base T (14/10);  total (86/68)
Fisher's exact test for biased strand distribution p-value = 8.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

Only 100 of 163 total aligned reads displayed.

ATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTA  >  NC_000913/3798136‑3798333
                                                                                                   |                                                                                                  
aTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaa                                                                                                    >  1:3260126/1‑100 (MQ=255)
aTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaa                                                                                                    <  1:1030604/100‑1 (MQ=255)
 ttCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                   >  1:3615980/1‑100 (MQ=255)
   ccAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAATAATAAAATTAATTAGGAATAAAGAAATTTAAGAAATGAGTTaaaa                                                                                                  >  2:3086176/1‑99 (MQ=255)
   ccAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAc                                                                                                 <  1:2315008/100‑1 (MQ=255)
   ccAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAc                                                                                                 >  2:1691882/1‑100 (MQ=255)
   ccAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAc                                                                                                 >  1:534674/1‑100 (MQ=255)
   ccAATTATCATTATATCCGCTATAGGCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAc                                                                                                 <  2:770167/100‑1 (MQ=255)
    cAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACt                                                                                                >  2:1317907/1‑100 (MQ=255)
     aaTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAATAAGTGAGTTAAAACTc                                                                                               >  2:844757/1‑100 (MQ=255)
     aaTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTc                                                                                               <  2:534672/100‑1 (MQ=255)
      aTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                              <  1:444216/100‑1 (MQ=255)
          tCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActtc                                                                                          <  1:1112762/100‑1 (MQ=255)
           cATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActtct                                                                                         <  1:1132260/100‑1 (MQ=255)
           cATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtct                                                                                         <  1:1822789/100‑1 (MQ=255)
            aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtctt                                                                                        >  2:1103374/1‑100 (MQ=255)
            aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtctt                                                                                        >  2:2278438/1‑100 (MQ=255)
            aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtctt                                                                                        >  1:2292204/1‑100 (MQ=255)
             ttATATCCGCTATAGTCCTCTTACTAGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTa                                                                                       <  2:2897219/100‑1 (MQ=39)
             ttATATCCACTATAGTCCACTTACTCATCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTCAGTTAAAACTCACTTCTTa                                                                                       <  2:923259/100‑1 (MQ=255)
               atatCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaa                                                                                     <  1:1758794/100‑3 (MQ=255)
                  tCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCta                                                                                   <  2:689085/99‑1 (MQ=255)
                   cctatataGTCCTCTTACTCATCATTAATATTCGTAACGATCCGATTAATTAAGAATAAACGAGTTTAAGAAGTGAGTTAAATCTCACTTCTTATc                                                                                     <  2:1617058/92‑1 (MQ=255)
                    ctcTATAGTCCTCTTACTCGTCATTAAGAATCGTAACAATACAATTAATTAAGAATAAACATGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATAc                                                                                <  2:3477513/98‑1 (MQ=255)
                    cGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATAc                                                                                <  1:3558338/100‑1 (MQ=255)
                    cGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATAc                                                                                >  1:1248636/1‑100 (MQ=255)
                       tataGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAAc                                                                             <  1:2042057/100‑1 (MQ=255)
                       tataGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAATTATTATTTATACaa                                                                              >  2:2105755/1‑99 (MQ=255)
                            tactctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAt                                                                        <  1:1879496/98‑1 (MQ=255)
                               tctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAtctc                                                                     >  2:3240558/1‑100 (MQ=255)
                                ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATctct                                                                    >  1:965607/1‑100 (MQ=255)
                                  tACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                  <  2:726263/100‑1 (MQ=255)
                                    ctctTCATTAATAATCGTAACAATACAATTAATTAATAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTAc                                                                <  2:2888797/100‑1 (MQ=255)
                                      cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTg                                                              <  2:1840135/100‑1 (MQ=255)
                                       gTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGa                                                             >  2:3424248/1‑100 (MQ=255)
                                        tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGaa                                                            <  1:3392312/100‑1 (MQ=255)
                                        tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGaa                                                            <  1:3023224/100‑1 (MQ=255)
                                         cATTAATAATCGTAACAATACAATTAATTAAGAGTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt                                                           >  2:2179099/1‑100 (MQ=255)
                                         cATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGa                                                             >  2:40405/1‑98 (MQ=255)
                                         cATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt                                                           >  1:236266/1‑100 (MQ=255)
                                                 atCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaaa                                                   <  1:1571288/100‑1 (MQ=255)
                                                 atCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaaa                                                   <  2:3315381/100‑1 (MQ=255)
                                                  tCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAg                                                  <  2:3253090/100‑1 (MQ=255)
                                                  tCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAg                                                  <  1:167603/100‑1 (MQ=255)
                                                  tCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAg                                                  <  1:2417790/100‑1 (MQ=255)
                                                   cGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGACCTGAATTAAAACTCAAGTATTAGCTATACAAATTAATCTCTTTATTGCAtaaatagaa                                                   >  2:3572117/1‑90 (MQ=255)
                                                   cGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGt                                                 <  1:427222/100‑1 (MQ=255)
                                                    gTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAATCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGt                                                 >  1:154201/1‑99 (MQ=255)
                                                     tAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGttt                                               >  1:3375362/1‑100 (MQ=255)
                                                       acaatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCa                                             <  1:2179113/100‑1 (MQ=255)
                                                       acaatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCa                                             <  1:3424273/100‑1 (MQ=255)
                                                       acaatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCa                                             <  2:2573403/100‑1 (MQ=255)
                                                       acaatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTATTATCTATACAACTTAATTTCTTTAATAAATTAGTTAAAGtttaa                                             >  2:3509564/1‑98 (MQ=255)
                                                        caatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAg                                            <  2:1537369/100‑1 (MQ=255)
                                                         aatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGt                                           <  2:1464147/100‑1 (MQ=255)
                                                         aatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGt                                           >  1:913952/1‑100 (MQ=255)
                                                             caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                       >  1:2363408/1‑100 (MQ=255)
                                                             caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                       >  2:3258121/1‑100 (MQ=255)
                                                             caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                       <  2:1428527/100‑1 (MQ=255)
                                                             caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                       <  1:2191895/100‑1 (MQ=255)
                                                                ttaattaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGata                                    <  1:3258356/100‑1 (MQ=255)
                                                                ttaattaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGata                                    <  2:43647/100‑1 (MQ=255)
                                                                 taattaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGatat                                   >  2:1036552/1‑100 (MQ=255)
                                                                   attaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCt                                 >  1:160471/1‑100 (MQ=255)
                                                                     taaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCttt                               >  1:156961/1‑100 (MQ=255)
                                                                      aaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCtttt                              >  1:2244192/1‑100 (MQ=255)
                                                                        gAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTAc                            >  1:3095589/1‑100 (MQ=255)
                                                                          aTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAGCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCg                          <  1:2793662/100‑1 (MQ=255)
                                                                            aaaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGta                        >  1:814738/1‑100 (MQ=255)
                                                                            aaaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGta                        <  1:3316897/100‑1 (MQ=255)
                                                                             tgtatagatAAGAAGTGAGTTTTAACTCACTTCTTATATATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtat                       >  2:1941028/9‑100 (MQ=37)
                                                                             aaCAAGTTTAAGAAGTGAGTTAAAACTAACTTCTTATCTATACAACTTAATCTCTTTACAGAATTAGTTAAAGTATAAGTATCGATATCTTTTa                             >  2:3122313/1‑94 (MQ=255)
                                                                              aCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTGATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtatn                      <  1:1342579/100‑2 (MQ=255)
                                                                               tatagatAAGAAGTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTatat                     >  1:1068972/7‑100 (MQ=38)
                                                                                atagatAAGAAGTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATg                    >  1:1928914/6‑100 (MQ=38)
                                                                                  gTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCTTATATGtt                  >  2:2430065/1‑100 (MQ=255)
                                                                                  gTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtt                  >  2:13176/1‑100 (MQ=255)
                                                                                   tttAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGCATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                 >  1:2716852/1‑100 (MQ=255)
                                                                                   gatAAGAAGTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                 >  2:140677/3‑100 (MQ=255)
                                                                                    ttAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgg                <  2:1248628/100‑1 (MQ=255)
                                                                                     tAAGAAGTGAGTTAAAACTCACTTCTTATTTATACAACTTAATCTCTTTACTGTATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg               <  2:226136/100‑1 (MQ=255)
                                                                                     tAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg               >  1:3121367/1‑100 (MQ=255)
                                                                                      aagaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg              <  1:1354000/100‑1 (MQ=255)
                                                                                      aagaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg              >  1:2003954/1‑100 (MQ=255)
                                                                                       agacgTGAGTTAAAAATCACTTCTTATCTATAAAAATTAATCTCTTTACTGAATTAGTTAAAGTTACAGTATCGATATTTTTTCTCATATATGtt                  >  2:1188023/1‑95 (MQ=255)
                                                                                       agaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa             >  2:2785992/1‑100 (MQ=255)
                                                                                       agaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa             >  2:3064882/1‑100 (MQ=255)
                                                                                       agaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa             >  1:746964/1‑100 (MQ=255)
                                                                                       agaagTGAGGTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTAATGAATTAGTTAAAGTTTCAGTATCGATATCTTTTAACGTATCTGTTgcgg              >  2:833732/1‑99 (MQ=255)
                                                                                        gaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt            >  1:1760599/1‑100 (MQ=255)
                                                                                        gaagTGAGTTTTAACGCACTGCTTATCTAGACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGAGAt            >  1:1584793/1‑100 (MQ=255)
                                                                                         aagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATa           <  1:1887285/100‑1 (MQ=255)
                                                                                           gTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACt         >  1:2578438/1‑100 (MQ=255)
                                                                                            tGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGTATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACt         <  2:1595901/99‑1 (MQ=255)
                                                                                            tGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTa        >  1:2285875/1‑100 (MQ=255)
                                                                                               gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt     <  1:2420945/100‑1 (MQ=255)
                                                                                               gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt     >  2:456842/1‑100 (MQ=255)
                                                                                               gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt     >  2:3081279/1‑100 (MQ=255)
                                                                                                 tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTTAATTAGTTCAAGTGTCAGTATCGATCTCTTTTACCGTATATGTTGGGTATACTATTTg    >  2:2212118/1‑99 (MQ=255)
                                                                                                  aaatCTCACTTCTTATCTATACTACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCCTATATGTTGCGGATACTATTTTta  <  2:3393879/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
ATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTA  >  NC_000913/3798136‑3798333

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