Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,239,564 A→T 5.9% 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 5.9% 474.7 / 16.1 136intergenic (+31/‑159)USA300HOU_RS06160/recGhypothetical protein/ATP‑dependent DNA helicase RecG
Reads supporting (aligned to +/- strand):  ref base A (86/42);  new base T (5/3);  total (91/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

Only 100 of 147 total aligned reads displayed.

GAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCAAC  >  USA300TCH1516_ALE/1239474‑1239654
                                                                                           |                                                                                          
gAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaaa                                                                                            >  1:3475587/1‑91 (MQ=255)
gAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaaa                                                                                            >  1:712511/1‑91 (MQ=255)
gAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaaa                                                                                            >  1:792997/1‑91 (MQ=255)
     ggAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatcatc                                                                                       <  1:1186354/91‑1 (MQ=255)
     ggAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatcatc                                                                                       >  1:1626649/1‑91 (MQ=255)
       aaGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAg                                                                                     >  1:1678458/1‑91 (MQ=255)
        aGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGt                                                                                    <  1:3602796/91‑1 (MQ=255)
         gTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGtt                                                                                   <  1:1265527/91‑1 (MQ=255)
         gTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGtt                                                                                   <  1:2881595/91‑1 (MQ=255)
          ttCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTg                                                                                  >  1:1719081/1‑91 (MQ=255)
           tCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTgg                                                                                 >  1:1315925/1‑91 (MQ=255)
            cATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGt                                                                                >  1:2413146/1‑91 (MQ=255)
            cATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGt                                                                                >  1:2413178/1‑91 (MQ=255)
             aTGAAGGTGAACAACCAATTTATCAATATTTCTTTTCAGTACAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTa                                                                               >  1:3592109/1‑91 (MQ=255)
              tGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAg                                                                              >  1:704084/1‑91 (MQ=255)
               gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                             >  1:1429516/1‑91 (MQ=255)
               gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                             >  1:1429424/1‑91 (MQ=255)
               gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                             >  1:671248/1‑91 (MQ=255)
               gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                             <  1:728717/91‑1 (MQ=255)
                aaGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtt                                                                            <  1:218390/91‑1 (MQ=255)
                 aGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGttt                                                                           >  1:268504/1‑91 (MQ=255)
                 aGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGttt                                                                           <  1:3229851/91‑1 (MQ=255)
                  ggTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtttt                                                                          >  1:1313631/1‑91 (MQ=255)
                     ggACAACCAATTTATCAATATTTCTTTTCAGCAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTa                                                                       >  1:3687503/1‑91 (MQ=255)
                      gACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAt                                                                      <  1:1996739/91‑1 (MQ=18)
                        cAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAttt                                                                    >  1:3140166/1‑91 (MQ=255)
                        cAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAttt                                                                    >  1:3140334/1‑91 (MQ=255)
                         aaCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAtttt                                                                   >  1:540937/1‑91 (MQ=255)
                          aCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTg                                                                  >  1:1121984/1‑91 (MQ=255)
                           ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAttttaa                                                                 <  1:631331/91‑3 (MQ=17)
                           ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                                 <  1:3866583/91‑1 (MQ=255)
                           ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                                 <  1:1022722/91‑1 (MQ=255)
                           ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                                 <  1:689652/91‑1 (MQ=255)
                            cAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCt                                                                >  1:2226611/1‑91 (MQ=255)
                            cAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCt                                                                >  1:2226571/1‑91 (MQ=255)
                               tttATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTCTGCTAtt                                                             >  1:2082718/1‑91 (MQ=255)
                                 tATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTAtttt                                                           >  1:2235962/1‑91 (MQ=255)
                                   tCAATATTTCTTTTCAGTAGAATAAAAATTTAACATACAAAACTACCAATGATAAATCACCAGTTGGTAGTTTTTTATTTTGCTATTTTAg                                                         >  1:3305490/1‑91 (MQ=255)
                                   tCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAg                                                         >  1:3049741/1‑91 (MQ=255)
                                   tCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAg                                                         >  1:373030/1‑91 (MQ=255)
                                   tCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAg                                                         >  1:616493/1‑91 (MQ=255)
                                    cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                        >  1:2018392/1‑91 (MQ=255)
                                    cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                        >  1:2018163/1‑91 (MQ=255)
                                    cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                        >  1:3982496/1‑91 (MQ=255)
                                    cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                        >  1:847651/1‑91 (MQ=255)
                                    cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                        >  1:847415/1‑91 (MQ=255)
                                    cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                        >  1:3143422/1‑91 (MQ=255)
                                    cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                        <  1:1200458/91‑1 (MQ=255)
                                     aaTATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTg                                                       >  1:1917732/1‑91 (MQ=255)
                                          ttCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATAtt                                                  <  1:1248143/91‑1 (MQ=255)
                                            cTTTTCCGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGc                                                >  1:3231847/1‑91 (MQ=255)
                                            cTTTTCAGTAGAATAAAAATTTAAAATAAAAACCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGc                                                <  1:3536855/91‑1 (MQ=255)
                                             ttttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCg                                               >  1:151469/1‑91 (MQ=255)
                                                tCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGt                                            <  1:3334599/91‑1 (MQ=255)
                                                tCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGt                                            <  1:420693/91‑1 (MQ=255)
                                                 cAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGtt                                           <  1:1467679/91‑1 (MQ=255)
                                                   gTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTaa                                         >  1:3874838/1‑91 (MQ=255)
                                                      gAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAg                                      >  1:3946994/1‑91 (MQ=255)
                                                          taaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATc                                  <  1:2779652/90‑1 (MQ=255)
                                                          aaaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATc                                  >  1:970445/1‑91 (MQ=255)
                                                          aaaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATc                                  <  1:1158159/91‑1 (MQ=255)
                                                            aaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGt                                >  1:2529614/1‑91 (MQ=255)
                                                            aaaCTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTGTCGt                                <  1:3270860/91‑1 (MQ=255)
                                                             aaTTTAAAATAACAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                               >  1:3175661/1‑91 (MQ=255)
                                                             aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                               <  1:610184/91‑1 (MQ=255)
                                                             aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                               <  1:2737954/91‑1 (MQ=255)
                                                                ttAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTtctc                            >  1:817509/1‑91 (MQ=255)
                                                                 taaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCg                           >  1:472664/1‑91 (MQ=255)
                                                                 taaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCg                           <  1:2688013/91‑1 (MQ=255)
                                                                  aaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGa                          >  1:533506/1‑91 (MQ=255)
                                                                  aaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGa                          >  1:533639/1‑91 (MQ=255)
                                                                  aaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGa                          >  1:920858/1‑91 (MQ=255)
                                                                   aaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAg                         >  1:3734291/1‑91 (MQ=255)
                                                                   aaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAg                         >  1:473042/1‑91 (MQ=255)
                                                                    aataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGt                        >  1:3472271/1‑91 (MQ=255)
                                                                     ataaaaAACTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGtt                       >  1:1228465/1‑91 (MQ=18)
                                                                     ataaaaAACTACCAACTGATGATCTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGtt                       >  1:3288385/1‑91 (MQ=18)
                                                                     ataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGtt                       >  1:207884/1‑91 (MQ=255)
                                                                     ataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGtt                       >  1:209948/1‑91 (MQ=255)
                                                                     ataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGtt                       >  1:572843/1‑91 (MQ=255)
                                                                       aaaaaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGc                     >  1:1056004/1‑91 (MQ=18)
                                                                       aaaaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGc                     <  1:1067389/91‑1 (MQ=255)
                                                                        aaaaaCTACCACTGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCt                    <  1:3764534/91‑1 (MQ=255)
                                                                        aaaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCt                    >  1:3721250/1‑91 (MQ=255)
                                                                        aaaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCt                    >  1:655976/1‑91 (MQ=255)
                                                                         aaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTa                   >  1:428530/1‑91 (MQ=255)
                                                                          aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                  >  1:494789/1‑91 (MQ=18)
                                                                          aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                  >  1:2004606/1‑91 (MQ=18)
                                                                          aaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa                  <  1:2046591/91‑1 (MQ=255)
                                                                           aaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAAc                 <  1:1121031/91‑1 (MQ=255)
                                                                            aCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACa                <  1:2448824/91‑1 (MQ=255)
                                                                             cTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACaa               >  1:639258/1‑91 (MQ=255)
                                                                              tACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAAt              <  1:3107292/91‑1 (MQ=255)
                                                                              tACCAATGATAAAGCATCAGTTGGTAGTTTTTTATTTTGCTATGTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAAt              <  1:4150861/91‑1 (MQ=255)
                                                                                ccAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGt            <  1:2494706/91‑1 (MQ=255)
                                                                                 cAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTc           <  1:2786600/91‑1 (MQ=255)
                                                                                      gATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGt       >  1:417278/1‑91 (MQ=255)
                                                                                      gATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGt       <  1:2167517/91‑1 (MQ=255)
                                                                                          aaTCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCaa   >  1:1679714/1‑91 (MQ=255)
                                                                                           atcatcAGTTGGTAGTTTTTTATTTTGCTATTATAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCAAc  <  1:3545854/91‑1 (MQ=255)
                                                                                           |                                                                                          
GAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTCAAC  >  USA300TCH1516_ALE/1239474‑1239654

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