Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,798,235 A→T 21.6% 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,2350AT21.6% 392.4 / 116.1 195intergenic (+28/+4)waaL/waaUO‑antigen ligase/lipopolysaccharide core biosynthesis
Reads supporting (aligned to +/- strand):  ref base A (80/72);  new base T (24/18);  total (104/91)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-03

Only 100 of 201 total aligned reads displayed.

ATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTA  >  NC_000913/3798136‑3798333
                                                                                                   |                                                                                                  
aTTCCAATTATCGTTATATCCGCTATAGTCCTCTTACTCGTCATCAAGAATAGTAAGAATACAATTAATTAAGAATAAACAAGTTTAGGACGTGAGTTaa                                                                                                    >  1:4178634/1‑100 (MQ=255)
aTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaa                                                                                                    <  2:2803069/100‑1 (MQ=255)
 ttCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                   <  2:1436976/100‑1 (MQ=255)
    cAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACt                                                                                                <  2:958841/100‑1 (MQ=255)
        tATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCAct                                                                                            <  1:2363497/100‑1 (MQ=255)
        tATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAct                                                                                            <  1:201065/100‑1 (MQ=255)
        tATCATTATATCCGCTATAGTCCCCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAc                                                                                             >  2:80335/1‑99 (MQ=255)
         aTCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActt                                                                                           <  2:252403/100‑1 (MQ=255)
          tCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtc                                                                                          >  1:473542/1‑100 (MQ=255)
           cATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActtct                                                                                         <  1:3840081/100‑1 (MQ=255)
           cATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActtct                                                                                         <  1:1083321/100‑1 (MQ=255)
            aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtctt                                                                                        >  2:2365461/1‑100 (MQ=255)
            aTTATATCCGATATAGTCCTCTTACTCGTCATTAAAACTCGAAAAAAAACACTTAATTAAAAATAAACAAGTTTAAGAAGTAAATTACAACtacctta                                                                                          <  2:346607/98‑7 (MQ=255)
             ttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTa                                                                                       <  1:833230/100‑1 (MQ=255)
             ttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTa                                                                                       <  1:1400577/100‑1 (MQ=255)
              tataTCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaa                                                                                      <  1:4139568/100‑2 (MQ=255)
              tataTCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaa                                                                                      <  2:3050054/100‑2 (MQ=255)
              tataTCCGATATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAATACTCActtct                                                                                         <  2:2952777/97‑1 (MQ=255)
               atatCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATc                                                                                     <  1:3972662/100‑1 (MQ=255)
               atatCCGCTATAGTCCTCTAACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaa                                                                                     <  2:1389280/100‑3 (MQ=38)
                 atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCta                                                                                   >  2:4106522/1‑100 (MQ=255)
                 atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCta                                                                                   >  1:1768746/1‑100 (MQ=255)
                 atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCta                                                                                   <  1:1639773/100‑1 (MQ=255)
                  tCCGATAAAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCta                                                                                   <  2:2600309/99‑1 (MQ=255)
                   ccGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtata                                                                                 <  1:1526431/100‑1 (MQ=255)
                   ccGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtata                                                                                 >  2:1715111/1‑100 (MQ=255)
                     gCTATAGTCCTCTTAATCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACa                                                                               <  2:1547860/100‑1 (MQ=255)
                      cTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaacttgttt                                                                              <  2:1153832/100‑10 (MQ=37)
                       tataGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAAc                                                                             >  1:2169946/1‑100 (MQ=255)
                       tataGTCCTCTGAGTCGTCATTAAGAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAAc                                                                             <  2:1494577/100‑1 (MQ=255)
                           gTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTaa                                                                         >  2:1550853/1‑100 (MQ=255)
                            tCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAt                                                                        <  2:4124746/100‑1 (MQ=255)
                                ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATctct                                                                    <  2:1742647/100‑1 (MQ=255)
                                 ttACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCtt                                                                   <  1:1295842/100‑1 (MQ=255)
                                 ttACTCATCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCtt                                                                   <  2:362775/100‑1 (MQ=255)
                                   aCTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTa                                                                 <  1:3822126/100‑1 (MQ=255)
                                    cTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTAc                                                                >  1:2789817/1‑100 (MQ=255)
                                      cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTg                                                              >  1:3237684/1‑100 (MQ=255)
                                      cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTg                                                              <  1:3022744/100‑1 (MQ=255)
                                      cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTg                                                              >  1:2946/1‑100 (MQ=255)
                                      cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTg                                                              >  2:1094519/1‑100 (MQ=255)
                                        tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGaa                                                            <  2:2199567/100‑1 (MQ=255)
                                        tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGaa                                                            <  1:1816840/100‑1 (MQ=255)
                                         cATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTATTCTCTTTACTGAAt                                                           >  2:944528/1‑100 (MQ=255)
                                           ttAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTa                                                         >  2:3759820/1‑100 (MQ=255)
                                            taataaTCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAg                                                        >  1:911290/1‑100 (MQ=255)
                                             aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTa                                                         >  2:2736331/1‑98 (MQ=255)
                                             aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGt                                                       >  1:261997/1‑100 (MQ=255)
                                                aatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaa                                                    >  2:1994888/1‑100 (MQ=255)
                                                    gTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGtt                                                <  1:3895256/100‑1 (MQ=255)
                                                     tAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGttt                                               >  1:1583465/1‑100 (MQ=255)
                                                      aaCAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                              <  2:3703857/100‑1 (MQ=255)
                                                        caatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAg                                            <  1:2321902/100‑1 (MQ=255)
                                                         aatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGt                                           >  1:3410790/1‑100 (MQ=255)
                                                           tacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTAt                                         >  1:4025039/1‑100 (MQ=255)
                                                             caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                       >  2:569340/1‑100 (MQ=255)
                                                             caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                       <  2:3084324/100‑1 (MQ=255)
                                                              aattaattAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGa                                      <  2:2128022/100‑1 (MQ=255)
                                                               attaattaAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGat                                     >  1:1432240/1‑100 (MQ=255)
                                                               attaattaAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGat                                     >  1:3441948/1‑100 (MQ=255)
                                                                       aGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTa                             >  1:2371835/1‑100 (MQ=255)
                                                                       aGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTa                             >  1:2487290/1‑100 (MQ=255)
                                                                        gAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTAc                            <  2:2371830/100‑1 (MQ=255)
                                                                         aaTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTAcc                           <  2:3868164/100‑1 (MQ=255)
                                                                          aTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCg                          <  1:2060237/100‑1 (MQ=255)
                                                                          aTAAAAAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATTTATACAACTTAATCTCGTTACTTAATTAGTTAAAGTTTAAGTAACGATATCTTTTACCg                          <  2:2387668/100‑1 (MQ=255)
                                                                           gttgtatagatAAGAAGTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGt                         >  2:3280562/11‑100 (MQ=25)
                                                                            aaaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTCAAAGTTTCAGTATCGATATCTTTTACCGta                        <  2:2554919/100‑1 (MQ=255)
                                                                             aaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtat                       >  2:1561985/1‑100 (MQ=255)
                                                                             aaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtat                       <  1:3135364/100‑1 (MQ=255)
                                                                             aaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtat                       >  1:597905/1‑100 (MQ=255)
                                                                              aCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGATAAAGTTTCAGTATCGATATCTTTTACCGtata                      <  2:4154638/100‑1 (MQ=255)
                                                                               cAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTatat                     <  1:2974428/100‑1 (MQ=255)
                                                                               cAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTatat                     >  1:3688542/1‑100 (MQ=255)
                                                                                  gcttAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtt                  >  1:1519835/3‑100 (MQ=255)
                                                                                   gatAAGAAGTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                 >  1:2433504/3‑100 (MQ=255)
                                                                                   gatAAGAAGTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                 >  1:3715215/3‑100 (MQ=255)
                                                                                     tAAGAAGTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                 >  2:102949/1‑98 (MQ=255)
                                                                                     tAAGAAGTGAGTTAAAACTCGCTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg               <  2:3435908/100‑1 (MQ=255)
                                                                                     tAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg               <  2:1264141/100‑1 (MQ=255)
                                                                                     tAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg               >  2:3963500/1‑100 (MQ=255)
                                                                                       agaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa             >  1:2081065/1‑100 (MQ=255)
                                                                                       agaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa             >  1:256768/1‑100 (MQ=255)
                                                                                       agaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa             >  2:33451/1‑100 (MQ=255)
                                                                                       agaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa             >  1:1212077/1‑100 (MQ=255)
                                                                                        gaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt            >  1:551894/1‑100 (MQ=255)
                                                                                        gaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt            >  1:3216786/1‑100 (MQ=255)
                                                                                        gaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt            >  2:3098893/1‑100 (MQ=255)
                                                                                        gaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt            >  2:3243005/1‑100 (MQ=255)
                                                                                        gaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt            >  1:3334304/1‑100 (MQ=255)
                                                                                          agTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATAc          >  1:1853793/1‑100 (MQ=255)
                                                                                           gTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACt         >  1:726720/1‑100 (MQ=255)
                                                                                              aGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt      >  2:2474351/1‑100 (MQ=255)
                                                                                               gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt     >  1:1757278/1‑100 (MQ=255)
                                                                                               gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt     >  2:3417429/1‑100 (MQ=255)
                                                                                               gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt     >  2:733082/1‑100 (MQ=255)
                                                                                                ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTg    <  1:301615/100‑1 (MQ=255)
                                                                                                ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTg    <  1:1722021/100‑1 (MQ=255)
                                                                                                 tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGt   >  2:2014032/1‑100 (MQ=255)
                                                                                                 tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGt   >  2:2992100/1‑100 (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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: