Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 3,343 Δ1 bp 5.5% intergenic (+27/‑363) dnaN → / → USA300HOU_RS00015 DNA polymerase III subunit beta/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE3,3430G. 5.5% 769.4 / 35.0 219intergenic (+27/‑363)dnaN/USA300HOU_RS00015DNA polymerase III subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (113/94);  new base . (8/4);  total (121/98)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 236 total aligned reads displayed.

TATTCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATA  >  USA300TCH1516_ALE/3253‑3433
                                                                                          |                                                                                           
tATTCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTg                                                                                             <  1:1465210/91‑1 (MQ=255)
  ttCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGat                                                                                           <  1:3585958/91‑1 (MQ=255)
        aaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGattaatta                                                                                     >  1:232567/1‑91 (MQ=255)
         aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGAttaattaa                                                                                    <  1:2384760/91‑1 (MQ=255)
         aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaat                                                                                   <  1:1070094/91‑2 (MQ=39)
         aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaat                                                                                   <  1:3746006/91‑2 (MQ=39)
           cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaaa                                                                                 >  1:1564962/1‑91 (MQ=255)
           cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaaa                                                                                 >  2:2258715/1‑91 (MQ=255)
           cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                 <  2:1734203/91‑4 (MQ=38)
           cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                 <  1:2103671/91‑4 (MQ=38)
            aaaaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAAc                                                                                >  2:3017895/1‑91 (MQ=255)
              aaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGt                                                                              >  2:143917/1‑91 (MQ=255)
                ggTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCa                                                                            >  2:2150731/1‑91 (MQ=255)
                ggTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCa                                                                            >  2:2241031/1‑91 (MQ=255)
                 gTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCAt                                                                           >  1:2539503/1‑91 (MQ=255)
                   gacgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATcc                                                                         >  2:1545437/1‑91 (MQ=255)
                    acgacTCGGTAACGCAATAAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                        <  2:1173415/91‑1 (MQ=255)
                     cgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCtt                                                                       >  2:3559367/1‑91 (MQ=255)
                     cgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCtt                                                                       <  2:3666082/91‑1 (MQ=255)
                     cgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCtt                                                                       >  2:1692625/1‑91 (MQ=255)
                      gacTCGGTAACGCAATTAATTTTACCAATCAGACCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCttt                                                                      <  1:1179703/91‑1 (MQ=255)
                       acTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTa                                                                     >  2:1023495/1‑91 (MQ=255)
                       acTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTa                                                                     <  1:34048/91‑1 (MQ=255)
                          cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                  >  2:780857/1‑91 (MQ=255)
                           ggTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTTTAATTAATCACGTCATCCTTTAttta                                                                 <  2:1114343/91‑2 (MQ=21)
                               aCGCAATTACTTCTGCCACTCAGATCTTTCGAACACTAAATATACATGAAGGATGACGTGATCAATTAAAGCGTCATCCTTTATTTTTTgg                                                             <  2:1670221/91‑1 (MQ=255)
                               aCGCAATTAATTTTACGAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                             <  2:3707865/91‑1 (MQ=255)
                               aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                             <  2:1306967/91‑1 (MQ=255)
                               aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                             >  2:3725768/1‑91 (MQ=255)
                                cGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGc                                                            <  1:1949216/91‑1 (MQ=255)
                                cGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGc                                                            <  2:2938456/91‑1 (MQ=255)
                                 gCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCa                                                           >  2:1807796/1‑91 (MQ=255)
                                  cAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaa                                                          >  2:1568593/1‑91 (MQ=255)
                                  cAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaa                                                          >  1:1130296/1‑91 (MQ=255)
                                  cAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaa                                                          >  1:3827644/1‑91 (MQ=255)
                                   aattaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaa                                                         <  1:2566207/91‑1 (MQ=255)
                                    attaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaa                                                        >  1:111750/1‑91 (MQ=255)
                                    attaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaa                                                        >  2:887060/1‑91 (MQ=255)
                                     ttaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaaa                                                       >  1:1114657/1‑91 (MQ=255)
                                        attTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaataa                                                    >  1:2641630/1‑91 (MQ=255)
                                          tttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAAtt                                                  <  2:382896/91‑1 (MQ=255)
                                          tttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAAtt                                                  >  1:2218954/1‑91 (MQ=255)
                                          tttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAAtt                                                  >  1:2885187/1‑91 (MQ=255)
                                           ttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTc                                                 <  1:1805511/91‑1 (MQ=255)
                                             aCCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCtt                                                                       <  1:2822517/67‑1 (MQ=255)
                                             aCCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTa                                               <  2:490857/91‑1 (MQ=255)
                                              ccAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAg                                              <  2:921760/91‑1 (MQ=255)
                                              ccAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAg                                              >  2:1117063/1‑91 (MQ=255)
                                               cAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAgg                                             <  2:2170527/91‑1 (MQ=255)
                                                 aTCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTg                                           >  2:440152/1‑91 (MQ=255)
                                                  tCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGc                                          >  1:1195786/1‑91 (MQ=255)
                                                      aaCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTa                                               <  1:3817704/82‑1 (MQ=255)
                                                      aaCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTa                                               >  2:1859494/1‑82 (MQ=255)
                                                       aCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATg                                     <  1:878075/91‑1 (MQ=255)
                                                       aCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATg                                     <  1:1245052/91‑1 (MQ=255)
                                                       aCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATg                                     >  2:2836448/1‑91 (MQ=255)
                                                          tACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaa                                  >  1:2692300/1‑91 (MQ=255)
                                                           aCTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaa                                 <  2:2072805/91‑1 (MQ=255)
                                                           aCTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaa                                 <  1:2481248/91‑1 (MQ=255)
                                                            cTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaaa                                <  2:2398153/91‑1 (MQ=255)
                                                            cTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaaa                                <  2:2435892/91‑1 (MQ=255)
                                                            cTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaaa                                <  1:187276/91‑1 (MQ=255)
                                                            cTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaaa                                <  1:1568647/91‑1 (MQ=255)
                                                            cTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaaa                                >  1:315520/1‑91 (MQ=255)
                                                              aaaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaata                              <  2:325331/91‑1 (MQ=255)
                                                               aaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataa                             >  1:434629/1‑91 (MQ=255)
                                                               aaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataa                             <  1:2012853/91‑1 (MQ=255)
                                                                aaataaatATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataaa                            >  1:1335029/1‑91 (MQ=255)
                                                                aaataaatATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataaa                            >  2:3774136/1‑91 (MQ=255)
                                                                aaataaatATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATATTTTGGCAAAAATAATTCTAGAAGCGAAAGTAaaataaa                            >  2:3376461/1‑91 (MQ=255)
                                                                 aataaataTAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataaat                           >  1:1443669/1‑91 (MQ=255)
                                                                  ataaatataaatAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAAtt                          >  1:2058033/1‑91 (MQ=255)
                                                                  ataaatataaatAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAAtt                          >  1:61621/1‑91 (MQ=255)
                                                                    aaatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTg                        <  1:2948453/91‑1 (MQ=255)
                                                                     aatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTgg                       <  2:614226/91‑1 (MQ=255)
                                                                     aatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTgg                       <  1:344949/91‑1 (MQ=255)
                                                                     aatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTgg                       >  2:857803/1‑91 (MQ=255)
                                                                      atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaa                                  >  1:2864121/1‑79 (MQ=255)
                                                                       tataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgca                     <  2:3799819/91‑1 (MQ=255)
                                                                        ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                    >  1:1973675/1‑91 (MQ=255)
                                                                        ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                    <  2:1345408/91‑1 (MQ=255)
                                                                        ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                    >  2:3636107/1‑91 (MQ=255)
                                                                        ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                    <  2:2219053/91‑1 (MQ=255)
                                                                           aataaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAt                 <  1:666641/91‑1 (MQ=255)
                                                                           aataaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAt                 >  2:2704723/1‑91 (MQ=255)
                                                                            ataaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAATATAAATTTGGCAGCAtt                <  1:1732430/91‑1 (MQ=255)
                                                                             taaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATCTAAAATAAATTTGGCACCAttt               >  2:3136154/1‑91 (MQ=255)
                                                                              aaaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAtttt              >  1:2263683/1‑91 (MQ=21)
                                                                               aaGGATGACGTTTTAATTAATCATGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa             >  2:1600033/1‑91 (MQ=25)
                                                                               aaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa             >  1:387791/1‑91 (MQ=25)
                                                                                aGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaa            <  2:909137/91‑1 (MQ=255)
                                                                                aGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaa            <  2:2549406/91‑1 (MQ=255)
                                                                                aGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaa            <  1:2222238/91‑1 (MQ=255)
                                                                                aGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaa            >  2:272411/1‑91 (MQ=255)
                                                                                    tGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAg        <  1:848609/91‑1 (MQ=255)
                                                                                     gACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGc       >  2:1163174/1‑91 (MQ=255)
                                                                                     gACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGc       >  2:1162925/1‑91 (MQ=255)
                                                                                     gACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGc       >  1:2720441/1‑91 (MQ=255)
                                                                                          gATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATa  >  1:3067461/1‑91 (MQ=255)
                                                                                          gATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATa  <  1:3698802/91‑1 (MQ=255)
                                                                                          |                                                                                           
TATTCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATA  >  USA300TCH1516_ALE/3253‑3433

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: