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

Only 100 of 219 total aligned reads displayed.

ATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTAT  >  NC_000913/3798136‑3798334
                                                                                                   |                                                                                                   
aTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaa                                                                                                     >  2:536067/1‑100 (MQ=255)
 ttCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                    <  2:3755024/100‑1 (MQ=255)
      aTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                               <  2:2023949/100‑1 (MQ=255)
         aTCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActt                                                                                            >  2:3840580/1‑100 (MQ=255)
          tCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtc                                                                                           >  2:2550922/1‑100 (MQ=255)
           cATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActtct                                                                                          <  2:3145910/100‑1 (MQ=255)
           cATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtct                                                                                          <  2:3443429/100‑1 (MQ=255)
            aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtctt                                                                                         >  2:941535/1‑100 (MQ=255)
            aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtctt                                                                                         >  1:636367/1‑100 (MQ=255)
             ttATATTCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTa                                                                                        <  2:3672236/100‑1 (MQ=255)
             ttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACTAGTTTAAGAAGTGAGTTTTAACTCACTTCTTa                                                                                        <  2:4174311/100‑1 (MQ=39)
              tataTCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaa                                                                                       <  2:3170988/100‑2 (MQ=255)
              tataTCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTAt                                                                                       >  2:3204485/1‑100 (MQ=255)
                 atCCGCTATAGTCCTCTTACTCGTCATTAATAATCTTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaact                                                                                    <  1:177619/100‑5 (MQ=38)
                 atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaact                                                                                    <  2:2857238/100‑5 (MQ=38)
                 atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaact                                                                                    <  1:3260657/100‑5 (MQ=38)
                 atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCta                                                                                    >  1:1455616/1‑100 (MQ=255)
                   ccGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtata                                                                                  >  1:3513278/1‑100 (MQ=255)
                   ccGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtata                                                                                  <  2:1482718/100‑1 (MQ=255)
                   ccGATATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTaaacttg                                                                                  <  1:4090211/100‑7 (MQ=37)
                    cGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATAc                                                                                 <  2:761649/100‑1 (MQ=255)
                    cGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTAACTTCTTATCta                                                                                    >  1:3645165/1‑97 (MQ=255)
                     gCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACa                                                                                >  1:1563367/1‑100 (MQ=255)
                       tataGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAAc                                                                              <  2:2209443/100‑1 (MQ=255)
                       tataGTCCTCTTACACATCATTAATAATCGTAACAATGCAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACaa                                                                               <  1:442088/99‑1 (MQ=255)
                        ataGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACt                                                                             <  2:165162/100‑1 (MQ=255)
                        ataGTCATCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACt                                                                             <  1:2764179/100‑1 (MQ=255)
                        ataGACGTCTTACTCGTCATTAATAATCGTAACAATTCCATTAATTAAGAATAAACTAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACt                                                                             <  1:3312792/100‑1 (MQ=255)
                          aGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTa                                                                           >  2:1410001/1‑100 (MQ=255)
                           gTCCTCTTACTCGTCATTAATAATCGGAACAATTCAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACtt                                                                            <  1:2072860/98‑1 (MQ=255)
                             ccTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAtc                                                                        <  1:1208923/100‑1 (MQ=255)
                             ccTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTACCTATACAACTTAAtc                                                                        <  2:1289221/100‑1 (MQ=255)
                                 ttACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCtt                                                                    <  1:941531/100‑1 (MQ=255)
                                  tACTCGTCATTAATAATCGTAATAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                   <  1:1851603/100‑1 (MQ=255)
                                  tACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                   >  2:3392772/1‑100 (MQ=255)
                                  tACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                   >  1:1388310/1‑100 (MQ=255)
                                   aCTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTa                                                                  <  1:3195367/100‑1 (MQ=255)
                                    cTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTAc                                                                 >  1:554681/1‑100 (MQ=255)
                                      cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACCACTTAATCTCTTTACTg                                                               >  1:2209868/1‑100 (MQ=255)
                                      cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTg                                                               <  1:683026/100‑1 (MQ=255)
                                      cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTg                                                               >  1:4390534/1‑100 (MQ=255)
                                       gTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGa                                                              >  2:1857576/1‑100 (MQ=255)
                                       gTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTAAGTTAAAACAAACTTCTTATCTATAAAAATTAATCTCTTTa                                                                  >  1:1075192/1‑96 (MQ=255)
                                            taataaTCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAg                                                         <  2:224675/100‑1 (MQ=255)
                                             aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGt                                                        <  1:2803234/100‑1 (MQ=255)
                                                aatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAATTAACTTATTATCTATACAACTTAATCTATTTACTGAATTAGTTa                                                      >  1:2603671/1‑99 (MQ=255)
                                                aatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaa                                                     >  2:1309192/1‑100 (MQ=255)
                                                   cGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGt                                                  >  1:4055569/1‑100 (MQ=255)
                                                    gTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGtt                                                 <  2:3615329/100‑1 (MQ=255)
                                                      aaCAATACAATTAATTAAGGATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                               <  2:2098339/100‑1 (MQ=255)
                                                      aaCAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTAAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                               <  1:1496344/100‑1 (MQ=255)
                                                        caatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAg                                             >  2:668951/1‑100 (MQ=255)
                                                        caatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAg                                             <  2:608872/100‑1 (MQ=255)
                                                        caatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAg                                             >  2:1184577/1‑100 (MQ=255)
                                                        caatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAg                                             <  2:1706310/100‑1 (MQ=255)
                                                          atacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGt                                            <  1:378314/99‑1 (MQ=255)
                                                           tacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATCAGTTAAAGTTTAAGTa                                           >  1:4105467/1‑99 (MQ=255)
                                                             caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGGATTAGTTAAAGTTTCAGTATCg                                        <  1:2550909/100‑1 (MQ=255)
                                                             caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                        >  2:1380258/1‑100 (MQ=255)
                                                             caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                        >  2:210119/1‑100 (MQ=255)
                                                             caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                        >  2:2489492/1‑100 (MQ=255)
                                                              aattaattAAGAATAAACAAGTTTAAGAAGTGAGTTAAGACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGa                                       >  2:1624398/1‑100 (MQ=255)
                                                              aattaattAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCATATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTa                                           >  1:4378350/1‑96 (MQ=255)
                                                               attaattaAGAATAAACAAGTTTAAGACGTTAGTTAAAACTCACGTCATATCTATACAACTTAATCCCTTTACTGAATTAGTTAAAGTATCAGTATCg                                        <  1:2478571/98‑1 (MQ=255)
                                                                   attaaGAATAAACAAGTTTAAGAAGTGAGTTACAACTCACTTATTATATATACAACTTAATCTCTTTACCAAATTAGTTAAATTTTCAGTATCGatatgt                                  >  1:859032/1‑98 (MQ=255)
                                                                    ttaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCtt                                 <  2:3528095/100‑1 (MQ=255)
                                                                     taaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCttt                                <  1:4189184/100‑1 (MQ=255)
                                                                     taaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTAACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATAGATATCttt                                >  1:3708514/1‑100 (MQ=255)
                                                                      aaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCtttt                               <  2:3070620/100‑1 (MQ=255)
                                                                      aaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCttta                               >  2:4189214/1‑99 (MQ=255)
                                                                       aGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTa                              <  1:496482/100‑1 (MQ=255)
                                                                         aaTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCCTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTAcc                            <  1:2393605/100‑1 (MQ=255)
                                                                          aTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCg                           >  1:4073572/1‑100 (MQ=255)
                                                                           tAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGt                          >  2:1103982/1‑100 (MQ=255)
                                                                            aaaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGta                         >  2:3439157/1‑100 (MQ=255)
                                                                             aaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtat                        >  2:2644298/1‑100 (MQ=255)
                                                                                atagatAAGAAGTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTatat                      >  2:3956620/6‑99 (MQ=38)
                                                                                  gTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtt                   >  1:2927675/1‑100 (MQ=255)
                                                                                     tAAGAAGTGAGTTTTAACTCGCTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg                >  2:2584721/1‑100 (MQ=255)
                                                                                     tAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg                <  2:4055598/100‑1 (MQ=255)
                                                                                      aagaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg               >  2:2192479/1‑100 (MQ=255)
                                                                                       agaagTGAGTTTTCACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa              >  1:2719717/1‑100 (MQ=255)
                                                                                       agaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg               >  2:3657978/1‑99 (MQ=255)
                                                                                       agaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa              >  2:2507929/1‑100 (MQ=255)
                                                                                       agaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGa              >  2:23508/1‑100 (MQ=255)
                                                                                        gaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt             >  2:3787490/1‑100 (MQ=255)
                                                                                        gaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt             >  1:298539/1‑100 (MQ=255)
                                                                                        gaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt             <  1:3840555/100‑1 (MQ=255)
                                                                                          agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATAc           >  1:2141575/1‑100 (MQ=255)
                                                                                          agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATAc           <  1:907823/100‑1 (MQ=255)
                                                                                          agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATAc           <  2:883882/100‑1 (MQ=255)
                                                                                           gTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACt          >  1:3528202/1‑100 (MQ=255)
                                                                                            taatttAAAACTCACGTCTTATCTAGACAACTGAAGCTCTTTACTGCATGAGTTAAAGTTTCATTATCGTTATCTTTGACCGTATATGTTGGGGATACTa         >  1:3673659/5‑100 (MQ=255)
                                                                                            tGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTa         >  1:2866569/1‑100 (MQ=255)
                                                                                              aGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt       <  2:3316511/100‑1 (MQ=255)
                                                                                               gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt      >  1:368695/1‑100 (MQ=255)
                                                                                               gTTAAAACTCACCTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt      >  2:1067817/1‑100 (MQ=255)
                                                                                                ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTg     <  2:2238537/100‑1 (MQ=255)
                                                                                                 tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTATTACCGTATATGTTGGGGATACTATTTGt    >  1:2638094/1‑100 (MQ=255)
                                                                                                   aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGtat  <  2:2063259/100‑1 (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: