Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,779,807 A→G 6.2% intergenic (‑113/+29) USA300TCH1516_01664 ← / ← rplT hypothetical protein/50S ribosomal protein L20

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,779,8070AG 6.2% 632.3 / 20.7 209intergenic (‑113/+29)USA300TCH1516_01664/rplThypothetical protein/50S ribosomal protein L20
Reads supporting (aligned to +/- strand):  ref base A (92/103);  new base G (3/10);  total (96/113)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 218 total aligned reads displayed.

TGTTTGTTCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTG  >  USA300TCH1516_ALE/1779717‑1779897
                                                                                          |                                                                                          
 gtttgttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGAc                                                                                           <  2:2508728/91‑1 (MQ=255)
  tttgttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACa                                                                                          <  2:1909167/91‑1 (MQ=255)
  tttgttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACa                                                                                          <  1:1156589/91‑1 (MQ=255)
  tttgttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACa                                                                                          >  2:1097978/1‑91 (MQ=255)
  tttgttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACa                                                                                          >  1:2220769/1‑91 (MQ=255)
   ttgttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACaa                                                                                         <  2:1969058/91‑1 (MQ=255)
   ttgttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACaa                                                                                         <  2:47872/91‑1 (MQ=255)
    tgttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAAt                                                                                        >  1:1143575/1‑91 (MQ=255)
    tgttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAAt                                                                                        <  1:1728091/91‑1 (MQ=255)
     gttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATg                                                                                       <  1:825824/91‑1 (MQ=255)
       tctcTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGcc                                                                                     <  2:1585506/91‑1 (MQ=255)
        ctctTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCa                                                                                    <  2:201088/91‑1 (MQ=255)
         tctTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCAt                                                                                   <  2:2856354/91‑1 (MQ=255)
         tctTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCAt                                                                                   <  1:776572/91‑1 (MQ=255)
          ctTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCAtc                                                                                  >  1:2412512/1‑91 (MQ=255)
          ctTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCAtc                                                                                  <  2:853547/91‑1 (MQ=255)
           ttAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCAtct                                                                                 >  1:2757020/1‑91 (MQ=255)
             aaaTGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATctct                                                                               <  1:3866357/91‑1 (MQ=255)
              aaTGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTTTTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCtt                                                                              <  1:3279093/91‑1 (MQ=38)
              aaTGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCtt                                                                              <  2:2413613/91‑1 (MQ=38)
               aTGCATTTTAGCAAATATGTCATATTTTCGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCttt                                                                             >  2:397853/1‑91 (MQ=255)
               aTGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCttt                                                                             <  1:751972/91‑1 (MQ=38)
                tGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCtttt                                                                            >  2:358556/1‑91 (MQ=255)
                  cATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCtttttt                                                                          <  2:2356791/91‑1 (MQ=38)
                   aTTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCTTTTTta                                                                         <  2:2937330/91‑1 (MQ=37)
                    ttttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCTTTTTtat                                                                        <  1:2603808/91‑1 (MQ=37)
                    ttttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTtat                                                                        <  1:3433377/91‑1 (MQ=255)
                    ttttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTtat                                                                        <  1:2365936/91‑1 (MQ=255)
                     tttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTtata                                                                       >  2:1175461/1‑91 (MQ=255)
                      ttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATaa                                                                      >  2:1205011/1‑91 (MQ=255)
                      ttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTCGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATaa                                                                      >  2:2948738/1‑91 (MQ=255)
                       tAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAg                                                                     >  2:1540237/1‑91 (MQ=255)
                         gCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAgtg                                                                   <  2:2395073/91‑1 (MQ=255)
                          cAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAgtgt                                                                  <  1:2578690/91‑1 (MQ=255)
                          cAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAgtgt                                                                  >  2:2220826/1‑91 (MQ=255)
                             atatGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGttat                                                               <  1:683916/91‑1 (MQ=255)
                             atatGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTATTTATAAGTGttat                                                               >  2:534857/1‑91 (MQ=255)
                             atatGTCATAGTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAGAAGGGATGACCATGCCATCTCTTTTTTATAAGTGttat                                                               <  1:1982803/91‑1 (MQ=255)
                               atGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGttatta                                                             <  2:1971166/91‑1 (MQ=255)
                               atGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGttatta                                                             <  1:2873960/91‑1 (MQ=255)
                               atGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGttatta                                                             >  2:1087569/1‑91 (MQ=255)
                                 gTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTattatt                                                           <  1:2808551/91‑1 (MQ=255)
                                 gTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTattatt                                                           <  2:2086426/91‑1 (MQ=255)
                                   cATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTAtttt                                                         >  2:3315582/1‑91 (MQ=255)
                                   cATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTAtttt                                                         >  2:781747/1‑91 (MQ=255)
                                    atatTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTa                                                        >  1:1242268/1‑91 (MQ=255)
                                    atatTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTa                                                        >  1:1001267/1‑91 (MQ=255)
                                       ttttAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAg                                                     <  1:2520417/91‑1 (MQ=255)
                                       ttttAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAg                                                     <  2:3014212/91‑1 (MQ=255)
                                        tttAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGc                                                    >  2:3269208/1‑91 (MQ=255)
                                         ttAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCa                                                   >  2:344602/1‑91 (MQ=255)
                                         ttAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCa                                                   <  2:2673098/91‑1 (MQ=255)
                                          tAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCAt                                                  >  1:1775881/1‑91 (MQ=255)
                                           aGCGATTAAATGTGGCATTTTGCTATCTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATc                                                 >  2:2205848/1‑91 (MQ=255)
                                             cGATTATATGTGGCATTTAGCTATTTAGTAAATACAAAAGCGATGACAATGCCATCTCTTTTTTATCAGTGTTATTATTTTAAAGCATctt                                               <  2:2778432/91‑1 (MQ=255)
                                             cGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctt                                               >  1:1705241/1‑91 (MQ=255)
                                             cGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctt                                               <  1:1910208/91‑1 (MQ=255)
                                               aTTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTAAAGCATctttag                                            <  2:911722/91‑1 (MQ=255)
                                                    aTGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttagcttt                                        >  2:829177/1‑91 (MQ=255)
                                                    aTGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttagcttt                                        >  2:1296125/1‑91 (MQ=255)
                                                      gtgGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttagctttag                                      <  2:3694308/91‑1 (MQ=255)
                                                      gtgGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttagctttag                                      >  2:3052447/1‑91 (MQ=255)
                                                        ggCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGtt                                    <  2:3319670/91‑1 (MQ=255)
                                                        ggCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGtt                                    <  1:489403/91‑1 (MQ=255)
                                                         gCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTa                                   <  1:1175461/91‑1 (MQ=255)
                                                          cATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTAc                                  <  1:2931725/91‑1 (MQ=255)
                                                            ttttGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTa                                >  1:974994/1‑91 (MQ=255)
                                                             tttGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTaa                               <  1:1515457/91‑1 (MQ=255)
                                                              ttGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAAt                              >  1:1762418/1‑91 (MQ=255)
                                                                   aTTTAAAAAATATAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTCAGc                         >  1:3002133/1‑91 (MQ=255)
                                                                    tttAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCa                        >  1:3692172/1‑91 (MQ=255)
                                                                    tttAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCa                        <  1:3592861/91‑1 (MQ=255)
                                                                      taaaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCaaa                      <  2:2056034/91‑1 (MQ=255)
                                                                        aaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAATCATCTTTAGCTTTAGTTACTAATTGAGCAAATg                    >  2:1368718/1‑91 (MQ=255)
                                                                        aaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATg                    <  2:2244308/91‑1 (MQ=255)
                                                                         aaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGc                   <  2:1421196/91‑1 (MQ=255)
                                                                         aaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGc                   >  1:795805/1‑91 (MQ=255)
                                                                         aaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGc                   >  1:652105/1‑91 (MQ=255)
                                                                          aaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTCTTTTACAGCACCTTTAGCCTTAGTTACTCATTGAGCAAATCCt                  >  2:3460400/1‑91 (MQ=255)
                                                                           aataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCtt                 <  1:977156/91‑1 (MQ=255)
                                                                            aataaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCttt                >  2:1444034/4‑91 (MQ=255)
                                                                             taaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCtttt               >  1:3684082/1‑91 (MQ=255)
                                                                             taaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCtttt               <  2:1205471/91‑1 (MQ=255)
                                                                             taaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCtttt               >  1:1361001/1‑91 (MQ=255)
                                                                             taaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCtttt               >  1:2972605/1‑91 (MQ=255)
                                                                              aaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCttttt              >  2:1633669/1‑91 (MQ=255)
                                                                                 aaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGt           >  2:1404288/1‑91 (MQ=255)
                                                                                 aaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGt           >  2:3735376/1‑91 (MQ=255)
                                                                                  aaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTc          >  2:677363/1‑91 (MQ=255)
                                                                                   agagATGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCa         >  2:2379464/1‑91 (MQ=38)
                                                                                     agaTGACTATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGa       <  2:3453245/91‑1 (MQ=255)
                                                                                     agaTGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGa       <  2:1725703/91‑1 (MQ=255)
                                                                                     agaTGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGa       <  2:1745169/91‑1 (MQ=255)
                                                                                        tGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGaatt    >  2:3574175/1‑89 (MQ=255)
                                                                                         gACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtt   <  2:3751239/91‑1 (MQ=255)
                                                                                         gACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtt   >  2:1311098/1‑91 (MQ=255)
                                                                                         gACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtt   <  2:3002133/91‑1 (MQ=255)
                                                                                         gACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtt   <  2:2694947/91‑1 (MQ=255)
                                                                                          aCAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtcg  <  1:1633669/91‑3 (MQ=255)
                                                                                          aCAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTg  >  1:566320/1‑91 (MQ=255)
                                                                                          |                                                                                          
TGTTTGTTCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTG  >  USA300TCH1516_ALE/1779717‑1779897

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: