Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,798,234 A→T 19.4% intergenic (+27/+5) 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,2340AT19.4% 473.9 / 109.3 212intergenic (+27/+5)waaL/waaUO‑antigen ligase/lipopolysaccharide core biosynthesis
Reads supporting (aligned to +/- strand):  ref base A (83/88);  new base T (25/16);  total (108/104)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

Only 100 of 220 total aligned reads displayed.

CATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTA  >  NC_000913/3798135‑3798333
                                                                                                   |                                                                                                   
caTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTa                                                                                                     <  1:109289/100‑1 (MQ=255)
caTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTa                                                                                                     <  2:935745/100‑1 (MQ=255)
 aTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaa                                                                                                    <  2:492888/100‑1 (MQ=255)
  ttCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                   <  2:3995110/100‑1 (MQ=255)
  ttCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                   >  2:1102765/1‑100 (MQ=255)
    ccAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAc                                                                                                 <  2:3923890/100‑1 (MQ=255)
     cAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACt                                                                                                <  2:407308/100‑1 (MQ=255)
     cAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTGTAAGAAGTGAGTTAAAACt                                                                                                <  1:4201998/100‑1 (MQ=255)
       aTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                              <  2:2023949/100‑1 (MQ=255)
         tATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAct                                                                                            <  1:394526/100‑1 (MQ=255)
         tATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAct                                                                                            <  2:1021496/100‑1 (MQ=255)
           tCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtc                                                                                          >  2:3660697/1‑100 (MQ=255)
           tCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtc                                                                                          <  2:3741712/100‑1 (MQ=255)
            cATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAATCTCActtct                                                                                         <  1:2011664/100‑1 (MQ=255)
            cATTATATCCGCTATAGTCATCTTACGCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActtct                                                                                         <  1:2280856/100‑1 (MQ=39)
             aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtctt                                                                                        >  2:941535/1‑100 (MQ=255)
             aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtctt                                                                                        <  1:618007/100‑1 (MQ=255)
             aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtctt                                                                                        >  1:3103598/1‑100 (MQ=255)
             aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTAActtctt                                                                                        >  1:3228514/1‑100 (MQ=255)
               tataTCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaa                                                                                      <  2:3155627/100‑2 (MQ=255)
                atatCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATc                                                                                     <  2:2130322/100‑1 (MQ=255)
                atatCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATc                                                                                     <  2:3194786/100‑1 (MQ=255)
                  atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaact                                                                                   <  2:3406296/100‑5 (MQ=38)
                  atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaact                                                                                   <  2:2857238/100‑5 (MQ=38)
                  atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCta                                                                                   >  1:2510178/1‑100 (MQ=255)
                  atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCta                                                                                   >  2:3800145/1‑100 (MQ=255)
                    ccGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtata                                                                                 <  2:1482718/100‑1 (MQ=255)
                    ccGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtata                                                                                 >  2:2530306/1‑100 (MQ=255)
                    ccGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtata                                                                                 >  1:3513278/1‑100 (MQ=255)
                    ccGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtata                                                                                 <  2:3072420/100‑1 (MQ=255)
                    ccGATATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaacttg                                                                                 <  1:4090211/100‑7 (MQ=37)
                     cGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaacttg                                                                                 <  1:1649841/99‑7 (MQ=37)
                      gCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACa                                                                               >  1:1563367/1‑100 (MQ=255)
                      gCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACa                                                                               <  2:564745/100‑1 (MQ=255)
                        tttaGTCCTCTTACTCGTCATTAATAATTGTAGCAATATAATTAATTAAGAATAAACAAGTATAAGAAGTGAGTTTTAACTCACTCCTTaaa                                                                                     <  1:438621/90‑3 (MQ=25)
                        tataGTCCTCTTACACATCATTAATAATCGTAACAATGCAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACaa                                                                              <  1:442088/99‑1 (MQ=255)
                         ataGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACt                                                                            <  2:165162/100‑1 (MQ=255)
                            gTCCTCTTACTCGTCAGTAATAATCGTAACAATTCAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTaa                                                                         >  1:2849183/1‑100 (MQ=255)
                              ccTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTACCTATACAACTTAAtc                                                                       <  2:1289221/100‑1 (MQ=255)
                                  ttACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCtt                                                                   <  1:941531/100‑1 (MQ=255)
                                   tACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                  <  2:173807/100‑1 (MQ=255)
                                     cTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTAc                                                                >  1:554681/1‑100 (MQ=255)
                                     cTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTAc                                                                >  2:2351418/1‑100 (MQ=255)
                                       cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACCACTTAATCTCTTTACTg                                                              >  1:2209868/1‑100 (MQ=255)
                                        gTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGa                                                             >  2:138941/1‑100 (MQ=255)
                                        gTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTAAGTTAAAACAAACTTCTTATCTATAAAAATTAATCTCTTTa                                                                 >  1:1075192/1‑96 (MQ=255)
                                         tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGaa                                                            >  1:2789682/1‑100 (MQ=255)
                                          cATTAATAATCGTAACAATACAATTAATTAAGAATAAACAATTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt                                                           >  1:3673702/1‑100 (MQ=255)
                                          cATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACGACTTAATCTCTTTACTGAAt                                                           <  1:2683785/100‑1 (MQ=255)
                                            ttAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTa                                                         <  1:2739273/100‑1 (MQ=255)
                                              aataatCTTAACTATACAATTAATTAAGAATAATCAGGTTTAAGCATTGATTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTa                                                         <  1:1412918/98‑1 (MQ=255)
                                                 aatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaa                                                    >  2:1309192/1‑100 (MQ=255)
                                                  atCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaaa                                                   >  2:3246819/1‑100 (MQ=255)
                                                  atCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaaa                                                   <  2:2764412/100‑1 (MQ=255)
                                                   tCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAg                                                  <  1:674024/100‑1 (MQ=255)
                                                        acaatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCa                                             <  2:3645187/100‑1 (MQ=255)
                                                         caatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAg                                            >  2:1184577/1‑100 (MQ=255)
                                                         caatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAg                                            <  1:3623507/100‑1 (MQ=255)
                                                         caatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAg                                            >  2:209156/1‑100 (MQ=255)
                                                         caatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAg                                            >  2:2570617/1‑100 (MQ=255)
                                                           atacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGt                                           <  1:378314/99‑1 (MQ=255)
                                                            tacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATCAGTTAAAGTTTAAGTa                                          >  1:4105467/1‑99 (MQ=255)
                                                              caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                       >  2:2489492/1‑100 (MQ=255)
                                                               aattaattAAGAATAAACAAGTTTAAGAAGTGAGTTAAGACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGa                                      >  2:1624398/1‑100 (MQ=255)
                                                                   aattaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATc                                  <  2:1706074/100‑1 (MQ=255)
                                                                    attaaGAATAAACAAGTTTAAGAAGTGAGTTACAACTCACTTATTATATATACAACTTAATCTCTTTACCAAATTAGTTAAATTTTCAGTATCGatatgt                                 >  1:859032/1‑98 (MQ=255)
                                                                      taaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTAACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATAGATATCttt                               >  1:3708514/1‑100 (MQ=255)
                                                                        aGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTa                             <  1:496482/100‑1 (MQ=255)
                                                                          aaTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAATTTTCAGTATCGATATCTTTTAc                            <  1:2553760/99‑1 (MQ=255)
                                                                              tgtatagatAAGAAGTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtat                       >  2:343980/9‑100 (MQ=37)
                                                                               aCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtata                      <  1:1841986/100‑1 (MQ=255)
                                                                                  aGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGt                   >  2:209849/1‑100 (MQ=255)
                                                                                   gTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtt                  >  1:2927675/1‑100 (MQ=255)
                                                                                   gTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtt                  >  2:1351421/1‑100 (MQ=255)
                                                                                   gTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtt                  >  1:3839338/1‑100 (MQ=255)
                                                                                   gTTTAAGAAGTGAGTTAAAAAACAATTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTAAGTATCGATATCTTTTACCGTATATgat                  >  1:4120495/1‑98 (MQ=255)
                                                                                      tAAGAAGTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg               >  2:635573/1‑100 (MQ=255)
                                                                                       aagaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg              >  2:2192479/1‑100 (MQ=255)
                                                                                       aagaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg              >  2:3588382/1‑100 (MQ=255)
                                                                                       aagaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg              >  2:729331/1‑100 (MQ=255)
                                                                                        agaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg              >  2:872752/1‑99 (MQ=255)
                                                                                        agaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg              >  2:3522399/1‑99 (MQ=255)
                                                                                        agaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg              >  2:417624/1‑99 (MQ=255)
                                                                                        agaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa             >  2:2240587/1‑100 (MQ=255)
                                                                                        agaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGATGGGGa             >  1:1323396/1‑100 (MQ=255)
                                                                                         gaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt            >  1:298539/1‑100 (MQ=255)
                                                                                         gaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACAGTATATGTTGGGGAt            >  1:4056794/1‑100 (MQ=255)
                                                                                         gaagTGAGTTTTAACTCACGTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt            >  1:2439624/1‑100 (MQ=255)
                                                                                          aagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACTGTATATGTTGGGGATa           >  1:2139050/1‑100 (MQ=255)
                                                                                           agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTTGGGATAc          <  2:2302455/100‑1 (MQ=255)
                                                                                           agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATAc          <  1:907823/100‑1 (MQ=255)
                                                                                            gTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACt         >  1:3528202/1‑100 (MQ=255)
                                                                                             tGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACt         <  1:3278018/99‑1 (MQ=255)
                                                                                              gAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAt       <  1:1490784/100‑1 (MQ=255)
                                                                                               aGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt      <  2:3316511/100‑1 (MQ=255)
                                                                                               aGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt      >  2:2623759/1‑100 (MQ=255)
                                                                                               aGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt      <  2:2212296/100‑1 (MQ=255)
                                                                                                gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt     >  2:3757278/1‑100 (MQ=255)
                                                                                                 ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTg    <  1:3660675/100‑1 (MQ=255)
                                                                                                  tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTATTACCGTATATGTTGGGGATACTATTTGt   >  1:2638094/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTA  >  NC_000913/3798135‑3798333

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: