Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,850,932 T→A 6.2% H224L (CAT→CTT)  acsA_1 ← Acetyl‑coenzyme A synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,850,9320TA 6.2% 650.6 / 23.7 193H224L (CAT→CTT) acsA_1Acetyl‑coenzyme A synthetase
Reads supporting (aligned to +/- strand):  ref base T (97/84);  new base A (8/4);  total (105/88)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 200 total aligned reads displayed.

GGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTCAAATTCAT  >  USA300TCH1516_ALE/1850842‑1851022
                                                                                          |                                                                                          
ggATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACAt                                                                                            <  2:1012846/91‑1 (MQ=255)
ggATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACAt                                                                                            >  1:3345913/1‑91 (MQ=255)
  aTCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCATGTTCAACAAGc                                                                                          <  1:939070/91‑1 (MQ=255)
   tCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCa                                                                                         >  2:1157671/1‑91 (MQ=255)
   tCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCATGTTCAACAAGCa                                                                                         <  2:204127/91‑1 (MQ=255)
    ctgctgTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCaa                                                                                        <  2:865591/91‑1 (MQ=255)
    ctgctgTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCaa                                                                                        <  2:3342026/91‑1 (MQ=255)
      gctgTACAGCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTTCACTAACATTGCTTGTTGAACATGCAATa                                                                                      <  2:3381439/91‑1 (MQ=255)
      gctgTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACAt                                                                                            <  1:1029967/85‑1 (MQ=255)
      gctgTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATa                                                                                      <  1:3490066/91‑1 (MQ=255)
        tgTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATaca                                                                                    >  1:771501/1‑91 (MQ=255)
          tACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACAcc                                                                                  <  2:877356/91‑1 (MQ=255)
          tACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACAcc                                                                                  >  2:1106168/1‑91 (MQ=255)
           acacCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCt                                                                                 >  2:2858416/1‑91 (MQ=255)
            cacCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCtt                                                                                >  1:3684734/1‑91 (MQ=255)
             acCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCttt                                                                               <  2:1451221/91‑1 (MQ=255)
               cAATCAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTCTAg                                                                             >  1:3063876/1‑91 (MQ=255)
               cAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAg                                                                             >  2:1032807/1‑91 (MQ=255)
                aaTAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAgg                                                                            >  1:1451686/1‑91 (MQ=255)
                 aTAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGt                                                                           <  2:3283507/91‑1 (MQ=255)
                 aTAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGt                                                                           >  1:1087694/1‑91 (MQ=255)
                  tAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGtt                                                                          >  1:686502/1‑91 (MQ=255)
                   aaaCATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTg                                                                         >  1:2073646/1‑91 (MQ=255)
                     aCATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGcc                                                                       >  2:1662375/1‑91 (MQ=255)
                     aCATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGcc                                                                       >  2:534391/1‑91 (MQ=255)
                      catcatCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGccc                                                                      >  2:2687368/1‑91 (MQ=255)
                        tcatcTTCTTGTAACCCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGGTGAACATGCAATACACCTTTAGGTTGCCCAg                                                                    <  1:2163768/91‑1 (MQ=255)
                        tcatcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAg                                                                    >  2:572871/1‑91 (MQ=255)
                         catcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGt                                                                   >  2:1072533/1‑91 (MQ=255)
                          atcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTa                                                                  >  2:1593818/1‑91 (MQ=255)
                          atcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTa                                                                  <  1:3216731/91‑1 (MQ=255)
                            cttcttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGa                                                                <  1:3767800/91‑1 (MQ=255)
                            cttcttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGa                                                                >  2:998659/1‑91 (MQ=255)
                              tcttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAAc                                                              >  2:468145/1‑91 (MQ=255)
                               cttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAAcc                                                             <  2:1774799/91‑1 (MQ=255)
                                ttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAAccc                                                            >  2:3838083/1‑90 (MQ=255)
                                   tAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGAt                                                         <  2:3042120/91‑1 (MQ=255)
                                   tAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGAt                                                         >  1:2560436/1‑91 (MQ=255)
                                    aaaTCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATg                                                        <  1:2143779/91‑1 (MQ=255)
                                     aaTCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGt                                                       >  1:668578/1‑91 (MQ=255)
                                     aaTCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGt                                                       >  1:1043251/1‑91 (MQ=255)
                                       tCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGtat                                                     >  1:3415446/1‑91 (MQ=255)
                                       tCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGtat                                                     <  1:1125376/91‑1 (MQ=255)
                                        cTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGtata                                                    <  1:1255461/91‑1 (MQ=255)
                                        cTAATACATAGTTTCCAGAAATATAGTGCACTAACATTGCTTGTGGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGtata                                                    <  2:525472/91‑1 (MQ=255)
                                         tAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATaa                                                   >  1:1553222/1‑91 (MQ=255)
                                          aaTACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAAt                                                  >  1:2033538/1‑91 (MQ=255)
                                          aaTACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAAt                                                  >  1:2396539/1‑91 (MQ=255)
                                           aTACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATg                                                 >  1:2659734/1‑91 (MQ=255)
                                            tACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGt                                                <  1:1593818/91‑1 (MQ=255)
                                            tACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGt                                                <  2:3254970/91‑1 (MQ=255)
                                             aCATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTa                                               >  2:3880493/1‑91 (MQ=255)
                                              cATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTaa                                              <  2:738643/91‑1 (MQ=255)
                                                  ttttCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATc                                          >  2:2283031/1‑91 (MQ=255)
                                                  ttttCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATc                                          >  1:1821859/1‑91 (MQ=255)
                                                   tttCCAGAAATATAGTGCACTAACATTGCTTGTTGATCATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCa                                         <  2:3830217/91‑1 (MQ=255)
                                                   tttCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCa                                         >  1:1289870/1‑91 (MQ=255)
                                                     tCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCaaa                                       <  1:2697510/91‑1 (MQ=255)
                                                       cAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAAcc                                     <  2:245688/91‑1 (MQ=255)
                                                       cAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAAcc                                     <  1:1218860/91‑1 (MQ=255)
                                                          aaaTATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACcatc                                  <  2:2921498/91‑1 (MQ=255)
                                                           aaTATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACcatca                                 >  2:1082854/1‑91 (MQ=255)
                                                              ataGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATcc                              <  2:686443/91‑1 (MQ=255)
                                                              ataGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATcc                              <  2:15039/91‑1 (MQ=255)
                                                                 gTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGAc                           <  2:2492051/91‑1 (MQ=255)
                                                                 gTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGAc                           >  1:818789/1‑91 (MQ=255)
                                                                  tGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACt                          <  1:1513318/91‑1 (MQ=255)
                                                                  tGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACt                          <  1:93803/91‑1 (MQ=255)
                                                                  tGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACt                          <  1:1688405/91‑1 (MQ=255)
                                                                  tGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACt                          <  1:2450704/91‑1 (MQ=255)
                                                                   gCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCatcatc                               <  2:1238332/85‑1 (MQ=255)
                                                                   gCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACtt                         >  1:2107539/1‑91 (MQ=255)
                                                                   gCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACtt                         >  2:1924160/1‑91 (MQ=255)
                                                                   gCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACtt                         >  2:1546917/1‑91 (MQ=255)
                                                                    cACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACttt                        <  1:588850/91‑1 (MQ=255)
                                                                    cACGAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACttt                        <  1:3007001/91‑1 (MQ=255)
                                                                      cTAACATTGCTTGTTGGACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTaa                      <  2:3707695/91‑1 (MQ=255)
                                                                       tAACATTGCTTGTTGAACATGCAATACACGTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAAc                     <  1:3115768/91‑1 (MQ=255)
                                                                       gtgtaTTGCATGTTCAACAAGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAAc                     >  1:3188382/5‑91 (MQ=38)
                                                                         aCATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCa                   <  2:2560436/91‑1 (MQ=255)
                                                                          cATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCAt                  <  1:1791827/91‑1 (MQ=255)
                                                                          cATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCAt                  >  2:3675388/1‑91 (MQ=255)
                                                                           aTTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCAtt                 >  1:400488/1‑91 (MQ=255)
                                                                             tGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCa               >  1:2548764/1‑91 (MQ=255)
                                                                              gCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCaa              >  2:1884805/1‑91 (MQ=255)
                                                                               catgttcAACAAGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAAt             >  1:1604968/3‑91 (MQ=255)
                                                                               catgttcAACAAGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAAt             >  2:186817/3‑91 (MQ=255)
                                                                               cTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAAt             >  2:3476779/1‑91 (MQ=255)
                                                                               cTTGTTGAAAATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACTATTCAAt             >  1:1950414/1‑91 (MQ=255)
                                                                                ttgttgAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATg            >  2:1401045/1‑91 (MQ=255)
                                                                                atgttcAACAAGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATg            >  1:2369956/2‑91 (MQ=255)
                                                                                atgttcAACAAGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATg            >  1:75259/2‑91 (MQ=255)
                                                                                 tgttcAACAAGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGt           >  2:801635/1‑91 (MQ=255)
                                                                                  gttgAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTc          >  1:869182/1‑91 (MQ=255)
                                                                                  gttgAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTc          <  1:105686/91‑1 (MQ=255)
                                                                                   ttgAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTCa         <  1:1716621/91‑1 (MQ=255)
                                                                                    tgAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTCaa        >  1:2860863/1‑91 (MQ=255)
                                                                                     gAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTCaaa       <  1:843487/91‑1 (MQ=255)
                                                                                       aCATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTCAAAtt     >  1:3798494/1‑91 (MQ=255)
                                                                                          tGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTCAAATTCAt  >  2:934120/1‑91 (MQ=255)
                                                                                          |                                                                                          
GGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTCAAATTCAT  >  USA300TCH1516_ALE/1850842‑1851022

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: