Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,798,235 A→T 17.4% 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,2350AT17.4% 507.1 / 95.4 214intergenic (+28/+4)waaL/waaUO‑antigen ligase/lipopolysaccharide core biosynthesis
Reads supporting (aligned to +/- strand):  ref base A (85/90);  new base T (16/21);  total (102/112)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

Only 100 of 223 total aligned reads displayed.

ATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTAT  >  NC_000913/3798136‑3798334
                                                                                                   |                                                                                                   
aTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaa                                                                                                     >  1:2660318/1‑100 (MQ=255)
   ccAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAc                                                                                                  <  1:552933/100‑1 (MQ=255)
   ccAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAc                                                                                                  <  2:1939815/100‑1 (MQ=255)
    cAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACt                                                                                                 >  1:2276191/1‑100 (MQ=255)
     aaTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTc                                                                                                <  2:2093394/100‑1 (MQ=255)
     aaTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTc                                                                                                >  1:2122536/1‑100 (MQ=255)
     aaTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTc                                                                                                >  1:1851744/1‑100 (MQ=255)
         aTCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActt                                                                                            >  2:350657/1‑100 (MQ=255)
          tCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtc                                                                                           <  1:2811618/100‑1 (MQ=255)
          tCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtc                                                                                           <  1:2040633/100‑1 (MQ=255)
          tCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtc                                                                                           <  1:3048771/100‑1 (MQ=255)
           cATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActtct                                                                                          <  1:2871616/100‑1 (MQ=255)
            gttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCTTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActtctt                                                                                         <  2:3891497/99‑1 (MQ=39)
            aTTATATCCGCTATATTCCTTTTACTCGTCATTATTAATCGTAACAATACATTTAATTAAGTTTAAACAAGTTTAAGAAGTAATTTataataaaaaac                                                                                           >  2:150320/1‑90 (MQ=255)
            aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATTAACAAGTTTAAGAAGTGAGTTAAAACTCActtct                                                                                          >  2:3926092/1‑99 (MQ=255)
            aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtctt                                                                                         >  1:4133634/1‑100 (MQ=255)
             ttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActtct                                                                                          <  2:445265/98‑1 (MQ=255)
             ttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTa                                                                                        <  1:474412/100‑1 (MQ=255)
             ttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTa                                                                                        <  1:3601003/100‑1 (MQ=255)
              tatatcatatataGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTAt                                                                                       <  2:3574166/91‑1 (MQ=255)
              tataTCCGCTATAGTCTTCTTACTCGTCATTAATAATCGTAAAAATTCAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaa                                                                                       <  2:3651873/100‑2 (MQ=38)
              tataTCCGCTATAGTCCTTTTACTCGTAATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaa                                                                                       <  2:3522090/100‑2 (MQ=38)
              tataTCCGCTATAGTCCTATACTTGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaa                                                                                      <  2:1998977/100‑3 (MQ=37)
               atatCCGCTATAGTCCTTTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaa                                                                                       <  2:3642862/99‑2 (MQ=38)
               atatCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATc                                                                                      <  1:2205788/100‑1 (MQ=255)
               atatCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATc                                                                                      <  1:1546784/100‑1 (MQ=255)
               atatCCGCTATAGTCCTATTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAGACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaa                                                                                       <  1:2904722/99‑2 (MQ=38)
                 atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaact                                                                                    <  2:2846357/100‑5 (MQ=38)
                 atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCta                                                                                    <  1:3283274/100‑1 (MQ=255)
                 atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCta                                                                                    <  2:1880828/100‑1 (MQ=255)
                  tCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtat                                                                                   >  2:3834563/1‑100 (MQ=255)
                  tCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtat                                                                                   >  1:1678215/1‑100 (MQ=255)
                  tCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtat                                                                                   <  1:1523745/100‑1 (MQ=255)
                   ccGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtata                                                                                  >  2:601807/1‑100 (MQ=255)
                   ccGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAATCACTTCTTATCtata                                                                                  >  2:2412369/1‑100 (MQ=255)
                     gCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaacttgtt                                                                                <  1:3802040/100‑9 (MQ=37)
                        ataGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACt                                                                             <  2:854998/100‑1 (MQ=255)
                           gatctctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTaa                                                                          <  2:1662284/97‑1 (MQ=255)
                              ctctTTCTCGTCATTAATAATCGTAACAATACATTTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAt                                                                         <  2:1078139/98‑1 (MQ=255)
                                ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTaaagatc                                                                     >  1:1662296/1‑95 (MQ=255)
                                ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATctct                                                                     >  1:2113783/1‑100 (MQ=255)
                                ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATctct                                                                     <  2:671113/100‑1 (MQ=255)
                                 ttACTCGTCATTAATAATCGTAACAATACTATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTTATctct                                                                     >  2:562200/1‑99 (MQ=255)
                                 ttACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGATTTAAAGCTCACTTCTTATCTATACAACTTAATCTCtt                                                                    <  2:1539723/100‑1 (MQ=255)
                                  tACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                   >  2:3300577/1‑100 (MQ=255)
                                  tACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                   <  2:1567898/100‑1 (MQ=255)
                                  tACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCtt                                                                    >  2:4019542/1‑99 (MQ=255)
                                    cTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTAc                                                                 >  1:4006574/1‑100 (MQ=255)
                                     tCGTAATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACt                                                                <  2:1216797/100‑1 (MQ=255)
                                      cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTg                                                               >  2:2887633/1‑100 (MQ=255)
                                      cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTg                                                               >  2:1124605/1‑100 (MQ=255)
                                         cATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt                                                            >  1:1906627/1‑100 (MQ=255)
                                         cATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt                                                            >  1:241360/1‑100 (MQ=255)
                                         cATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt                                                            >  1:429800/1‑100 (MQ=255)
                                         cATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt                                                            >  1:469635/1‑100 (MQ=255)
                                           ttAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTa                                                          <  1:715285/100‑1 (MQ=255)
                                            taataaTCGTAACAATACAATTAATTAAGAATTAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAg                                                         <  2:1533679/100‑1 (MQ=255)
                                            taataaTCGTAACAATACAATTAATTAAGAATAAATAAGTTTAATAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAg                                                         <  2:1718013/100‑1 (MQ=255)
                                             aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGt                                                        >  1:3184035/1‑100 (MQ=255)
                                             aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGt                                                        >  1:3434816/1‑100 (MQ=255)
                                                  tCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAg                                                   >  2:302261/1‑100 (MQ=255)
                                                     tAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAATTAAAACTCACTTCTTATATATACAAATTAATCTATTTACTGAATTATTTAAAGtt                                                 >  2:1022373/1‑99 (MQ=255)
                                                      aaCAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                               <  2:3966941/100‑1 (MQ=255)
                                                        caatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAg                                             <  1:1740827/100‑1 (MQ=255)
                                                          atacaatTAATTAAGAATAAACAAGTTTAAGAACCTAGTTAAAACTCTCTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGt                                            <  2:2240937/99‑1 (MQ=255)
                                                          atacaatTAATTAAAAATAAACAAGTTTAATAAGTGATTTATAACTCAATTTTTATCTATAGAACTTAATCTCTTTACTGAATGAGTTAAAGTTTCAGTa                                           <  2:2547790/100‑1 (MQ=255)
                                                            acaatTAATTAAGAATAAACAAGTTTAAGAAGTTAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATc                                         <  2:1672960/100‑1 (MQ=255)
                                                             caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                        >  1:3390572/1‑100 (MQ=255)
                                                             caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCGTTACTGAATTAGTTAAAGTTTCAGTATCg                                        >  2:2486774/1‑100 (MQ=255)
                                                              aattaattAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGa                                       <  2:3412139/100‑1 (MQ=255)
                                                                ttaattaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGata                                     <  1:3718348/100‑1 (MQ=255)
                                                                     taaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCttt                                >  1:1343060/1‑100 (MQ=255)
                                                                     taaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCttt                                <  1:757768/100‑1 (MQ=255)
                                                                      aaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCCGTATCGATATCtttt                               >  2:50626/1‑100 (MQ=255)
                                                                       aGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTa                              <  2:4120751/100‑1 (MQ=255)
                                                                         aaTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTAcc                            >  1:4042096/1‑100 (MQ=255)
                                                                          aTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCg                           >  2:1773548/1‑100 (MQ=255)
                                                                              aCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtata                       >  1:1460269/1‑100 (MQ=255)
                                                                              aCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtata                       <  1:1007598/100‑1 (MQ=255)
                                                                                    atAAGAAGTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgg                 >  1:1089040/2‑100 (MQ=255)
                                                                                      aagaagTGAGTTAGCACTCACTTTTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAACAAACAGTATAGATATCTTTTACCGCATATGTTggg                >  2:3530872/1‑99 (MQ=255)
                                                                                      aagaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg               <  2:1542106/100‑1 (MQ=255)
                                                                                      aagaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCATATATGTTgggg               >  1:1196868/1‑100 (MQ=255)
                                                                                       agaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa              >  1:651410/1‑100 (MQ=255)
                                                                                       agaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa              >  2:2354382/1‑100 (MQ=255)
                                                                                       aaaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa              >  1:2804249/3‑100 (MQ=255)
                                                                                        gaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt             >  1:2524706/1‑100 (MQ=255)
                                                                                        gaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt             >  2:758842/1‑100 (MQ=255)
                                                                                        gaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt             <  2:2658494/100‑1 (MQ=255)
                                                                                        gaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt             >  1:3864428/1‑100 (MQ=255)
                                                                                          agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATAc           <  2:1851731/100‑1 (MQ=255)
                                                                                          agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATAc           <  1:125263/100‑1 (MQ=255)
                                                                                           gTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACt          >  1:1152258/1‑100 (MQ=255)
                                                                                            tGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTa         >  1:1830467/1‑100 (MQ=255)
                                                                                              aGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt       >  2:2750205/1‑100 (MQ=255)
                                                                                              aGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt       >  1:3462409/1‑100 (MQ=255)
                                                                                               gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATTTTTTACCGTATATGATGGGGATACTAttt      >  2:2592355/1‑100 (MQ=255)
                                                                                                ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTg     <  2:3681459/100‑1 (MQ=255)
                                                                                                ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTg     <  2:3184016/100‑1 (MQ=255)
                                                                                                   aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGtat  >  1:3165127/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTAT  >  NC_000913/3798136‑3798334

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: