Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 536,868 T→A 9.9% intergenic (+61/‑250) rplY → / → pth 50S ribosomal protein L25/Peptidyl‑tRNA hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE536,8680TA 9.9% 708.7 / 54.2 223intergenic (+61/‑250)rplY/pth50S ribosomal protein L25/Peptidyl‑tRNA hydrolase
Reads supporting (aligned to +/- strand):  ref base T (99/101);  new base A (11/11);  total (111/112)
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.48e-01

Only 100 of 244 total aligned reads displayed.

AAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGAC  >  USA300TCH1516_ALE/536778‑536958
                                                                                          |                                                                                          
aaaGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTATATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAat                                                                                            >  2:3058304/1‑91 (MQ=255)
aaaGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAat                                                                                            <  1:526621/91‑1 (MQ=255)
  agaagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAatat                                                                                          >  1:3321874/1‑91 (MQ=255)
  agaagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAatat                                                                                          <  1:3460875/91‑1 (MQ=255)
  agaagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAatat                                                                                          >  1:2991971/1‑91 (MQ=255)
   gaagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAATGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAAtata                                                                                         >  2:3542715/1‑91 (MQ=255)
      gaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGc                                                                                      >  1:777561/1‑91 (MQ=255)
      gaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGc                                                                                      >  2:2799970/1‑91 (MQ=255)
      gaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGc                                                                                      <  2:2262650/91‑1 (MQ=38)
       aCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCa                                                                                     >  1:2094922/1‑91 (MQ=255)
        cGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCac                                                                                    <  1:3457893/91‑2 (MQ=38)
         gaagaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATg                                                                                   >  1:334463/1‑91 (MQ=255)
          aagaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATgg                                                                                  >  1:918312/1‑91 (MQ=255)
           agaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCacagt                                                                                 <  1:1612152/91‑5 (MQ=25)
           agaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCacagt                                                                                 <  1:3764421/91‑5 (MQ=25)
            gaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTg                                                                                >  1:1748014/1‑91 (MQ=255)
            gaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTg                                                                                <  1:2928751/91‑1 (MQ=255)
              aaaaaCTGAAGAGTAATTTTAATCTGTTACATTAACGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCt                                                                              <  2:332059/91‑1 (MQ=255)
               aaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtt                                                                             >  2:754634/1‑91 (MQ=255)
               aaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtt                                                                             <  1:3198328/91‑1 (MQ=255)
               aaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtt                                                                             <  1:465554/91‑1 (MQ=255)
                aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTCTACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                            >  2:3254132/1‑91 (MQ=255)
                 aaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAGGCACTGTGCTTATTTTAATATAAGCATGGTGCtttt                                                                           <  1:1460402/91‑1 (MQ=255)
                 aaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtttt                                                                           >  1:556218/1‑91 (MQ=255)
                 aaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtttt                                                                           <  1:1365107/91‑1 (MQ=255)
                  aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                          >  1:673563/1‑91 (MQ=255)
                  aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                          >  1:169330/1‑91 (MQ=255)
                   cTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtg                                                                         <  1:1399425/91‑1 (MQ=255)
                    tGAAGAGTTATTTTAATCTGTTACATGATAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgt                                                                        <  2:1095152/91‑1 (MQ=255)
                    tGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgt                                                                        >  1:2591248/1‑91 (MQ=255)
                     gAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgtg                                                                       <  1:1674560/91‑1 (MQ=255)
                     gAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgtg                                                                       <  2:424144/91‑1 (MQ=255)
                      atgagTAATTTTAATCTGTTACATTACAGTTTTGATACTTTGTTTGACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTgtgt                                                                      <  2:955316/89‑1 (MQ=255)
                      aaGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTgtgt                                                                      >  1:824001/1‑91 (MQ=255)
                      aaGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTgtgt                                                                      <  1:644717/91‑1 (MQ=255)
                      aaGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTgtgt                                                                      <  1:2568176/91‑1 (MQ=255)
                       agagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtt                                                                     <  1:3740823/91‑1 (MQ=255)
                       agagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtt                                                                     <  2:535496/91‑1 (MQ=255)
                         agTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttat                                                                   <  1:3206731/91‑1 (MQ=255)
                           tAATTTTAATCTGTTACATTAAAGATTTTAGACTTTGATTAACAAGCACTGTGCTTCTTTTAATATACGCCTGCTGCTCTTTGTGttattc                                                                 >  2:920672/1‑90 (MQ=255)
                                 tAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGc                                                           <  2:216278/91‑1 (MQ=255)
                                 tAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGc                                                           <  2:1484529/91‑1 (MQ=255)
                                 tAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAACATAAGCATGGTGCTTTTTGTGTTATTATAAAGc                                                           >  1:3239450/1‑91 (MQ=255)
                                  aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCt                                                          >  1:1618303/1‑91 (MQ=255)
                                    tCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCtta                                                        <  2:1062657/91‑1 (MQ=255)
                                      tGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaat                                                      <  1:3180548/91‑1 (MQ=255)
                                      tGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaat                                                      >  2:1551880/1‑91 (MQ=255)
                                      tGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTACTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaat                                                      >  1:3665665/1‑91 (MQ=255)
                                         tACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaattaa                                                   <  1:3595799/91‑1 (MQ=255)
                                         tACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaattaa                                                   <  2:381986/91‑1 (MQ=255)
                                         tACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaattaa                                                   <  1:1025988/91‑1 (MQ=255)
                                           cATTAAAGTTTTTATTCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAAc                                                 >  2:2777931/1‑91 (MQ=255)
                                           cATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAAc                                                 <  1:1086088/91‑1 (MQ=255)
                                            aTTAAAGTTTTTATACTTTGTTTAACAAGCTCTGTGCTTATTTTAATATAAGCATGTTGCTTTTTGTGGTATTATAAAGCTTAATTAAACt                                                >  2:265418/1‑91 (MQ=255)
                                               aaaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTtta                                             <  1:190075/91‑1 (MQ=255)
                                               aaaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTtta                                             <  2:2752089/91‑1 (MQ=255)
                                                 aGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttatt                                           >  1:697615/1‑91 (MQ=255)
                                                 aGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTGTTTGTGCTATTATAAAGCTTAATTAAACTttatt                                           >  2:1279842/1‑91 (MQ=255)
                                                     tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                       >  2:3370453/1‑91 (MQ=255)
                                                     tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                       >  1:2865093/1‑91 (MQ=255)
                                                     tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                       <  1:3936096/91‑1 (MQ=255)
                                                     tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                       >  1:2865735/1‑91 (MQ=255)
                                                      ttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACttt                                      >  2:1693281/1‑91 (MQ=255)
                                                       tataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTg                                     <  2:3265986/91‑1 (MQ=255)
                                                       tataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTg                                     >  2:1864271/1‑91 (MQ=255)
                                                        ataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGt                                    >  1:2496317/1‑91 (MQ=255)
                                                         taCTTTGTTTAACAAGCACTGTGCTTATTTTCATATAAGCACGGTGCTTTTTGTGTTATTATAAAGCACAATTAAACTTCATTACTTTGTa                                   >  2:3013548/1‑91 (MQ=255)
                                                         taCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTa                                   >  2:1625266/1‑91 (MQ=255)
                                                         taCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTa                                   >  2:621950/1‑91 (MQ=255)
                                                           cTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCGTAATTAAACGTTATTACTTTGTACt                                 >  2:3422040/1‑91 (MQ=255)
                                                             ttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTaa                               <  2:3108477/91‑1 (MQ=255)
                                                              tGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTaaa                              >  1:1763028/1‑91 (MQ=255)
                                                               gTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAg                             <  1:150636/91‑1 (MQ=255)
                                                                   aaCAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTtta                         >  1:2371183/1‑91 (MQ=255)
                                                                   aaCAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTAGTACAAAGTAATAAAGTtta                         <  2:3428436/91‑1 (MQ=25)
                                                                     cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                       <  2:2525636/91‑1 (MQ=255)
                                                                     cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                       <  2:1043095/91‑1 (MQ=255)
                                                                      aaGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaatt                      <  1:3304182/91‑1 (MQ=255)
                                                                       aGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaatta                     >  2:3128514/1‑91 (MQ=255)
                                                                        gCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaattaa                    <  1:1627624/91‑1 (MQ=255)
                                                                         cACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTtaattaat                   >  2:3698717/1‑91 (MQ=255)
                                                                          accatgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTaattaatt                  >  2:1969655/5‑91 (MQ=25)
                                                                            tgtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAtttt                <  1:3699479/91‑1 (MQ=255)
                                                                             gtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTa               >  1:2410377/1‑91 (MQ=255)
                                                                              tgCTTATTTCAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAg              >  1:1971645/1‑91 (MQ=255)
                                                                               gCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGt             >  2:823211/1‑91 (MQ=255)
                                                                                 ttATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGa           >  2:2824231/1‑91 (MQ=255)
                                                                                 ttATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGa           >  1:594068/1‑91 (MQ=255)
                                                                                  tatattAAAATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAg          >  1:3053827/5‑91 (MQ=255)
                                                                                    ttttAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTa        >  2:2006109/1‑91 (MQ=255)
                                                                                    tattAAAATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTa        >  2:2211662/3‑91 (MQ=255)
                                                                                    tattAAAATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTa        >  1:669728/3‑91 (MQ=255)
                                                                                    tattAAAATAAGCATGGTGCTTTTTGTGTTATTATAAAGATTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTa        >  2:358969/3‑91 (MQ=255)
                                                                                     attAAAATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaa       >  2:1750534/2‑91 (MQ=255)
                                                                                       tAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaaa     <  1:3654024/91‑1 (MQ=255)
                                                                                       tAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaaa     <  2:1148789/91‑1 (MQ=255)
                                                                                        aaTATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAg    <  1:1010114/91‑1 (MQ=255)
                                                                                          tataAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGAc  <  1:2182647/91‑1 (MQ=255)
                                                                                          tataAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGAc  >  2:2084593/1‑91 (MQ=255)
                                                                                          tataAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGAc  <  2:3227084/91‑1 (MQ=255)
                                                                                          |                                                                                          
AAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGAC  >  USA300TCH1516_ALE/536778‑536958

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