Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,309,800 Δ1 bp 5.9% coding (69/1032 nt) yfiZ_2 ← putative siderophore transport system permease protein YfiZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,309,8000C. 5.9% 1042.2 / 62.0 287coding (69/1032 nt)yfiZ_2putative siderophore transport system permease protein YfiZ
Reads supporting (aligned to +/- strand):  ref base C (143/127);  new base . (9/8);  total (152/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 294 total aligned reads displayed.

GTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTT  >  USA300TCH1516_ALE/2309710‑2309890
                                                                                            |                                                                                           
gTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATAc                                                                                             >  1:1169564/1‑91 (MQ=255)
gTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATAc                                                                                             <  1:4471539/91‑1 (MQ=255)
gTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATAc                                                                                             >  1:4959238/1‑91 (MQ=255)
 tGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACg                                                                                            <  2:2317359/91‑1 (MQ=255)
 tGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACg                                                                                            >  1:5407105/1‑91 (MQ=255)
  gAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGt                                                                                           <  1:4540576/91‑1 (MQ=255)
      gATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATAAACAAATAAAAAGACGCACAACACTCACGTATATAGTGAgt                                                                                       <  2:1565338/91‑1 (MQ=21)
        tatCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAgtgt                                                                                     >  2:5506875/1‑91 (MQ=255)
          tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttg                                                                                   <  1:1742133/91‑1 (MQ=255)
          tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttg                                                                                   >  2:3785778/1‑91 (MQ=255)
          tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttg                                                                                   <  2:4457222/91‑1 (MQ=255)
            gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATAGTGAGTTtgtg                                                                                 <  1:3504933/91‑1 (MQ=25)
            gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATAGTGAGTTtgtg                                                                                 <  2:4564216/91‑1 (MQ=25)
            gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTtgtg                                                                                 >  1:2113409/1‑91 (MQ=255)
            gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTtgtg                                                                                 <  1:4654348/91‑1 (MQ=255)
            gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTtgtg                                                                                 >  1:3043542/1‑91 (MQ=255)
            gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTtgtg                                                                                 >  2:495725/1‑91 (MQ=255)
             cTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGc                                                                                <  2:33122/91‑1 (MQ=255)
             cTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGACGCACAACACTCACGTATATAGTGAGTTTGTGct                                                                               <  1:5298333/91‑2 (MQ=17)
              tAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCg                                                                              <  1:2030781/91‑1 (MQ=255)
               aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGt                                                                             >  2:2494540/1‑91 (MQ=255)
               aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGACGCACAACACTCACGTATATAGTGAGTTTGTGcttt                                                                             <  2:2576281/91‑4 (MQ=12)
                 aattaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtct                                                                           <  2:593095/91‑1 (MQ=255)
                 aattaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtct                                                                           <  1:4368534/91‑1 (MQ=255)
                   ttaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtcttt                                                                         <  2:1228664/91‑1 (MQ=255)
                     aattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtcttttt                                                                       >  1:3667384/1‑91 (MQ=255)
                      cttttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGCGTGTTGTTAGtctttttt                                                                      >  2:2138588/2‑90 (MQ=255)
                      attTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAATGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttc                                                                      <  1:4251260/91‑1 (MQ=255)
                      attTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttc                                                                      >  1:1070597/1‑91 (MQ=255)
                       ttttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttct                                                                     >  1:406840/1‑91 (MQ=255)
                        tttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttctt                                                                    >  1:689134/1‑91 (MQ=255)
                         ttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttcttt                                                                   >  1:819351/1‑91 (MQ=255)
                          tCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttctttt                                                                  <  2:3175863/91‑1 (MQ=255)
                          tCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttctttt                                                                  >  2:1986738/1‑91 (MQ=255)
                          tCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttctttt                                                                  <  1:1660563/91‑1 (MQ=255)
                           cGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTa                                                                 >  2:3902351/1‑91 (MQ=255)
                           cGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTa                                                                 >  1:4855573/1‑91 (MQ=255)
                            gaagaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTAt                                                                <  2:117239/91‑1 (MQ=255)
                            gaagaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTAt                                                                >  2:4763275/1‑91 (MQ=255)
                              agaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGa                                                              <  2:4165483/91‑1 (MQ=255)
                                aaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATc                                                            >  1:133396/1‑91 (MQ=255)
                                 aCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCt                                                           <  2:1674588/91‑1 (MQ=255)
                                 aCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCt                                                           <  1:3999171/91‑1 (MQ=255)
                                    aaTCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATg                                                        <  2:129927/91‑1 (MQ=255)
                                    aaTCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATg                                                        <  1:3118356/91‑1 (MQ=255)
                                     aTCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATgg                                                       >  1:2621373/1‑91 (MQ=255)
                                     aTCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATgg                                                       >  2:312263/1‑91 (MQ=255)
                                       cGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCa                                                     <  1:552879/91‑1 (MQ=255)
                                           atatTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGAtt                                                 >  2:806399/1‑91 (MQ=255)
                                                  aaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGTCAGATTGGCTACt                                          >  2:4286004/1‑91 (MQ=255)
                                                  aaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                          <  1:3267959/91‑1 (MQ=255)
                                                  aaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                          >  2:1712996/1‑91 (MQ=255)
                                                   atatatACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACtt                                         >  1:3751041/1‑91 (MQ=255)
                                                    tatataCAAATAAAAAGAAAGCACAATCTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTa                                        <  2:5134129/91‑1 (MQ=255)
                                                     atataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAt                                       <  2:20838/91‑1 (MQ=255)
                                                      tataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAtt                                      <  1:2960700/91‑1 (MQ=255)
                                                       ataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAttt                                     >  2:1375752/1‑91 (MQ=255)
                                                        taCAAATAAAACGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGCCTACTTAtttt                                    >  2:5620159/1‑90 (MQ=255)
                                                         aCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCt                                   <  2:5217145/91‑1 (MQ=255)
                                                          cAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCtt                                  >  1:3607070/1‑91 (MQ=255)
                                                          cAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCtt                                  >  2:4675447/1‑91 (MQ=255)
                                                            aaTAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCtttt                                <  2:610230/91‑1 (MQ=255)
                                                            aaTAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCtttt                                <  1:1535411/91‑1 (MQ=255)
                                                              tAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTtat                              <  2:1865492/91‑1 (MQ=255)
                                                              tAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTtat                              <  2:3410425/91‑1 (MQ=255)
                                                                 aaagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGc                           >  1:3307390/1‑91 (MQ=255)
                                                                  aagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                          >  1:2251989/1‑91 (MQ=255)
                                                                  aagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                          >  1:108254/1‑91 (MQ=255)
                                                                  aagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                          <  2:2892950/91‑1 (MQ=255)
                                                                  aagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                          <  2:217649/91‑1 (MQ=255)
                                                                  aagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                          <  2:1481140/91‑1 (MQ=255)
                                                                    gaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAt                        >  2:2627382/1‑91 (MQ=255)
                                                                    gaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAt                        >  1:3192866/1‑91 (MQ=255)
                                                                    gaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAt                        >  2:4221106/1‑91 (MQ=255)
                                                                     aaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAtt                       <  1:3408039/91‑1 (MQ=255)
                                                                     aaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAtt                       >  2:4072550/1‑91 (MQ=255)
                                                                     aaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAtt                       <  1:2877907/91‑1 (MQ=255)
                                                                      aagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAttt                      >  2:4018884/1‑91 (MQ=255)
                                                                       cgCACAACACTCACGTATATAGTGAGTTTGTGCTTTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAttt                      >  1:1981090/2‑91 (MQ=17)
                                                                       agCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTa                     >  2:1924926/1‑91 (MQ=255)
                                                                       agCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTa                     >  2:3498575/1‑91 (MQ=255)
                                                                        gCACAACACTCACGTATATAGTGAGTTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc                    >  1:3065506/1‑91 (MQ=25)
                                                                         cacaAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACt                   <  1:2529466/91‑1 (MQ=255)
                                                                          acaAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACtt                  >  1:2384523/1‑91 (MQ=255)
                                                                           caAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTc                 >  1:5137472/1‑91 (MQ=255)
                                                                             acaCTCACGTATATAGTGAGTTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCa               >  2:2359091/3‑91 (MQ=37)
                                                                              caCTCACGTATATAGTGAGTTTGTGCGTCTTTTTCTTTTATTATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAc              >  2:4843801/2‑91 (MQ=37)
                                                                               aCTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAc              <  1:4970698/91‑1 (MQ=255)
                                                                                cTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAcc             >  2:2980235/1‑91 (MQ=255)
                                                                                 tCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCt            <  2:3751373/91‑1 (MQ=255)
                                                                                  cACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCtt           <  1:1286204/91‑1 (MQ=255)
                                                                                  cACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCtt           <  1:4754586/91‑1 (MQ=255)
                                                                                    cTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc         >  1:5081254/1‑91 (MQ=255)
                                                                                    cTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc         >  2:5300295/1‑91 (MQ=255)
                                                                                    cTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc         >  2:2727967/1‑91 (MQ=255)
                                                                                    cTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc         >  1:965765/1‑91 (MQ=255)
                                                                                         ataCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACtatt     >  2:4997054/1‑91 (MQ=255)
                                                                                          taCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTg    <  1:57374/91‑1 (MQ=255)
                                                                                          taCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTg    <  2:3427093/91‑1 (MQ=255)
                                                                                           aCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGt   <  1:880140/91‑1 (MQ=255)
                                                                                            |                                                                                           
GTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTT  >  USA300TCH1516_ALE/2309710‑2309890

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: