Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,685,149 T→A 6.5% intergenic (‑161/+59) 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,1490TA 6.5% 351.4 / 10.9 104intergenic (‑161/+59)USA300HOU_RS08395/rpsUhypothetical protein/30S ribosomal protein S21
Reads supporting (aligned to +/- strand):  ref base T (41/56);  new base A (5/2);  total (46/58)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 109 total aligned reads displayed.

ATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCT  >  USA300TCH1516_ALE/1685059‑1685238
                                                                                          |                                                                                         
atatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACt                                                                                           <  1:174728/91‑1 (MQ=255)
 tatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACtt                                                                                          <  1:3405482/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)
  atCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtg                                                                                         >  1:2211800/1‑91 (MQ=255)
   tCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgt                                                                                        <  1:3255911/91‑1 (MQ=255)
   tCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgt                                                                                        >  1:952626/1‑91 (MQ=255)
    cAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgtg                                                                                       <  1:2405683/91‑1 (MQ=255)
    cAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgtg                                                                                       <  1:2405671/91‑1 (MQ=255)
     aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                      <  1:2363161/91‑1 (MQ=255)
     aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                      <  1:1320412/91‑1 (MQ=255)
        ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                                   <  1:3158172/91‑1 (MQ=255)
        ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                                   <  1:3157578/91‑1 (MQ=255)
        ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                                   <  1:2346219/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:3748480/91‑1 (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:1977798/91‑1 (MQ=255)
                     aaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaa                                                                      >  1:161365/1‑91 (MQ=255)
                       cTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatt                                                                    >  1:3526980/1‑91 (MQ=255)
                        tatataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatta                                                                   >  1:1747358/1‑91 (MQ=255)
                            taACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAatat                                                               <  1:2027035/91‑1 (MQ=255)
                             aaCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATAtt                                                              <  1:2506881/91‑1 (MQ=255)
                              aCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTc                                                             >  1:696906/1‑91 (MQ=255)
                              aCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTc                                                             <  1:295790/91‑1 (MQ=255)
                              aCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTc                                                             <  1:4123010/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:2750862/1‑91 (MQ=255)
                                       aattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAgg                                                    <  1:1861415/91‑1 (MQ=255)
                                       aattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAgg                                                    <  1:1861458/91‑1 (MQ=255)
                                        attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGactc                                                   <  1:547922/91‑4 (MQ=255)
                                        attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGactc                                                   <  1:2423734/91‑4 (MQ=255)
                                        attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                                   <  1:919776/91‑1 (MQ=255)
                                        attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                                   <  1:917238/91‑1 (MQ=255)
                                          taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:3359035/91‑1 (MQ=255)
                                          taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:3767394/91‑1 (MQ=255)
                                          taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:1395753/91‑1 (MQ=255)
                                          taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:1396111/91‑1 (MQ=255)
                                          taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:1969832/91‑1 (MQ=255)
                                          taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:1945390/91‑1 (MQ=255)
                                          taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                                 <  1:4084396/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:3243117/1‑91 (MQ=255)
                                                           aaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaat                                >  1:3407968/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)
                                                               cATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt                            <  1:267979/91‑1 (MQ=255)
                                                                  cacaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGa                         >  1:3840833/1‑91 (MQ=255)
                                                                   acaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGaa                        >  1:402342/1‑91 (MQ=255)
                                                                    caagttAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAt                       >  1:1068581/5‑91 (MQ=37)
                                                                     aagttAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtt                      >  1:1051079/4‑91 (MQ=37)
                                                                     aagttAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtt                      >  1:3088842/4‑91 (MQ=37)
                                                                      tatTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAttt                     >  1:2912638/1‑91 (MQ=255)
                                                                         tAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacg                  <  1:1087670/91‑1 (MQ=255)
                                                                         tAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacg                  <  1:870196/91‑1 (MQ=255)
                                                                         tAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacg                  <  1:564389/91‑1 (MQ=255)
                                                                         tAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacg                  <  1:563751/91‑1 (MQ=255)
                                                                           gtgtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtt                >  1:3105576/1‑91 (MQ=255)
                                                                           gtgtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtt                >  1:553867/1‑91 (MQ=255)
                                                                            tgtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgttt               >  1:2289535/1‑91 (MQ=255)
                                                                                 gACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTg          <  1:1947939/91‑1 (MQ=255)
                                                                                  acacTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGc         >  1:3271341/1‑91 (MQ=255)
                                                                                  acacTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGc         >  1:1690736/1‑91 (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)
                                                                                          |                                                                                         
ATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCT  >  USA300TCH1516_ALE/1685059‑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: