Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,432,246 A→G 8.9% intergenic (+593/+54) birA_2 → / ← USA300HOU_RS12525 Bifunctional ligase/repressor BirA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,432,2460AG 8.9% 601.6 / 33.9 199intergenic (+593/+54)birA_2/USA300HOU_RS12525Bifunctional ligase/repressor BirA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (88/93);  new base G (5/13);  total (93/106)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 216 total aligned reads displayed.

AATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACGAAAGA  >  USA300TCH1516_ALE/2432156‑2432335
                                                                                          |                                                                                         
aaTCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaatta                                                                                           >  1:1969634/1‑91 (MQ=255)
aaTCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaatta                                                                                           <  1:2607100/91‑1 (MQ=255)
aaTCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaatta                                                                                           <  1:475902/91‑1 (MQ=255)
 aTCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaattaa                                                                                          <  1:2659424/91‑1 (MQ=255)
  tCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                         >  1:165032/1‑91 (MQ=255)
  tCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                         <  1:435218/91‑1 (MQ=255)
   ccTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATa                                                                                        >  2:3035090/1‑91 (MQ=255)
   ccTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCACAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATa                                                                                        <  1:3088009/91‑1 (MQ=255)
     tAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGa                                                                                      >  2:3518365/1‑91 (MQ=255)
     tAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGa                                                                                      >  1:107643/1‑91 (MQ=255)
      aaTACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGaa                                                                                     <  1:2857653/91‑1 (MQ=37)
      aaTACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGaa                                                                                     <  1:979545/91‑1 (MQ=25)
      aaTACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGaa                                                                                     >  2:770791/1‑91 (MQ=255)
      aaTACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGaa                                                                                     >  1:2180554/1‑91 (MQ=255)
        tACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAAt                                                                                   <  1:561678/91‑1 (MQ=25)
        tACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAAt                                                                                   <  1:1809088/91‑1 (MQ=37)
        tACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAAt                                                                                   <  1:513666/91‑1 (MQ=25)
        tACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAAt                                                                                   <  1:597543/91‑1 (MQ=255)
        tACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAAt                                                                                   <  1:278695/91‑1 (MQ=255)
         aCTGATGCGCTAATTTAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAACAGCTGATTTCTATCAAATAATTAATAGAAATc                                                                                  >  1:3306102/1‑91 (MQ=255)
          cTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATCa                                                                                 <  1:1883575/91‑1 (MQ=25)
          cTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCa                                                                                 >  2:2456477/1‑91 (MQ=255)
            gATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGc                                                                               >  1:1867126/1‑91 (MQ=255)
             aTGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCt                                                                              <  1:1959510/91‑1 (MQ=255)
             aTGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCt                                                                              >  1:2374691/1‑91 (MQ=255)
             aTGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCt                                                                              >  2:1392796/1‑91 (MQ=255)
              tGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtt                                                                             <  2:2259375/91‑1 (MQ=255)
              tGCGCTAATTGAATCTTTATTTATTCAGTTCCATAATAAAATGCACATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtt                                                                             >  2:339835/1‑91 (MQ=255)
               gcgcTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttt                                                                            <  2:2354778/91‑1 (MQ=255)
               gcgcTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttt                                                                            <  2:3260602/91‑1 (MQ=255)
               gcgcTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttt                                                                            <  1:3038596/91‑1 (MQ=255)
                  cTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtttttt                                                                         >  1:1222033/1‑91 (MQ=255)
                    aaTTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCa                                                                       >  1:2251267/1‑91 (MQ=255)
                    aaTTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCa                                                                       <  1:458152/91‑1 (MQ=255)
                      ttGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCAtt                                                                     >  1:276066/1‑91 (MQ=255)
                        gAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGc                                                                   >  2:2596420/1‑91 (MQ=255)
                        gAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGc                                                                   <  2:753210/91‑1 (MQ=255)
                         aaTCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGcc                                                                  >  2:1786651/1‑91 (MQ=255)
                         aaTCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGcc                                                                  >  2:896527/1‑91 (MQ=255)
                         aaTCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGcc                                                                  <  2:2780405/91‑1 (MQ=255)
                          aTCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCt                                                                 >  1:2109141/1‑91 (MQ=255)
                          aTCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCt                                                                 <  1:789552/91‑1 (MQ=255)
                           tCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTa                                                                >  2:1902088/1‑91 (MQ=255)
                           tCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTa                                                                <  1:1395120/91‑1 (MQ=255)
                              ttatttatTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTGTCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTaaaa                                                             <  1:611043/91‑1 (MQ=255)
                              ttatttatTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTaaaa                                                             >  1:1324582/1‑91 (MQ=255)
                              ttatttatTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTaaaa                                                             <  1:2516682/91‑1 (MQ=255)
                               tatttattCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTaaaaa                                                            >  2:2911516/1‑91 (MQ=255)
                                atttattCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAAc                                                           >  1:3394389/1‑91 (MQ=255)
                                atttattCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAAc                                                           <  1:978717/91‑1 (MQ=255)
                                   tattCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTa                                                        <  2:1557876/91‑1 (MQ=255)
                                    attCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTaa                                                       >  1:2673554/1‑91 (MQ=255)
                                    attCAGTTTCAGAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGGAATCAGCTTTTTTCATTGCCTAAAAACTTaa                                                       <  2:2966952/91‑1 (MQ=255)
                                     ttCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAAt                                                      >  1:1460255/1‑91 (MQ=255)
                                      tCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATg                                                     >  1:1601539/1‑91 (MQ=255)
                                       cAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAAGTTAATGt                                                    >  1:1951409/1‑91 (MQ=255)
                                       cAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGt                                                    >  1:1244398/1‑91 (MQ=255)
                                         gTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTcc                                                  <  2:1295565/91‑1 (MQ=255)
                                         gTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTcc                                                  <  2:1345564/91‑1 (MQ=255)
                                            tCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCACCTTTTTTCATTGCCTAAAAACTTAATGTCCCGa                                               >  1:3177852/1‑91 (MQ=255)
                                             cATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGAc                                              >  1:97271/1‑91 (MQ=255)
                                             cATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGAc                                              >  1:2196528/1‑91 (MQ=255)
                                              ataataAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGAcc                                             <  2:1601539/91‑1 (MQ=255)
                                               taataaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACct                                            <  1:2948759/91‑1 (MQ=255)
                                                aataaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACctc                                           <  2:629162/91‑1 (MQ=255)
                                                aataaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACctc                                           <  1:2491720/91‑1 (MQ=255)
                                                 ataaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACctct                                          <  2:2011633/91‑1 (MQ=255)
                                                  taaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCtt                                         <  2:284010/91‑1 (MQ=255)
                                                  taaAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCtt                                         <  1:3123800/91‑1 (MQ=255)
                                                      aTGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGACCCGACCTCTTAATc                                     >  2:2509613/1‑91 (MQ=255)
                                                       tGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCt                                    <  1:2606807/91‑1 (MQ=255)
                                                       tGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCt                                    <  1:3554311/91‑1 (MQ=255)
                                                       tGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCt                                    >  2:243402/1‑91 (MQ=255)
                                                           tatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGc                                >  2:2311772/1‑91 (MQ=255)
                                                                aaaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATaaa                           <  1:3035090/91‑1 (MQ=255)
                                                                 aaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAAt                          <  1:1702966/91‑1 (MQ=255)
                                                                 aaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAAt                          <  1:292287/91‑1 (MQ=255)
                                                                 aaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAAt                          <  2:1867126/91‑1 (MQ=255)
                                                                  aaaGCTGATTTCTATTAATTATTTGATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATa                         >  1:1808353/1‑91 (MQ=21)
                                                                   aaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATAc                        <  2:1022132/91‑1 (MQ=255)
                                                                   aaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATAc                        <  2:2767354/91‑1 (MQ=255)
                                                                   aaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATAc                        >  2:2865946/1‑91 (MQ=255)
                                                                    aGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACt                       >  1:2511678/1‑91 (MQ=255)
                                                                     gCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACtt                      <  2:1117623/91‑1 (MQ=255)
                                                                         aTTTCTATTAATTATTTGATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttatta                  >  2:3100402/1‑91 (MQ=25)
                                                                         aTTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttatta                  >  2:3035636/1‑91 (MQ=255)
                                                                         aTTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACttatta                  <  1:2439322/91‑1 (MQ=255)
                                                                          tttCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTAc                 <  1:1765216/91‑1 (MQ=255)
                                                                          tttCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTAc                 >  2:267198/1‑91 (MQ=255)
                                                                          tttCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTAc                 >  2:2136130/1‑91 (MQ=255)
                                                                           ttCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACt                <  2:3719665/91‑1 (MQ=255)
                                                                              tATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGat             <  1:2604366/91‑1 (MQ=255)
                                                                              tATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGat             <  2:3321288/91‑1 (MQ=255)
                                                                              tATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTAACTACGCATAAATACTTATTACTGat             >  1:3620887/1‑91 (MQ=255)
                                                                               aTCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGata            >  1:1110728/1‑91 (MQ=255)
                                                                                tCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGatat           <  1:3816708/91‑1 (MQ=255)
                                                                                 cAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTACAAACTTAATGTCCCGACCTCTTTATCTGCGCATAAATACTTATTACTCAtata          >  1:3612916/1‑91 (MQ=255)
                                                                                  aaaTAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATaa         >  1:3795709/1‑91 (MQ=255)
                                                                                     taattaatAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACga      <  2:2627724/91‑1 (MQ=255)
                                                                                      aattaatAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACgaa     >  1:2184408/1‑91 (MQ=255)
                                                                                          |                                                                                         
AATCCTAATACTGATGCGCTAATTGAATCTTTATTTATTCAGTTTCATAATAAAATGCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTCATTGCCTAAAAACTTAATGTCCCGACCTCTTTATCTACGCATAAATACTTATTACTGATATAACGAAAGA  >  USA300TCH1516_ALE/2432156‑2432335

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: