Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,366,121 A→T 12.2% intergenic (+420/‑34) rpmG2_1 → / → rpsN2 50S ribosomal protein L33 2/Alternate 30S ribosomal protein S14

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,366,1210AT12.2% 654.6 / 70.0 221intergenic (+420/‑34)rpmG2_1/rpsN250S ribosomal protein L33 2/Alternate 30S ribosomal protein S14
Reads supporting (aligned to +/- strand):  ref base A (96/98);  new base T (14/13);  total (110/111)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 249 total aligned reads displayed.

TAAATATGATGGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGTAAATAAA  >  USA300TCH1516_ALE/1366031‑1366211
                                                                                          |                                                                                          
tAAATATGATGGTTTAGAAAATGCTTTAACAAAGAAATGAAACGTTACTGGTGAGTTATGTGAGGATTGTGTTATTATATAAAGCGtaata                                                                                            <  2:3607595/91‑1 (MQ=255)
 aaaTATGATGGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGtaataa                                                                                           >  2:2648252/1‑91 (MQ=255)
    tatGATGGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTa                                                                                        <  1:716132/91‑1 (MQ=255)
    tatGATGGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTa                                                                                        <  1:3231804/91‑1 (MQ=255)
     atgatgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATTATTAc                                                                                       <  2:870889/91‑1 (MQ=255)
     atgatgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATTATTAc                                                                                       <  1:1525314/91‑1 (MQ=255)
       gatgGTTTAGAACATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGa                                                                                     >  1:2975494/1‑91 (MQ=255)
       gatgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATTATTACGa                                                                                     >  1:1335145/1‑91 (MQ=255)
       gatgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGa                                                                                     >  1:2574997/1‑91 (MQ=255)
        atgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGAt                                                                                    >  2:1773909/1‑91 (MQ=255)
        atgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGAt                                                                                    >  1:1467003/1‑91 (MQ=255)
         tgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGAtt                                                                                   <  2:1756787/91‑1 (MQ=255)
         tgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGAtt                                                                                   >  2:2286907/1‑91 (MQ=255)
           gTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATTATTACGAttta                                                                                 <  1:59675/91‑2 (MQ=255)
           gTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTg                                                                                 >  1:1432856/1‑91 (MQ=255)
            tttAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGa                                                                                <  1:2102564/91‑1 (MQ=255)
            tttAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGa                                                                                <  1:1078296/91‑1 (MQ=255)
             ttAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATTATTACGAtttata                                                                               <  2:2998505/91‑4 (MQ=255)
             ttAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATTATTACGAtttata                                                                               <  1:3385988/91‑4 (MQ=255)
             ttAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGAt                                                                               <  1:3815458/91‑1 (MQ=255)
              tAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATa                                                                              >  1:355069/1‑91 (MQ=255)
              tAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATa                                                                              >  1:1111089/1‑91 (MQ=255)
               aGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATaa                                                                             >  1:2636086/1‑91 (MQ=255)
                gAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATaaa                                                                            <  2:574590/91‑1 (MQ=255)
                 aaaaTGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATaaaa                                                                           <  2:1179580/91‑1 (MQ=255)
                      gATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGa                                                                      <  1:2228071/91‑1 (MQ=255)
                      gATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGa                                                                      <  2:2429895/91‑1 (MQ=255)
                      gATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAAGCGTAATAATTACGATTTGATAAAAAGTGa                                                                      <  1:1644586/91‑1 (MQ=255)
                       aTTTAACAAAGAAATGAAACTTTACTCTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAg                                                                     >  2:2208346/1‑91 (MQ=255)
                        tttAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAgg                                                                    <  2:130637/91‑1 (MQ=255)
                         ttAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGt                                                                   >  2:2890962/1‑91 (MQ=255)
                         ttAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGt                                                                   >  2:2035604/1‑91 (MQ=255)
                         ttAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGt                                                                   >  2:2167178/1‑91 (MQ=255)
                          tAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTa                                                                  >  2:1643629/1‑91 (MQ=255)
                          tAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTa                                                                  <  2:3212378/91‑1 (MQ=255)
                            aCAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAAc                                                                <  1:2956953/91‑1 (MQ=255)
                            aCAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAAc                                                                >  1:77389/1‑91 (MQ=255)
                            aCAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAAc                                                                <  2:1866299/91‑1 (MQ=255)
                                aGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACtata                                                            >  2:2598321/1‑91 (MQ=255)
                                aGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACtata                                                            <  2:3252581/91‑1 (MQ=255)
                                 gAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACtatat                                                           >  1:3493065/1‑91 (MQ=255)
                                 gAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACtatat                                                           <  1:3637425/91‑1 (MQ=255)
                                   aaTGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTatatat                                                         <  2:2330624/91‑1 (MQ=255)
                                    aTGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATg                                                        <  2:206074/91‑1 (MQ=255)
                                     tGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATgg                                                       >  2:706562/1‑91 (MQ=255)
                                     tGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATgg                                                       <  1:1799904/91‑1 (MQ=255)
                                     tGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATgg                                                       >  1:3167371/1‑91 (MQ=255)
                                       aaaCTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCt                                                     >  1:2049031/1‑91 (MQ=255)
                                       aaaCTTTACTGTTGAATTATGTGAGGATAGTGTTATTATATAAATCGTAATAATTACGCTTTGATAAAAAGTGAGGTAACTATATATGGCt                                                     >  2:492111/1‑91 (MQ=255)
                                        aaCTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATTGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTa                                                    <  2:3480326/91‑1 (MQ=255)
                                        aaCTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTa                                                    >  2:2194043/1‑91 (MQ=255)
                                        aaCTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTa                                                    <  2:2379161/91‑1 (MQ=255)
                                         aCTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTaa                                                   >  2:3862605/1‑91 (MQ=255)
                                          cTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAg                                                  >  1:3862246/1‑91 (MQ=255)
                                          cTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAg                                                  <  2:344256/91‑1 (MQ=255)
                                          cTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAg                                                  <  2:1467003/91‑1 (MQ=255)
                                           tttACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGa                                                 <  1:463774/91‑1 (MQ=255)
                                             tACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGaaa                                               >  1:2557792/1‑91 (MQ=255)
                                             tACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGaaa                                               >  2:272795/1‑91 (MQ=255)
                                             tACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGaaa                                               <  1:1603208/91‑1 (MQ=255)
                                                tGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCt                                            <  2:1819510/91‑1 (MQ=255)
                                                      aTTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAATGAGGTAACTATATATGGCTAAGAAATCTAAAATa                                      <  2:3650734/91‑1 (MQ=255)
                                                       ttATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAg                                     >  2:890865/1‑91 (MQ=255)
                                                       ttATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAg                                     <  2:1975504/91‑1 (MQ=255)
                                                        tATGTGAGGATTGTGTTATTCTATAGATCGTAATATTTACGATTTGATAAAAAGTGAGGTAACTATATATGTCTAAGAAATCTAAAATAGc                                    <  2:2108713/91‑1 (MQ=255)
                                                           gtgAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCaaa                                 >  2:3258459/1‑91 (MQ=255)
                                                             gAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGGGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAag                               <  1:1773892/91‑1 (MQ=255)
                                                             gAGGATTCTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAag                               <  1:1941436/91‑1 (MQ=255)
                                                              aGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGAaaaaaa                              <  2:3814291/91‑3 (MQ=255)
                                                                      gtTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGaaaaa                      <  1:3352773/91‑1 (MQ=255)
                                                                      gtTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGaaaaa                      <  1:119847/91‑1 (MQ=255)
                                                                      gtTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGaaaaa                      >  1:3257126/1‑91 (MQ=255)
                                                                           tatataAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAaagagaa                 <  1:1095423/91‑1 (MQ=255)
                                                                           tatataAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAaagagaa                 <  1:2320227/91‑1 (MQ=255)
                                                                               taAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGt             <  1:1931757/91‑1 (MQ=255)
                                                                               caaaTCGTAATTATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGt             >  1:656157/2‑91 (MQ=255)
                                                                                aaaTCGTAATTATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGtt            >  2:2340072/1‑91 (MQ=255)
                                                                                aaaTCGTAATTATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGtt            >  2:3782612/1‑91 (MQ=255)
                                                                                aaaTCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGtt            >  1:1422593/1‑91 (MQ=255)
                                                                                aaaTCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGtt            <  1:1713343/91‑1 (MQ=255)
                                                                                aaaTCGTAATAATTACGATATGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGtt            <  2:2254720/91‑1 (MQ=255)
                                                                                 aaTCGTAATTATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTa           >  1:3041631/1‑91 (MQ=255)
                                                                                  aTCGTAATTATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAg          >  1:1388537/1‑91 (MQ=255)
                                                                                  aTCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAg          >  1:461891/1‑91 (MQ=255)
                                                                                  aTCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAg          >  1:416062/1‑91 (MQ=255)
                                                                                  aTCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAg          >  1:808436/1‑91 (MQ=255)
                                                                                   tCGTAATTATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGt         >  1:2723429/1‑91 (MQ=255)
                                                                                   tCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGt         >  2:2572770/1‑91 (MQ=255)
                                                                                     gTAATTATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAGAAGAGAAGAGTTAGtaa       >  1:3329593/1‑91 (MQ=255)
                                                                                     gTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaa       <  2:1001918/91‑1 (MQ=255)
                                                                                     gTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaa       >  1:216705/1‑91 (MQ=255)
                                                                                      taattaTTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaaa      >  2:197000/1‑91 (MQ=255)
                                                                                      taattaTTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaaa      >  1:527787/1‑91 (MQ=255)
                                                                                      taattaTTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaaa      >  2:1290815/1‑91 (MQ=255)
                                                                                      taattaTTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaaa      >  2:3681171/1‑91 (MQ=255)
                                                                                      taattaTTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaaa      >  1:2839874/1‑91 (MQ=255)
                                                                                      taattaTTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaaa      >  1:416709/1‑91 (MQ=255)
                                                                                      taataaTTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaaa      >  1:791199/1‑91 (MQ=255)
                                                                                       aataatTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaaat     >  2:281223/1‑91 (MQ=255)
                                                                                       aataatTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaaat     <  2:3318219/91‑1 (MQ=255)
                                                                                          |                                                                                          
TAAATATGATGGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGTAAATAAA  >  USA300TCH1516_ALE/1366031‑1366211

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