Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,685,148 C→T 6.3% intergenic (‑160/+60) USA300HOU_RS08395 ← / ← rpsU hypothetical protein/30S ribosomal protein S21

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,685,1480CT 6.3% 350.8 / 7.9 106intergenic (‑160/+60)USA300HOU_RS08395/rpsUhypothetical protein/30S ribosomal protein S21
Reads supporting (aligned to +/- strand):  ref base C (42/57);  new base T (5/2);  total (47/59)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

AATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCT  >  USA300TCH1516_ALE/1685058‑1685238
                                                                                          |                                                                                          
aaTATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAAc                                                                                            <  1:2265386/91‑1 (MQ=255)
 atatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACt                                                                                           <  1:174728/91‑1 (MQ=255)
  tatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACtt                                                                                          >  1:3549548/1‑91 (MQ=255)
  tatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACtt                                                                                          >  1:1477040/1‑91 (MQ=255)
  tatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACtt                                                                                          <  1:1470563/91‑1 (MQ=255)
  tatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACtt                                                                                          <  1:3405482/91‑1 (MQ=255)
   atCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtg                                                                                         >  1:2211800/1‑91 (MQ=255)
    tCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgt                                                                                        <  1:3255911/91‑1 (MQ=255)
    tCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgt                                                                                        >  1:952626/1‑91 (MQ=255)
     cAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgtg                                                                                       <  1:2405671/91‑1 (MQ=255)
     cAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgtg                                                                                       <  1:2405683/91‑1 (MQ=255)
      aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                      <  1:2363161/91‑1 (MQ=255)
      aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                      <  1:1320412/91‑1 (MQ=255)
         ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                                   <  1:3157578/91‑1 (MQ=255)
         ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                                   <  1:3158172/91‑1 (MQ=255)
          tGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtt                                                                                  <  1:3135036/91‑1 (MQ=255)
            attattTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGACACTAATATGTGAtgtttt                                                                                <  1:606586/91‑2 (MQ=21)
              tatttatAATCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAagtttgtt                                                                              >  1:2079173/1‑91 (MQ=255)
              tatttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtttgtt                                                                              >  1:268289/1‑91 (MQ=255)
              tatttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtttgtt                                                                              <  1:1511555/91‑1 (MQ=255)
               atttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATgtttgttt                                                                             >  1:3046366/1‑91 (MQ=255)
                tttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGtttt                                                                            >  1:249153/1‑91 (MQ=255)
                tttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGtttt                                                                            >  1:1373046/1‑91 (MQ=255)
                  tataAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtat                                                                          <  1:3949025/91‑1 (MQ=255)
                    taAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatat                                                                        >  1:258224/1‑91 (MQ=255)
                     aaaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACAGTAACTTGTGATGTTTGTTTtatata                                                                       >  1:3853081/1‑91 (MQ=255)
                     aaaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatata                                                                       >  1:3558969/1‑91 (MQ=255)
                      aaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaa                                                                      >  1:161365/1‑91 (MQ=255)
                      aaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaa                                                                      <  1:1977798/91‑1 (MQ=255)
                        cTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatt                                                                    >  1:3526980/1‑91 (MQ=255)
                         tatataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatta                                                                   >  1:1747358/1‑91 (MQ=255)
                         tatataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatta                                                                   >  1:3241507/1‑91 (MQ=255)
                              aaCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATAtt                                                              <  1:2506881/91‑1 (MQ=255)
                               aCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTc                                                             <  1:4123010/91‑1 (MQ=255)
                               aCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTc                                                             >  1:696906/1‑91 (MQ=255)
                               aCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTc                                                             <  1:295790/91‑1 (MQ=255)
                                 tttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgt                                                           <  1:2895298/91‑1 (MQ=255)
                                 tttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgt                                                           <  1:132180/91‑1 (MQ=255)
                                  ttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtg                                                          <  1:2286181/91‑1 (MQ=255)
                                       taattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAg                                                     >  1:2735437/1‑91 (MQ=255)
                                        aattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAgg                                                    <  1:1861458/91‑1 (MQ=255)
                                        aattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAgg                                                    <  1:1861415/91‑1 (MQ=255)
                                        aattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAgg                                                    >  1:2750862/1‑91 (MQ=255)
                                         attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAATTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGactc                                                   <  1:3588420/91‑4 (MQ=255)
                                         attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGactc                                                   <  1:547922/91‑4 (MQ=255)
                                         attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGactc                                                   <  1:2423734/91‑4 (MQ=255)
                                         attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGactc                                                   <  1:549862/91‑4 (MQ=255)
                                         attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                                   <  1:919776/91‑1 (MQ=255)
                                         attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                                   <  1:917238/91‑1 (MQ=255)
                                          ttaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGa                                                  >  1:2481924/1‑91 (MQ=255)
                                           taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:3359035/91‑1 (MQ=255)
                                           taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:3767394/91‑1 (MQ=255)
                                           taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:1396111/91‑1 (MQ=255)
                                           taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:1969832/91‑1 (MQ=255)
                                           taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:4084396/91‑1 (MQ=255)
                                           taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:1945390/91‑1 (MQ=255)
                                           taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:1831974/91‑1 (MQ=255)
                                             aaaaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTc                                               <  1:1615717/91‑1 (MQ=255)
                                               aaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCaa                                             <  1:2989565/91‑1 (MQ=255)
                                                 ttGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACa                                           <  1:1849190/91‑1 (MQ=255)
                                                   gACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACaga                                         >  1:2489626/1‑91 (MQ=255)
                                                   gACAACAAAAAAACATCACATATTACTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACaga                                         <  1:3845626/91‑1 (MQ=255)
                                                    acaacaAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACagag                                        >  1:1366303/1‑91 (MQ=255)
                                                     caacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAgg                                       >  1:2230887/1‑91 (MQ=255)
                                                      aacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGt                                      <  1:2439611/91‑1 (MQ=255)
                                                      aacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGt                                      <  1:2189850/91‑1 (MQ=255)
                                                       acaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTa                                     <  1:1214900/91‑1 (MQ=255)
                                                        caaAAAAACATCGCATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTat                                    <  1:3148607/91‑1 (MQ=255)
                                                        caaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTat                                    <  1:3148722/91‑1 (MQ=255)
                                                           aaaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACATAGGTagtaa                                 >  1:3732698/1‑91 (MQ=255)
                                                           aaaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaa                                 <  1:973544/91‑1 (MQ=255)
                                                           aaaaaCATCACATATTAGTGTCAGACACTAACATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGGCAACAGAGGTattaa                                 <  1:3536771/91‑1 (MQ=255)
                                                            aaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaat                                >  1:3407968/1‑91 (MQ=255)
                                                            aaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaat                                >  1:3243117/1‑91 (MQ=255)
                                                              aaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaatta                              <  1:3566546/91‑1 (MQ=255)
                                                              aaCATCACAAGTTAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaatta                              >  1:1371798/1‑91 (MQ=21)
                                                              aaCATCACAAGTTAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaatta                              >  1:1372666/1‑91 (MQ=21)
                                                                   cacaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGa                         >  1:3840833/1‑91 (MQ=255)
                                                                    acaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGaa                        >  1:402342/1‑91 (MQ=255)
                                                                    acaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGaa                        >  1:2208282/1‑91 (MQ=255)
                                                                     caagttAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAt                       >  1:1068581/5‑91 (MQ=37)
                                                                      aagttAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtt                      >  1:1051079/4‑91 (MQ=37)
                                                                       tatTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAttt                     >  1:2912638/1‑91 (MQ=255)
                                                                       tatTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAttt                     >  1:2912750/1‑91 (MQ=255)
                                                                          tAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacg                  <  1:563751/91‑1 (MQ=255)
                                                                          tAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacg                  <  1:1087670/91‑1 (MQ=255)
                                                                          tAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacg                  <  1:870196/91‑1 (MQ=255)
                                                                          tAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacg                  <  1:564389/91‑1 (MQ=255)
                                                                            gtgtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtt                >  1:3105576/1‑91 (MQ=255)
                                                                             tgtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgttt               >  1:2289535/1‑91 (MQ=255)
                                                                              gtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgttta              >  1:1435331/1‑91 (MQ=255)
                                                                                  gACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTg          <  1:1947939/91‑1 (MQ=255)
                                                                                     acTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAg       >  1:1504819/1‑91 (MQ=255)
                                                                                      cTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTATGTGCAGc      >  1:2081679/1‑91 (MQ=255)
                                                                                      cTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGc      <  1:2303661/91‑1 (MQ=255)
                                                                                        aaCTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCtt    >  1:1658928/1‑91 (MQ=255)
                                                                                        aaCTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCtt    >  1:1661685/1‑91 (MQ=255)
                                                                                          cTTGTGATGTTTGTTTTATATAATTACGGTTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGCATTTACGTTTACGTGCAGCTTCt  <  1:931494/91‑1 (MQ=255)
                                                                                          cTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCt  >  1:2851595/1‑91 (MQ=255)
                                                                                          cTTCTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCt  <  1:3222386/91‑1 (MQ=255)
                                                                                          |                                                                                          
AATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCT  >  USA300TCH1516_ALE/1685058‑1685238

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: