Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 536,863 T→A 8.5% intergenic (+56/‑255) 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,8630TA 8.5% 754.2 / 44.8 228intergenic (+56/‑255)rplY/pth50S ribosomal protein L25/Peptidyl‑tRNA hydrolase
Reads supporting (aligned to +/- strand):  ref base T (111/97);  new base A (4/15);  total (116/112)
Fisher's exact test for biased strand distribution p-value = 7.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

Only 100 of 235 total aligned reads displayed.

AAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAA  >  USA300TCH1516_ALE/536773‑536953
                                                                                          |                                                                                          
 aaGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTAttt                                                                                           >  2:3473567/1‑91 (MQ=255)
   gCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTa                                                                                         >  1:1231799/1‑91 (MQ=255)
    cAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATATTaa                                                                                        <  1:2368314/91‑1 (MQ=255)
     aaaGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAat                                                                                       >  1:1789297/1‑91 (MQ=255)
     aaaGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATATTaaa                                                                                       <  2:3113625/91‑2 (MQ=255)
     aaaGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGGTTAACAAGCACTGTGCTTATATTaaa                                                                                       <  2:2248790/91‑2 (MQ=255)
      aagaagACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAata                                                                                      >  2:1538189/1‑91 (MQ=255)
      aagaagACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAata                                                                                      >  1:2519407/1‑91 (MQ=255)
      aagaagACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAata                                                                                      >  1:168061/1‑91 (MQ=255)
       agaagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAatat                                                                                     <  1:3460875/91‑1 (MQ=255)
       agaagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAatat                                                                                     >  2:3389963/1‑91 (MQ=255)
       agaagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAatat                                                                                     >  1:2991971/1‑91 (MQ=255)
       agaagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAatat                                                                                     >  1:3321874/1‑91 (MQ=255)
        gaagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAATGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAAtata                                                                                    >  2:3542715/1‑91 (MQ=255)
         aagaCAAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATaa                                                                                   <  2:3050040/91‑1 (MQ=255)
           gaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGc                                                                                 >  2:1575713/1‑91 (MQ=255)
           gaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGc                                                                                 >  1:1445659/1‑91 (MQ=255)
           gaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGc                                                                                 >  1:777561/1‑91 (MQ=255)
            aCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCa                                                                                <  1:3510055/91‑1 (MQ=38)
             cGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCac                                                                               <  1:3457893/91‑2 (MQ=38)
              gaagaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATg                                                                              <  2:1920797/91‑1 (MQ=255)
              gaagaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATg                                                                              >  1:2491163/1‑91 (MQ=255)
               aagaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATgg                                                                             >  1:3255875/1‑91 (MQ=255)
               aagaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATgg                                                                             >  1:2091997/1‑91 (MQ=255)
               aagaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATgg                                                                             <  2:429161/91‑1 (MQ=255)
               aagaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATgg                                                                             <  1:359855/91‑1 (MQ=255)
                agaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCacagt                                                                            <  1:1612152/91‑5 (MQ=25)
                 gaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTg                                                                           <  1:2928751/91‑1 (MQ=255)
                 gaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTg                                                                           >  1:1748014/1‑91 (MQ=255)
                   aaaaaCTGAAGAGTAATTTTAATCTGTTACATTAACGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCt                                                                         <  2:332059/91‑1 (MQ=255)
                   aaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCt                                                                         >  1:2112376/1‑91 (MQ=255)
                    aaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtt                                                                        >  2:754634/1‑91 (MQ=255)
                     acaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTGTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                       <  2:3122077/89‑1 (MQ=255)
                      aaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtttt                                                                      >  1:556218/1‑91 (MQ=255)
                       aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                     >  2:3078190/1‑91 (MQ=255)
                       aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                     >  1:673537/1‑91 (MQ=255)
                        cTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtg                                                                    <  2:253219/91‑1 (MQ=255)
                        cTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtg                                                                    <  1:1399425/91‑1 (MQ=255)
                          gAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgtg                                                                  <  1:1674560/91‑1 (MQ=255)
                          gAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgtg                                                                  <  2:424117/91‑1 (MQ=255)
                           atgagTAATTTTAATCTGTTACATTACAGTTTTGATACTTTGTTTGACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTgtgt                                                                 <  2:955316/89‑1 (MQ=255)
                           aaGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTgtgt                                                                 <  1:2568176/91‑1 (MQ=255)
                           aaGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTgtgt                                                                 <  1:644717/91‑1 (MQ=255)
                            agagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtt                                                                <  1:3740823/91‑1 (MQ=255)
                            agagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtt                                                                <  2:535496/91‑1 (MQ=255)
                            agagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtt                                                                <  1:533516/91‑1 (MQ=255)
                             gagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtta                                                               >  1:3691832/1‑91 (MQ=255)
                              agTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttat                                                              <  2:1127044/91‑1 (MQ=255)
                              agTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttat                                                              <  1:3206731/91‑1 (MQ=255)
                                tAATTTTAATCTGTTACATTAAAGATTTTAGACTTTGATTAACAAGCACTGTGCTTCTTTTAATATACGCCTGCTGCTCTTTGTGttattc                                                            >  2:920672/1‑90 (MQ=255)
                                 aaTTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTtattat                                                           <  2:169617/91‑1 (MQ=255)
                                  aTTTTAGTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATtata                                                          <  1:821660/91‑1 (MQ=255)
                                   ttttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaa                                                         >  2:1839510/1‑91 (MQ=255)
                                    tttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaaa                                                        >  1:895191/1‑91 (MQ=255)
                                       aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCt                                                     >  1:1618303/1‑91 (MQ=255)
                                           tGTTACATTAAAGTTTTTATACTTTGTTTAACAGGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaat                                                 <  1:1511792/91‑1 (MQ=255)
                                           tGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaat                                                 <  2:3297062/91‑1 (MQ=255)
                                            gTTACATTAAAGTTTGTATACTTTGGGTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaatt                                                <  2:3425673/91‑1 (MQ=255)
                                              tACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaattaa                                              <  1:1025988/91‑1 (MQ=255)
                                                cATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAAc                                            <  1:1086088/91‑1 (MQ=255)
                                                 aTTAAAGTTTTTATACTTTGTTTAACAAGCTCTGTGCTTATTTTAATATAAGCATGTTGCTTTTTGTGGTATTATAAAGCTTAATTAAACt                                           >  2:265418/1‑91 (MQ=255)
                                                  ttAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACtt                                          >  1:3128895/1‑91 (MQ=255)
                                                     aaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttat                                       >  1:418682/1‑91 (MQ=255)
                                                      aGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttatt                                      <  1:1657552/91‑1 (MQ=255)
                                                      aGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttatt                                      >  1:697412/1‑91 (MQ=255)
                                                      aGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTGTTTGTGCTATTATAAAGCTTAATTAAACTttatt                                      >  2:1279842/1‑91 (MQ=255)
                                                       gTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttatta                                     <  2:279000/91‑1 (MQ=255)
                                                        tttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTAc                                    >  1:318997/1‑91 (MQ=255)
                                                          tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                  >  2:1074852/1‑91 (MQ=255)
                                                          tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                  >  2:3370453/1‑91 (MQ=255)
                                                          tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                  >  1:2865093/1‑91 (MQ=255)
                                                          tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                  >  1:2865735/1‑91 (MQ=255)
                                                           ttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACttt                                 >  2:140348/1‑91 (MQ=255)
                                                           ttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACttt                                 >  1:3733948/1‑91 (MQ=255)
                                                           ttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACttt                                 >  1:3934650/1‑91 (MQ=255)
                                                           ttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACttt                                 >  2:2390625/1‑91 (MQ=255)
                                                              taCTTTGTTTAACAAGCACTGTGCCTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTa                              <  2:3226447/91‑1 (MQ=255)
                                                                cTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACt                            >  1:515305/1‑91 (MQ=255)
                                                                 tttGTTTAACAAGCACTGTGCTGATTTTAATATCAGCATGGTGCTATTTGTGTGATTATATAGCTTCATTAGACTTTATTACGTTgtagta                           >  2:649690/1‑91 (MQ=255)
                                                                  ttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTaa                          >  2:2169474/1‑91 (MQ=255)
                                                                  ttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTaa                          >  1:1479471/1‑91 (MQ=255)
                                                                  ttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTaa                          >  1:3351772/1‑91 (MQ=255)
                                                                   tGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTaaa                         >  1:1763028/1‑91 (MQ=255)
                                                                       tAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGttt                     >  2:2153917/1‑91 (MQ=255)
                                                                        aaCAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTtta                    >  1:2371218/1‑91 (MQ=255)
                                                                          cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                  >  2:3691881/1‑91 (MQ=255)
                                                                          cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                  <  2:1045171/91‑1 (MQ=255)
                                                                           aaGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaatt                 <  2:3255875/91‑1 (MQ=255)
                                                                           aaGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaatt                 <  1:3304182/91‑1 (MQ=255)
                                                                           aaGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaatt                 >  2:1633572/1‑91 (MQ=255)
                                                                               accatgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTaattaatt             >  2:1970350/5‑91 (MQ=25)
                                                                                cTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTAGTACAAAGTAATAAAGTttaattaagct            <  1:957965/91‑4 (MQ=21)
                                                                                 catgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAtttt           >  1:2812917/3‑91 (MQ=37)
                                                                                   tgCTTATTTTAATATAAGCATGGTGCTTGTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAg         >  2:2765543/1‑91 (MQ=255)
                                                                                   tgCTTATTTCAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAg         >  1:1971645/1‑91 (MQ=255)
                                                                                    gCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGt        >  2:823211/1‑91 (MQ=255)
                                                                                    gCTTATATTAAAATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGt        >  1:3204641/1‑91 (MQ=255)
                                                                                      ttATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGa      >  1:594068/1‑91 (MQ=255)
                                                                                      ttATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGa      >  2:2824231/1‑91 (MQ=255)
                                                                                          tttAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaa  <  2:2026199/91‑1 (MQ=255)
                                                                                          |                                                                                          
AAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAA  >  USA300TCH1516_ALE/536773‑536953

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