Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,685,139 A→T 8.4% intergenic (‑151/+69) 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,1390AT 8.4% 626.9 / 34.1 191intergenic (‑151/+69)USA300HOU_RS08395/rpsUhypothetical protein/30S ribosomal protein S21
Reads supporting (aligned to +/- strand):  ref base A (87/88);  new base T (9/7);  total (96/95)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 197 total aligned reads displayed.

ACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGT  >  USA300TCH1516_ALE/1685050‑1685229
                                                                                         |                                                                                          
aCCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAGATTAGTGTCAg                                                                                           >  1:3515032/1‑91 (MQ=255)
 ccTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGa                                                                                          >  2:3374847/1‑91 (MQ=255)
 ccTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGa                                                                                          >  1:662352/1‑91 (MQ=255)
 ccTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGa                                                                                          <  1:2057163/91‑1 (MQ=255)
  cTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGac                                                                                         <  2:1218254/91‑1 (MQ=255)
    ataGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAATATTGACAACAAAAAAACATCACATATTAGTGTCAGacac                                                                                       >  2:1152518/1‑91 (MQ=255)
     taGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACt                                                                                      >  1:60109/1‑91 (MQ=255)
     taGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACt                                                                                      >  1:2811547/1‑91 (MQ=255)
      aGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTa                                                                                     <  1:2360510/91‑1 (MQ=255)
       gAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTaa                                                                                    >  2:1973016/1‑91 (MQ=255)
       gAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTaa                                                                                    <  1:391593/91‑1 (MQ=255)
       gAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTaa                                                                                    >  2:1972530/1‑91 (MQ=255)
        aaTATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAAc                                                                                   >  1:1876637/1‑91 (MQ=255)
        aaTATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAAc                                                                                   >  2:1329175/1‑91 (MQ=255)
         atatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACt                                                                                  <  1:2198415/91‑1 (MQ=255)
          tatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACtt                                                                                 <  1:3406252/91‑1 (MQ=255)
            tCAATTTGACTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgt                                                                               >  1:1955380/1‑91 (MQ=255)
              aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                             <  2:2686234/91‑1 (MQ=255)
               aTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAt                                                                            >  2:2412331/1‑91 (MQ=255)
               aTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAt                                                                            <  1:1546654/91‑1 (MQ=255)
               aTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAt                                                                            <  1:1973016/91‑1 (MQ=255)
                tttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtg                                                                           >  2:1078387/1‑91 (MQ=255)
                 ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                          >  1:2646473/1‑91 (MQ=255)
                 ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                          >  1:1103001/1‑91 (MQ=255)
                 ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                          >  2:3179935/1‑91 (MQ=255)
                 ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGACACTAATATGTGAtgt                                                                          <  2:1621450/91‑1 (MQ=21)
                  tGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtt                                                                         <  1:2484165/91‑1 (MQ=255)
                    attattTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtttg                                                                       >  2:2201339/1‑91 (MQ=255)
                      tatttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAACAAACATCACAAGTTAGTGTCTGACACTAATATGTGAtgtttttt                                                                     <  2:2372221/91‑4 (MQ=17)
                      tatttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtttgtt                                                                     <  2:100752/91‑1 (MQ=255)
                      tatttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtttgtt                                                                     >  1:424290/1‑91 (MQ=255)
                       atttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATgtttgttt                                                                    <  2:2945187/91‑1 (MQ=255)
                         ttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTta                                                                  <  1:3520036/91‑1 (MQ=255)
                            taAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatat                                                               >  2:3623954/1‑91 (MQ=255)
                             aaaCTATATAACTTTTTTTAATTATATATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatata                                                              >  2:3630635/1‑91 (MQ=255)
                             aaaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatata                                                              >  1:821142/1‑91 (MQ=255)
                             aaaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatata                                                              >  1:1028962/1‑91 (MQ=255)
                               aCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaat                                                            >  2:3689793/1‑91 (MQ=255)
                                cTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatt                                                           <  2:2865002/91‑1 (MQ=255)
                                cTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatt                                                           >  1:959862/1‑91 (MQ=255)
                                  atataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaattaa                                                         <  2:1058499/91‑1 (MQ=255)
                                   tataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaattaat                                                        <  2:1758067/91‑1 (MQ=255)
                                      aaCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATAtt                                                     >  2:1506254/1‑91 (MQ=255)
                                      aaCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATAtt                                                     <  1:3759663/91‑1 (MQ=255)
                                         tttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgt                                                  >  1:2492539/1‑91 (MQ=255)
                                          ttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtg                                                 >  2:677981/1‑91 (MQ=255)
                                           tttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtgt                                                <  2:731120/91‑1 (MQ=255)
                                           tttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtgt                                                <  2:3025846/91‑1 (MQ=255)
                                               taattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAg                                            >  2:2745943/1‑91 (MQ=255)
                                                 attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGactc                                          <  1:1457154/91‑4 (MQ=255)
                                                  ttaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGa                                         >  1:3318758/1‑91 (MQ=255)
                                                  ttaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGa                                         >  1:2342545/1‑91 (MQ=255)
                                                   taaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAg                                        <  2:1537002/91‑1 (MQ=255)
                                                    aaaaaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGt                                       <  2:179959/91‑1 (MQ=255)
                                                         ttGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACa                                  <  2:200169/91‑1 (MQ=255)
                                                         ttGACAACAAAAAAACATCACATACTAGTGCCAGACACTAACTTGTGATGCTCGTTTTATATAATTAATATTCGTGTTGTGGGACTCAAcc                                  >  2:1876374/1‑90 (MQ=255)
                                                          tGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACag                                 <  1:3646827/91‑1 (MQ=255)
                                                           gACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACaga                                >  1:3401858/1‑91 (MQ=255)
                                                            acaacaAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACagag                               <  1:2412331/91‑1 (MQ=255)
                                                             caacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAgg                              >  1:352911/1‑91 (MQ=255)
                                                             caacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAgg                              >  1:351527/1‑91 (MQ=255)
                                                              aacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCGACAGAGGt                             >  2:943697/1‑91 (MQ=255)
                                                               acaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTa                            >  1:2512042/1‑91 (MQ=255)
                                                               acaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTa                            >  2:2852698/1‑91 (MQ=255)
                                                                caaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTat                           <  2:211369/91‑1 (MQ=255)
                                                                 aaaaaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTatt                          >  2:2980108/1‑91 (MQ=255)
                                                                  aacaaaCATCACAAGTTAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTatta                         >  2:3477440/4‑91 (MQ=18)
                                                                  aaaaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTatta                         <  1:2243172/91‑1 (MQ=255)
                                                                    aaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaat                       <  2:3316796/91‑1 (MQ=255)
                                                                     aaaCATCACAACTTAGTGCCACACACTAACTTGTGACGATTGTTCTATATAAGTAATATTCGTGTTGAGGGAGGCAACAGAGGTattaatt                      >  2:3275759/1‑91 (MQ=255)
                                                                       aCATCACAAGTTAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattat                    >  2:2344484/1‑91 (MQ=25)
                                                                        cATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt                   <  1:3401009/91‑1 (MQ=255)
                                                                        cATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt                   <  1:169656/91‑1 (MQ=255)
                                                                        cATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt                   <  1:243586/91‑1 (MQ=255)
                                                                         aTCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTAttt                  <  2:3804284/91‑1 (MQ=255)
                                                                         aTCACAAGTTAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTAttt                  >  1:1345593/1‑91 (MQ=25)
                                                                          tCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTg                 >  1:2149400/1‑91 (MQ=255)
                                                                          tCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTg                 >  2:2854926/1‑91 (MQ=255)
                                                                          tCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTg                 >  1:580621/1‑91 (MQ=255)
                                                                          tCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTg                 >  2:1042001/1‑91 (MQ=255)
                                                                          tCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTg                 >  1:1595153/1‑91 (MQ=255)
                                                                           cacaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGa                >  2:766189/1‑91 (MQ=255)
                                                                            acaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGaa               <  2:2627262/91‑1 (MQ=255)
                                                                            acaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGaa               <  1:3707914/91‑1 (MQ=255)
                                                                             caagttAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAt              >  1:2538557/5‑91 (MQ=37)
                                                                             caTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAt              <  2:2797635/91‑1 (MQ=255)
                                                                              atatTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtt             <  2:512812/91‑1 (MQ=255)
                                                                               tatTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAttt            >  2:3763088/1‑91 (MQ=255)
                                                                               agttAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAttt            >  1:775972/3‑91 (MQ=37)
                                                                                gttAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTTAACAGAGGTATTATTTATTTGAAttta           >  2:2475829/2‑91 (MQ=37)
                                                                                 ttAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttac          <  1:1022267/91‑1 (MQ=255)
                                                                                  tAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacg         <  1:3832727/91‑1 (MQ=255)
                                                                                  tAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacg         >  1:3542647/1‑91 (MQ=255)
                                                                                  tAGCGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacg         <  1:3132631/91‑1 (MQ=255)
                                                                                     tgtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttatgttt      <  1:3782131/91‑1 (MQ=255)
                                                                                     tgtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGGGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgttt      <  1:2897779/91‑1 (MQ=255)
                                                                                      gtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgttta     >  1:2705389/1‑91 (MQ=255)
                                                                                       tCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtttac    <  1:428333/91‑1 (MQ=255)
                                                                                         aGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATttacgtttacgt  >  2:1400111/1‑91 (MQ=255)
                                                                                         aGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATttacgtttacgt  >  1:1833643/1‑91 (MQ=255)
                                                                                         |                                                                                          
ACCTATAGAATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGT  >  USA300TCH1516_ALE/1685050‑1685229

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: