Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,200,692 G→T 7.6% L169F (TTG→TTT 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,6920GT 7.6% 448.6 / 22.7 131L169F (TTG→TTTdivIVASeptum site‑determining protein DivIVA
Reads supporting (aligned to +/- strand):  ref base G (51/70);  new base T (5/5);  total (56/75)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 133 total aligned reads displayed.

ATTAGACTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCAATGAT  >  USA300TCH1516_ALE/1200602‑1200782
                                                                                          |                                                                                          
aTTAGACTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTg                                                                                            <  1:4019465/91‑1 (MQ=255)
 ttAGACTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGc                                                                                           >  1:4084579/1‑91 (MQ=255)
   aGACTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCAt                                                                                         >  1:3954716/1‑91 (MQ=255)
    gACTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATg                                                                                        >  1:1051812/1‑91 (MQ=255)
    gACTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATg                                                                                        <  1:1617716/91‑1 (MQ=255)
     aCTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGa                                                                                       <  1:1301893/91‑1 (MQ=255)
     aCTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGa                                                                                       >  1:522872/1‑91 (MQ=255)
      cTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTTCATgca                                                                                      <  1:1172255/91‑3 (MQ=255)
      cTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGaa                                                                                      >  1:4002025/1‑91 (MQ=255)
       ttattaAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGaaa                                                                                     >  1:1591643/1‑91 (MQ=255)
        tattaAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGaaaa                                                                                    <  1:3468835/91‑1 (MQ=255)
         attaAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTTCATGCAAAt                                                                                   <  1:713491/91‑1 (MQ=255)
          ttaAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATg                                                                                  >  1:605879/1‑91 (MQ=255)
             aaaaaCGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTTCATGCAAATGAtg                                                                               <  1:1445842/91‑2 (MQ=255)
             aaaaaCGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTTCATGCAAATGAtg                                                                               <  1:1747684/91‑2 (MQ=255)
             aaaaaCGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGAtt                                                                               <  1:921170/91‑1 (MQ=255)
             aaaaaCGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGAtt                                                                               <  1:634599/91‑1 (MQ=255)
              aaaaCGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGAttt                                                                              <  1:3584340/91‑1 (MQ=255)
              aaaaCGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGAttt                                                                              <  1:2743220/91‑1 (MQ=255)
              aaaaCGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGAttt                                                                              >  1:2417941/1‑91 (MQ=255)
                  cGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAg                                                                          <  1:4043683/91‑1 (MQ=255)
                  cGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAg                                                                          <  1:2520897/91‑1 (MQ=255)
                   gAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGc                                                                         >  1:156278/1‑91 (MQ=255)
                   gAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGc                                                                         <  1:134295/91‑1 (MQ=255)
                      gATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCAAg                                                                      >  1:1016248/1‑91 (MQ=255)
                       aTTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGa                                                                     <  1:649276/91‑1 (MQ=255)
                        ttGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGAt                                                                    <  1:4061505/91‑1 (MQ=255)
                            gATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAg                                                                >  1:2593397/1‑91 (MQ=255)
                             aTTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGt                                                               <  1:2039056/91‑1 (MQ=255)
                             aTTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGt                                                               <  1:1432453/91‑1 (MQ=255)
                             aTTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTCAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGt                                                               <  1:3307534/91‑1 (MQ=255)
                               tACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTg                                                             <  1:2333484/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)
                                      tgAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCaaa                                                      >  1:3207685/1‑91 (MQ=255)
                                          ttATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGca                                                  <  1:2140140/91‑1 (MQ=255)
                                           tATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGcac                                                 >  1:2211204/1‑91 (MQ=255)
                                           tATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGcac                                                 >  1:2601946/1‑91 (MQ=255)
                                           tATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGcac                                                 <  1:1695884/91‑1 (MQ=255)
                                           tATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGcac                                                 <  1:2255840/91‑1 (MQ=255)
                                           tATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGcac                                                 >  1:712563/1‑91 (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)
                                                    gACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATg                                        >  1:3825351/1‑91 (MQ=255)
                                                    gACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATg                                        >  1:2161254/1‑91 (MQ=255)
                                                     aCGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGc                                       <  1:1059581/91‑1 (MQ=255)
                                                       gCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCAGAAAATAATGCAt                                     <  1:3205402/91‑1 (MQ=255)
                                                       gCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTGAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCAt                                     <  1:3593156/91‑1 (MQ=255)
                                                        cTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATc                                    <  1:1337008/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)
                                                         tGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCa                                   <  1:1494611/91‑1 (MQ=255)
                                                          gAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCaa                                  >  1:3738195/1‑91 (MQ=255)
                                                           aaCAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCaaa                                 <  1:3953659/91‑1 (MQ=255)
                                                              aaGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATac                              >  1:3918497/1‑91 (MQ=255)
                                                              aaGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATac                              >  1:2422805/1‑91 (MQ=255)
                                                                gTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATacac                            >  1:3456468/1‑91 (MQ=255)
                                                                  gACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCa                          <  1:3920732/91‑1 (MQ=255)
                                                                  gACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCa                          >  1:2438562/1‑91 (MQ=255)
                                                                   aCGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAg                         >  1:2097278/1‑91 (MQ=255)
                                                                    cGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGa                        >  1:1577791/1‑91 (MQ=255)
                                                                    cGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGa                        <  1:1750425/91‑1 (MQ=255)
                                                                    cGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGa                        >  1:2849731/1‑91 (MQ=255)
                                                                     gCTTGAAAATATTCATCATTTGCATGAAAATGAGTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGAc                       <  1:3357969/91‑1 (MQ=255)
                                                                        tGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACaat                    >  1:1280447/1‑91 (MQ=255)
                                                                         gAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACaata                   >  1:898905/1‑91 (MQ=255)
                                                                         gAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACaata                   <  1:1813610/91‑1 (MQ=255)
                                                                         gAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACaata                   <  1:1813326/91‑1 (MQ=255)
                                                                          aaaaTATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACaataa                  >  1:3514181/1‑91 (MQ=255)
                                                                          aaaaTATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACaataa                  <  1:1773055/91‑1 (MQ=255)
                                                                          aaaaTATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACaataa                  <  1:1081390/91‑1 (MQ=255)
                                                                            aaTATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATc                <  1:1709004/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:3241369/91‑1 (MQ=255)
                                                                              tatTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATcaa              <  1:2635767/91‑1 (MQ=255)
                                                                                ttCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACTCCAGACAATAATcaaca            >  1:1153911/1‑91 (MQ=255)
                                                                                ttCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTCCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATcaaca            >  1:3597711/1‑91 (MQ=255)
                                                                                 tcatcaTTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATcaacaa           >  1:1614778/1‑91 (MQ=255)
                                                                                  catcatTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAAt          >  1:421427/1‑91 (MQ=255)
                                                                                  catcatTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAAt          <  1:863503/91‑1 (MQ=255)
                                                                                  catcatTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAAt          >  1:695517/1‑91 (MQ=255)
                                                                                  catcatTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAAt          <  1:551114/91‑1 (MQ=255)
                                                                                   atcatTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATc         >  1:1858270/1‑91 (MQ=255)
                                                                                    tcatTTTCATGCAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATcc        >  1:3098322/1‑91 (MQ=255)
                                                                                    tcatTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATcc        >  1:1304421/1‑91 (MQ=255)
                                                                                    tcatTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATcc        >  1:76122/1‑91 (MQ=255)
                                                                                    tcatTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATcc        <  1:2895149/91‑1 (MQ=255)
                                                                                     catTTTCATGCAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCa       >  1:1205453/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)
                                                                                          gCACGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCAATGAt  <  1:3344508/91‑1 (MQ=255)
                                                                                          |                                                                                          
ATTAGACTTATTAAAAAACGAAGATTGGGATTACTTGTTGAATTATGATTTAGACGCTGAACAAGTGACGCTTGAAAATATTCATCATTTGCATGAAAATGATTTAAAGCCAGATGAAGTTGCAGCAAATGCACAAAATAATGCATCAAATACACCAGACAATAATCAACAATCCAATGAT  >  USA300TCH1516_ALE/1200602‑1200782

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