Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,239,564 A→T 7.7% intergenic (+31/‑159) USA300HOU_RS06160 → / → recG hypothetical protein/ATP‑dependent DNA helicase RecG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,239,5640AT 7.7% 737.5 / 36.8 226intergenic (+31/‑159)USA300HOU_RS06160/recGhypothetical protein/ATP‑dependent DNA helicase RecG
Reads supporting (aligned to +/- strand):  ref base A (110/98);  new base T (9/9);  total (119/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 242 total aligned reads displayed.

GAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCAAC  >  USA300TCH1516_ALE/1239474‑1239654
                                                                                           |                                                                                          
gAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaaa                                                                                            >  2:3063259/1‑91 (MQ=255)
  tgTGGAAGTTCATGAAGATTGACAACCAATCTATCAATATTTCTTTTCAGAAAAACAAAAATTTAAAATAAATACATCCCCATGATAAatc                                                                                          >  2:3042023/2‑91 (MQ=255)
  aGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatc                                                                                          >  2:1616792/1‑91 (MQ=255)
  aGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTATTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATGAatc                                                                                          <  1:3070944/91‑1 (MQ=255)
   gTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatca                                                                                         <  2:1837036/91‑1 (MQ=255)
       aaGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAg                                                                                     >  2:3480544/1‑91 (MQ=255)
       aaGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAg                                                                                     >  2:1906368/1‑91 (MQ=255)
        aGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGt                                                                                    >  1:3171533/1‑91 (MQ=255)
        aGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGt                                                                                    >  2:2398602/1‑91 (MQ=255)
           tCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTgg                                                                                 <  1:2539066/91‑1 (MQ=255)
            cATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGt                                                                                >  2:2313709/1‑91 (MQ=255)
            cATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGt                                                                                >  2:2019848/1‑91 (MQ=255)
             aTGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTa                                                                               >  2:536903/1‑91 (MQ=255)
             aTGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTa                                                                               >  1:293052/1‑91 (MQ=255)
              tGAAGGTGGACAACAAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAg                                                                              <  2:1441271/91‑1 (MQ=255)
               gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                             >  1:3469647/1‑91 (MQ=255)
               gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                             <  1:1268199/91‑1 (MQ=255)
                aaGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtt                                                                            >  1:3225600/1‑91 (MQ=255)
                 aGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGttt                                                                           <  1:3694659/91‑1 (MQ=255)
                   ttGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGttttt                                                                         <  2:3853740/90‑1 (MQ=18)
                   gTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGttttt                                                                         <  2:321157/91‑1 (MQ=255)
                    tGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtttttt                                                                        >  2:2777059/1‑91 (MQ=255)
                      gACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAt                                                                      <  2:2632439/91‑1 (MQ=18)
                      gACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAt                                                                      >  1:2146653/1‑91 (MQ=255)
                       aCAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAtt                                                                     <  1:2681907/91‑1 (MQ=18)
                       aCAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAtt                                                                     <  1:1072607/91‑1 (MQ=17)
                       aCAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAtt                                                                     <  1:576772/91‑1 (MQ=255)
                         aaCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAtttt                                                                   <  2:1660879/91‑1 (MQ=17)
                         aaCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAtttt                                                                   >  1:442885/1‑91 (MQ=255)
                           ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                                 <  1:1959414/91‑1 (MQ=255)
                           ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                                 <  2:2001241/91‑1 (MQ=255)
                           ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                                 <  1:968107/91‑1 (MQ=255)
                            cAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAttttaaa                                                                <  1:2697882/91‑4 (MQ=14)
                               tttATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTAtt                                                             <  2:3276818/91‑1 (MQ=255)
                                    cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                        >  2:2636004/1‑91 (MQ=255)
                                    cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                        >  2:3376188/1‑91 (MQ=255)
                                    cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                        <  2:826824/91‑1 (MQ=255)
                                     aaTATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTg                                                       <  2:436151/91‑1 (MQ=255)
                                      atatTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGa                                                      <  1:163005/91‑1 (MQ=255)
                                        atTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGata                                                    >  1:904673/1‑91 (MQ=255)
                                         tttCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGatat                                                   <  2:274377/91‑1 (MQ=255)
                                         tttCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGatat                                                   <  2:273440/91‑1 (MQ=255)
                                           tCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTg                                                 >  2:3492851/1‑91 (MQ=255)
                                           tCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTg                                                 >  2:2117992/1‑91 (MQ=255)
                                            cTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGc                                                >  1:3783622/1‑91 (MQ=255)
                                              tttCAGTAGAATAAAAGTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGctg                                              <  2:3712769/91‑3 (MQ=255)
                                              tttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCgg                                              >  2:379735/1‑91 (MQ=255)
                                              tttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCgg                                              >  2:477638/1‑91 (MQ=255)
                                               ttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCggg                                             <  2:2582502/91‑1 (MQ=255)
                                               ttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCggg                                             <  2:2585536/91‑1 (MQ=255)
                                               ttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCggg                                             >  1:3696674/1‑91 (MQ=255)
                                               ttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCggg                                             <  1:2777059/91‑1 (MQ=255)
                                                tCAGTAGAATAGAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGt                                            <  1:3513735/91‑1 (MQ=255)
                                                  aGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTa                                          >  2:2488359/1‑91 (MQ=255)
                                                    tAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTaaa                                        <  1:3855038/91‑1 (MQ=255)
                                                     aGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTaaaa                                       <  2:2316014/91‑1 (MQ=255)
                                                      gAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAg                                      <  1:3376188/91‑1 (MQ=255)
                                                      gAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAg                                      <  2:1288235/91‑1 (MQ=255)
                                                       aaTAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGt                                     <  1:512276/91‑1 (MQ=255)
                                                       aaTAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGt                                     <  1:1691001/91‑1 (MQ=255)
                                                          aaaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATc                                  >  1:2210410/1‑91 (MQ=255)
                                                          aaaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATc                                  <  2:119237/91‑1 (MQ=255)
                                                           aaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCg                                 >  1:341898/1‑91 (MQ=255)
                                                             aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                               <  1:983487/91‑1 (MQ=255)
                                                             aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                               >  2:926744/1‑91 (MQ=255)
                                                             aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                               <  1:3456434/91‑1 (MQ=255)
                                                               tttAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGctct                             >  1:3123436/1‑91 (MQ=255)
                                                               tttAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTtct                             >  2:3107932/1‑91 (MQ=255)
                                                                 taaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCg                           >  1:554669/1‑91 (MQ=255)
                                                                  aaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTGTTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGa                          <  1:3271713/91‑1 (MQ=255)
                                                                  aaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGa                          >  1:1578494/1‑91 (MQ=255)
                                                                    aataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGt                        >  2:1758932/1‑91 (MQ=255)
                                                                          aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                  >  1:2024157/1‑91 (MQ=18)
                                                                          aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                  >  1:3143115/1‑91 (MQ=18)
                                                                          aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                  >  2:1879964/1‑91 (MQ=18)
                                                                          aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTGGCTaa                  >  2:3581352/1‑91 (MQ=18)
                                                                          aaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                  <  1:1875001/91‑1 (MQ=255)
                                                                          aaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                  <  2:86586/91‑1 (MQ=255)
                                                                           aaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAAc                 >  1:761564/1‑91 (MQ=255)
                                                                            aCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACa                >  2:115424/1‑91 (MQ=255)
                                                                            aCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACa                <  2:2936261/91‑1 (MQ=255)
                                                                            aCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACa                <  1:2769600/91‑1 (MQ=255)
                                                                             cTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACaa               >  1:3253683/1‑91 (MQ=255)
                                                                             cTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACaa               >  2:1731621/1‑91 (MQ=255)
                                                                             cTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACaa               >  2:164814/1‑91 (MQ=255)
                                                                             cTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACaa               >  1:710420/1‑91 (MQ=255)
                                                                              tgccAATGATAAATCATCAGTTGGGAGTTTTTTATTTTGCTATTTTAGTGTTATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAAt              <  2:3599512/89‑1 (MQ=255)
                                                                              tACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAAt              <  2:2265623/91‑1 (MQ=255)
                                                                               aCCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATg             >  1:64980/1‑91 (MQ=255)
                                                                                ccAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGt            <  1:3151271/91‑1 (MQ=255)
                                                                                ccAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGt            <  1:1828550/91‑1 (MQ=255)
                                                                                   aTGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCAt         >  1:2500144/1‑91 (MQ=255)
                                                                                     tGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATg        <  1:423071/91‑1 (MQ=255)
                                                                                      gATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGt       >  1:2294076/1‑91 (MQ=255)
                                                                                       aTAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGtt      <  2:1380683/91‑1 (MQ=255)
                                                                                       aTAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGtt      <  2:3011240/91‑1 (MQ=255)
                                                                                        tAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTc     <  1:1714955/91‑1 (MQ=255)
                                                                                        tAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTc     <  1:1389930/91‑1 (MQ=255)
                                                                                          aaTCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCaa   <  1:1188027/91‑1 (MQ=255)
                                                                                           atcatcAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCAAc  >  2:2008902/1‑91 (MQ=255)
                                                                                           |                                                                                          
GAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCAAC  >  USA300TCH1516_ALE/1239474‑1239654

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