Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,850,927 G→C 6.9% Q226E (CAA→GAA)  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,9270GC 6.9% 648.9 / 31.0 187Q226E (CAA→GAA) acsA_1Acetyl‑coenzyme A synthetase
Reads supporting (aligned to +/- strand):  ref base G (93/81);  new base C (2/11);  total (95/92)
Fisher's exact test for biased strand distribution p-value = 9.22e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 209 total aligned reads displayed.

AACCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTCAAA  >  USA300TCH1516_ALE/1850837‑1851017
                                                                                          |                                                                                          
aaCCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCttgttg                                                                                            <  1:392748/91‑1 (MQ=255)
aaCCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCttgttg                                                                                            >  2:2048256/1‑91 (MQ=255)
 aCCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTTCTTGTTGa                                                                                           >  2:1285858/1‑91 (MQ=255)
 aCCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGa                                                                                           >  1:3871600/1‑91 (MQ=255)
  ccTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGaa                                                                                          >  2:1716621/1‑91 (MQ=255)
     ggATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACAt                                                                                       >  1:3345913/1‑91 (MQ=255)
      gATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCATGTTCAACaag                                                                                      <  1:2591006/91‑3 (MQ=255)
      gATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCATGTTCAACaag                                                                                      <  1:2815149/91‑3 (MQ=255)
      gATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCATGTTCAACaag                                                                                      <  2:539977/91‑3 (MQ=255)
       aTCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTGATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGc                                                                                     <  2:1180588/91‑1 (MQ=255)
       aTCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCATGTTCAACaagc                                                                                     <  2:2282130/91‑4 (MQ=255)
       aTCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCATGTTCAACaagc                                                                                     <  1:29934/91‑4 (MQ=255)
       aTCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCATGTTCAACAAGc                                                                                     <  1:939070/91‑1 (MQ=255)
         ctgctgTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCaa                                                                                   <  2:865591/91‑1 (MQ=255)
         ctgctgTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCaa                                                                                   <  2:3342026/91‑1 (MQ=255)
          tgctgTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAAt                                                                                  <  2:3623478/91‑1 (MQ=255)
           gctgTACAGCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTTCACTAACATTGCTTGTTGAACATGCAATa                                                                                 <  2:3381439/91‑1 (MQ=255)
           gctgTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACAt                                                                                       >  2:1029967/1‑85 (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)
                acacCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCt                                                                            <  1:3728869/91‑1 (MQ=255)
                 cacCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCtt                                                                           <  2:113158/91‑1 (MQ=255)
                  acCAATAAACATCATCTTCTTGTAAATCTAATACATCTTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCttt                                                                          <  2:3836540/91‑1 (MQ=255)
                  acCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCttt                                                                          >  1:489574/1‑91 (MQ=255)
                  acCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCttt                                                                          <  1:1807020/91‑1 (MQ=255)
                     aaTAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAgg                                                                       >  1:2137832/1‑91 (MQ=255)
                     aaTAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAgg                                                                       >  1:1451686/1‑91 (MQ=255)
                     aaTAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAgg                                                                       <  2:3731879/91‑1 (MQ=255)
                      aTAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGt                                                                      <  2:3283507/91‑1 (MQ=255)
                      aTAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGt                                                                      >  1:1087694/1‑91 (MQ=255)
                        aaaCATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTg                                                                    >  1:2073647/1‑91 (MQ=255)
                          aCATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGcc                                                                  >  2:1051672/1‑91 (MQ=255)
                          aCATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGcc                                                                  <  2:1905052/91‑1 (MQ=255)
                          aCATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGcc                                                                  >  2:534378/1‑91 (MQ=255)
                           catcatCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGccc                                                                 >  2:2687368/1‑91 (MQ=255)
                             tcatcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAg                                                               >  1:1980980/1‑91 (MQ=255)
                             tcatcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAg                                                               <  1:2163247/91‑1 (MQ=255)
                             tcatcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAg                                                               >  2:572871/1‑91 (MQ=255)
                               atcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTa                                                             <  1:3216731/91‑1 (MQ=255)
                               atcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTa                                                             >  2:1590833/1‑91 (MQ=255)
                                tcttctTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAg                                                            <  2:2297134/91‑1 (MQ=255)
                                   tcttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAAc                                                         >  2:468145/1‑91 (MQ=255)
                                   tcttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAAc                                                         >  2:3696807/1‑91 (MQ=255)
                                    cttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAAcc                                                        >  2:994107/1‑91 (MQ=255)
                                        tAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGAt                                                    >  2:1882301/1‑91 (MQ=255)
                                        tAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGAt                                                    >  1:2560436/1‑91 (MQ=255)
                                        tAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGAt                                                    <  2:3042120/91‑1 (MQ=255)
                                          aaTCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGt                                                  >  1:1043251/1‑91 (MQ=255)
                                          aaTCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGt                                                  >  1:1712123/1‑91 (MQ=255)
                                          aaTCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGt                                                  >  1:668578/1‑91 (MQ=255)
                                            tCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGtat                                                <  1:1125376/91‑1 (MQ=255)
                                            tCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGtat                                                <  2:161068/91‑1 (MQ=255)
                                            tCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGtat                                                >  1:3415446/1‑91 (MQ=255)
                                             cTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGtata                                               <  1:1255461/91‑1 (MQ=255)
                                              tAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATaa                                              >  1:1553222/1‑91 (MQ=255)
                                               aaTACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAAt                                             >  1:2396539/1‑91 (MQ=255)
                                               aaTACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAAt                                             >  1:2036823/1‑91 (MQ=255)
                                               aaTACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAAt                                             >  2:2474929/1‑91 (MQ=255)
                                                 tACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGt                                           <  2:3254970/91‑1 (MQ=255)
                                                 tACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGt                                           <  1:1593818/91‑1 (MQ=255)
                                                  aCATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTa                                          >  2:3880493/1‑91 (MQ=255)
                                                  aCATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTa                                          <  1:3596829/91‑1 (MQ=255)
                                                   cATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTaa                                         <  2:738643/91‑1 (MQ=255)
                                                   cATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTaa                                         >  2:1442384/1‑91 (MQ=255)
                                                        tttCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCa                                    <  1:1730429/91‑1 (MQ=255)
                                                        tttCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCa                                    >  1:1289870/1‑91 (MQ=255)
                                                          tCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCaaa                                  <  1:2697510/91‑1 (MQ=255)
                                                            cAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAAcc                                <  1:1218860/91‑1 (MQ=255)
                                                              gAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACcat                              >  1:1887128/1‑91 (MQ=255)
                                                                aaTATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACcatca                            >  2:1082854/1‑91 (MQ=255)
                                                                   ataGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATcc                         <  2:15217/91‑1 (MQ=255)
                                                                   ataGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATcc                         <  2:686443/91‑1 (MQ=255)
                                                                      gTGCCCTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGAc                      >  1:3059359/1‑91 (MQ=255)
                                                                      gTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGAc                      <  2:2492069/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:2450704/91‑1 (MQ=255)
                                                                        gCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCatcatc                          <  2:1238332/85‑1 (MQ=255)
                                                                        gCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACtt                    >  1:236178/1‑91 (MQ=255)
                                                                         cACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACttt                   <  2:1822002/91‑1 (MQ=255)
                                                                           cTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTaa                 <  1:726418/91‑1 (MQ=255)
                                                                            gtgtaTTGCATGTTCAACAAGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAAc                >  1:3188382/5‑91 (MQ=38)
                                                                              gtaTTGCATGTTCAACAAGCAATACACCTTTAGGTTGCCCAGTAGACCCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCa              >  1:3097835/3‑91 (MQ=255)
                                                                              aCATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCa              <  2:2560436/91‑1 (MQ=255)
                                                                                aTTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCAtt            >  1:3772822/1‑91 (MQ=255)
                                                                                  tGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCa          >  1:2548764/1‑91 (MQ=255)
                                                                                    cTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAAt        >  2:3476779/1‑91 (MQ=255)
                                                                                      tgttcAACAAGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGt      >  2:801635/1‑91 (MQ=255)
                                                                                       gttgAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTc     <  1:105686/91‑1 (MQ=255)
                                                                                       gttgAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTc     >  1:869182/1‑91 (MQ=255)
                                                                                       gttgAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTc     >  1:98850/1‑91 (MQ=255)
                                                                                       gttgAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTc     <  2:1458571/91‑1 (MQ=255)
                                                                                       gttgAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTc     >  2:2859766/1‑91 (MQ=255)
                                                                                       gttgAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTc     >  2:1984178/1‑91 (MQ=255)
                                                                                         tgAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTCaa   <  1:1633947/91‑1 (MQ=255)
                                                                                         tgAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTCaa   >  1:2860863/1‑91 (MQ=255)
                                                                                          gAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTCaaa  <  1:843487/91‑1 (MQ=255)
                                                                                          |                                                                                          
AACCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATGTCAAA  >  USA300TCH1516_ALE/1850837‑1851017

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: