Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,200,697 A→C 7.3% E171A (GAA→GCA)  divIVA → Septum site‑determining protein DivIVA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,200,6970AC 7.3% 424.9 / 21.3 123E171A (GAA→GCA) divIVASeptum site‑determining protein DivIVA
Reads supporting (aligned to +/- strand):  ref base A (49/65);  new base C (5/4);  total (54/69)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 130 total aligned reads displayed.

ACTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCAATGATTCAGA  >  USA300TCH1516_ALE/1200607‑1200787
                                                                                          |                                                                                          
aCTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGa                                                                                            >  1:522872/1‑91 (MQ=255)
aCTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGa                                                                                            <  1:1301893/91‑1 (MQ=255)
 cTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGaa                                                                                           >  1:4002025/1‑91 (MQ=255)
  ttattaAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGaaa                                                                                          >  1:1591643/1‑91 (MQ=255)
   tattaAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGaaaa                                                                                         <  1:3468876/91‑1 (MQ=255)
    attaAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTTCATGCAAAt                                                                                        <  1:713491/91‑1 (MQ=255)
     ttaAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATg                                                                                       >  1:605879/1‑91 (MQ=255)
     ttaAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATg                                                                                       <  1:2758659/91‑1 (MQ=255)
        aaaaaCGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTTCATGCAAATGAtg                                                                                    <  1:1747684/91‑2 (MQ=255)
        aaaaaCGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTTCATGCAAATGAtg                                                                                    <  1:1747667/91‑2 (MQ=255)
        aaaaaCGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTTCATGCAAATGAtg                                                                                    <  1:1445842/91‑2 (MQ=255)
        aaaaaCGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGAtt                                                                                    <  1:921170/91‑1 (MQ=255)
        aaaaaCGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGAtt                                                                                    <  1:634599/91‑1 (MQ=255)
         aaaaCGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGAttt                                                                                   >  1:2417941/1‑91 (MQ=255)
             cGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAg                                                                               <  1:4043700/91‑1 (MQ=255)
             cGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAg                                                                               <  1:4043683/91‑1 (MQ=255)
             cGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAg                                                                               >  1:2606589/1‑91 (MQ=255)
             cGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAg                                                                               <  1:2520897/91‑1 (MQ=255)
              gAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGc                                                                              <  1:134295/91‑1 (MQ=255)
                 gATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAg                                                                           <  1:2388225/91‑1 (MQ=255)
                 gATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCAAg                                                                           >  1:1016248/1‑91 (MQ=255)
                   ttGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGAt                                                                         <  1:4061505/91‑1 (MQ=255)
                   ttGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGAt                                                                         <  1:1796031/91‑1 (MQ=255)
                       gATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAg                                                                     >  1:2593397/1‑91 (MQ=255)
                        aTTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGt                                                                    <  1:1432453/91‑1 (MQ=255)
                        aTTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGt                                                                    <  1:2039056/91‑1 (MQ=255)
                          tACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTg                                                                  <  1:2333501/91‑1 (MQ=255)
                           aCTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTgc                                                                 <  1:3828276/91‑1 (MQ=255)
                               gttgAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTgcagca                                                             >  1:2923127/1‑91 (MQ=255)
                               gttgAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTgcagca                                                             >  1:3591207/1‑91 (MQ=255)
                               gttgAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTgcagca                                                             >  1:2472806/1‑91 (MQ=255)
                                 tgAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCaaa                                                           >  1:3207685/1‑91 (MQ=255)
                                   aaTTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATg                                                         <  1:1250541/91‑1 (MQ=255)
                                     ttATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGca                                                       <  1:1945210/91‑1 (MQ=255)
                                     ttATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGca                                                       <  1:2140140/91‑1 (MQ=255)
                                      tATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGcac                                                      >  1:712563/1‑91 (MQ=255)
                                      tATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGcac                                                      <  1:1695884/91‑1 (MQ=255)
                                      tATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGcac                                                      >  1:2211204/1‑91 (MQ=255)
                                      tATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGcac                                                      <  1:2255840/91‑1 (MQ=255)
                                         gATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACaaa                                                   >  1:2720583/1‑91 (MQ=255)
                                         gATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACaaa                                                   >  1:1691196/1‑91 (MQ=255)
                                          aTTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACaaaa                                                  <  1:3446448/91‑1 (MQ=255)
                                           tttAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAaat                                                 <  1:3713228/91‑1 (MQ=255)
                                            ttAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAaata                                                >  1:1686413/1‑91 (MQ=255)
                                             tAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATTAAAATGATTTAAAGCCAGATGAGGTTGCAGCAAATGCACAAaataa                                               <  1:3288133/91‑1 (MQ=255)
                                               gACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATg                                             >  1:3825351/1‑91 (MQ=255)
                                                aCGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGc                                            <  1:1059581/91‑1 (MQ=255)
                                                  gCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTGAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCAt                                          <  1:3593156/91‑1 (MQ=255)
                                                  gCTGAACAAGTCACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCAt                                          <  1:3762377/91‑1 (MQ=255)
                                                   cTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACACAATAATGCATc                                         <  1:1336797/91‑1 (MQ=255)
                                                   cTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATc                                         >  1:91605/1‑91 (MQ=255)
                                                    tGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCa                                        <  1:1656790/91‑1 (MQ=255)
                                                    tGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCa                                        <  1:3480455/91‑1 (MQ=255)
                                                    tGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCa                                        <  1:2811816/91‑1 (MQ=255)
                                                      aaCAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCaaa                                      <  1:1596905/91‑1 (MQ=255)
                                                      aaCAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCaaa                                      <  1:3953659/91‑1 (MQ=255)
                                                         aaGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATac                                   >  1:2422805/1‑91 (MQ=255)
                                                           gTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATacac                                 >  1:3456468/1‑91 (MQ=255)
                                                           gTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATacac                                 <  1:43227/91‑1 (MQ=255)
                                                             gACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCa                               <  1:3920348/91‑1 (MQ=255)
                                                             gACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCa                               <  1:3920732/91‑1 (MQ=255)
                                                              aCGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAg                              >  1:2097278/1‑91 (MQ=255)
                                                               cGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGa                             >  1:1577791/1‑91 (MQ=255)
                                                               cGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGa                             >  1:2849731/1‑91 (MQ=255)
                                                               cGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGa                             <  1:1750425/91‑1 (MQ=255)
                                                                gCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGAc                            <  1:2242662/91‑1 (MQ=255)
                                                                gCTTGAAAATATTCATCATTTGCATGAAAATGAGTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGAc                            <  1:3357969/91‑1 (MQ=255)
                                                                  ttGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACaa                          >  1:1545156/1‑91 (MQ=255)
                                                                    gAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACaata                        <  1:1813326/91‑1 (MQ=255)
                                                                    gAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACaata                        >  1:898905/1‑91 (MQ=255)
                                                                    gAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACaata                        <  1:1813610/91‑1 (MQ=255)
                                                                     aaaaTATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACaataa                       >  1:3514181/1‑91 (MQ=255)
                                                                     aaaaTATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACaataa                       <  1:1081390/91‑1 (MQ=255)
                                                                       aaTATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATc                     >  1:49719/1‑91 (MQ=255)
                                                                       aaTATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATc                     >  1:3042998/1‑91 (MQ=255)
                                                                        atatTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATca                    <  1:1529627/91‑1 (MQ=255)
                                                                        atatTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATca                    <  1:3654185/91‑1 (MQ=255)
                                                                        atatTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATca                    <  1:3654548/91‑1 (MQ=255)
                                                                        atatTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATca                    <  1:3241369/91‑1 (MQ=255)
                                                                         tatTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATcaa                   <  1:2635767/91‑1 (MQ=255)
                                                                           ttCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATcaaca                 >  1:1153843/1‑91 (MQ=255)
                                                                             catcatTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAAt               <  1:863503/91‑1 (MQ=255)
                                                                             catcatTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAAt               >  1:695517/1‑91 (MQ=255)
                                                                               tcatTTTCATGCAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATcc             >  1:3098322/1‑91 (MQ=255)
                                                                               tcatTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATcc             <  1:2895149/91‑1 (MQ=255)
                                                                               tcatTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATcc             >  1:1304421/1‑91 (MQ=255)
                                                                               tcatTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATcc             >  1:76122/1‑91 (MQ=255)
                                                                                catTTTCATGCAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCa            >  1:1205453/1‑91 (MQ=255)
                                                                                catTTTCATGCAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCa            >  1:574491/1‑91 (MQ=255)
                                                                                catTTTCATGCAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCa            >  1:683416/1‑91 (MQ=255)
                                                                                catTTTCATGCAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCa            >  1:1985501/1‑91 (MQ=255)
                                                                                 atTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCaa           <  1:772278/91‑1 (MQ=255)
                                                                                   ttGTATGAATATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCAATg         <  1:3726861/91‑1 (MQ=255)
                                                                                    tGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCAATGa        >  1:2136842/1‑91 (MQ=255)
                                                                                     gCACGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCAATGAt       <  1:3344508/91‑1 (MQ=255)
                                                                                      cATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCTCAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCAATGAtt      >  1:3590025/1‑91 (MQ=255)
                                                                                       aTGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCAATGATTc     >  1:1443250/1‑91 (MQ=255)
                                                                                        tGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCAATGATTCa    >  1:1207846/1‑91 (MQ=255)
                                                                                         gAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCAATGATTCAg   >  1:2122153/1‑91 (MQ=255)
                                                                                          aaaaTGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCAATGATTCAGa  <  1:1506245/91‑1 (MQ=255)
                                                                                          |                                                                                          
ACTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCAATGATTCAGA  >  USA300TCH1516_ALE/1200607‑1200787

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