Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,850,922 T→A 6.9% Q227H (CAA→CAT 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,9220TA 6.9% 624.0 / 25.5 187Q227H (CAA→CATacsA_1Acetyl‑coenzyme A synthetase
Reads supporting (aligned to +/- strand):  ref base T (93/81);  new base A (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 = 9.99e-01

Only 100 of 207 total aligned reads displayed.

ACCCAACCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATG  >  USA300TCH1516_ALE/1850833‑1851012
                                                                                         |                                                                                          
aCCCAACCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCtt                                                                                           >  1:3736613/1‑91 (MQ=255)
 cccAACCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATGTTTTCCAGAAATATAGTGCACTAACATTGCttg                                                                                          <  1:3335191/91‑1 (MQ=255)
 cccAACCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCttg                                                                                          <  2:232998/91‑1 (MQ=255)
 cccAACCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCttg                                                                                          <  2:3293995/91‑1 (MQ=255)
   cAACCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCttgtt                                                                                        >  1:3510235/1‑91 (MQ=255)
    aaCCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCttgttg                                                                                       <  2:2520793/91‑1 (MQ=255)
     aCCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTTCTTGTTGa                                                                                      >  2:1285858/1‑91 (MQ=255)
     aCCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGa                                                                                      >  1:3871600/1‑91 (MQ=255)
       cTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAAc                                                                                    >  1:965213/1‑91 (MQ=255)
        tGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCATGTTCAACa                                                                                   <  2:352456/91‑1 (MQ=255)
         ggATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCttg                                                                                          >  2:1786291/1‑83 (MQ=255)
         ggATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCttg                                                                                          <  1:1786291/83‑1 (MQ=255)
         ggATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACAt                                                                                  <  1:892/91‑1 (MQ=255)
          gATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCATGTTCAACaag                                                                                 <  1:2591006/91‑3 (MQ=255)
           aTCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCATGTTCAACAAGc                                                                                <  1:939039/91‑1 (MQ=255)
            tCTGCTGTACACCAATAAACATCATCTTCTTGTCAATCTAATACATATTTTCCAGAAATCTAGTGCACTAACATTGCATGTTCAACAAGCa                                                                               <  1:3445819/91‑1 (MQ=255)
            tCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCa                                                                               >  2:1157671/1‑91 (MQ=255)
            tCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCa                                                                               >  2:402685/1‑91 (MQ=255)
            tCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCATGTTCAACAAGCa                                                                               <  2:204127/91‑1 (MQ=255)
             ctgctgTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCaa                                                                              <  2:865693/91‑1 (MQ=255)
              tgctgTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAAt                                                                             <  2:3623478/91‑1 (MQ=255)
               gctgTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACAt                                                                                  <  1:1029967/85‑1 (MQ=255)
               gctgTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACAt                                                                                  >  2:1029967/1‑85 (MQ=255)
                   tACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACAcc                                                                        >  2:1106168/1‑91 (MQ=255)
                    acacCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCt                                                                       <  1:2764682/91‑1 (MQ=255)
                     cacCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCtt                                                                      >  1:3684734/1‑91 (MQ=255)
                      acCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCttt                                                                     <  1:1807020/91‑1 (MQ=255)
                        cAATCAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTCTAg                                                                   >  1:3063876/1‑91 (MQ=255)
                        cAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAg                                                                   >  2:767089/1‑91 (MQ=255)
                         aaTAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAgg                                                                  >  1:1451686/1‑91 (MQ=255)
                            aaaCATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTg                                                               >  1:2073646/1‑91 (MQ=255)
                              aCATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGcc                                                             <  2:1905052/91‑1 (MQ=255)
                              aCATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGcc                                                             >  2:1662375/1‑91 (MQ=255)
                              aCATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGcc                                                             >  2:1051672/1‑91 (MQ=255)
                              aCATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGcc                                                             >  2:534391/1‑91 (MQ=255)
                               catcatCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGccc                                                            >  2:805268/1‑91 (MQ=255)
                                atcatcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCa                                                           >  2:3400633/1‑91 (MQ=255)
                                 tcatcTTCTTGTAACCCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGGTGAACATGCAATACACCTTTAGGTTGCCCAg                                                          <  1:2163768/91‑1 (MQ=255)
                                 tcatcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAg                                                          >  2:679323/1‑91 (MQ=255)
                                 tcatcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAg                                                          >  2:89226/1‑91 (MQ=255)
                                 tcatcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAg                                                          >  2:572057/1‑91 (MQ=255)
                                   atcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTa                                                        <  2:326421/91‑1 (MQ=255)
                                   atcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTa                                                        >  2:1593818/1‑91 (MQ=255)
                                   atcTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTa                                                        <  1:3216731/91‑1 (MQ=255)
                                    tcttctTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAg                                                       <  2:177054/91‑1 (MQ=255)
                                     cttcttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGa                                                      <  2:1207193/91‑1 (MQ=255)
                                      ttcttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGaa                                                     <  2:2681403/91‑1 (MQ=255)
                                       tcttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAAc                                                    >  2:726418/1‑91 (MQ=255)
                                       tcttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAAc                                                    >  2:3696807/1‑91 (MQ=255)
                                        cttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAAcc                                                   >  2:994107/1‑91 (MQ=255)
                                        cttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAAcc                                                   >  1:113211/1‑91 (MQ=255)
                                         ttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAAccc                                                  >  2:3838083/1‑90 (MQ=255)
                                         ttGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCt                                                  >  2:3106358/1‑91 (MQ=255)
                                            tAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGAt                                               >  1:2560436/1‑91 (MQ=255)
                                            tAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGAt                                               >  2:1882301/1‑91 (MQ=255)
                                             aaaTCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATg                                              <  1:2143779/91‑1 (MQ=255)
                                              aaTCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGt                                             >  1:1043251/1‑91 (MQ=255)
                                              aaTCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGt                                             >  1:668578/1‑91 (MQ=255)
                                              aaTCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGt                                             >  1:1712123/1‑91 (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                                      <  1:824336/91‑1 (MQ=255)
                                                      aCATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTa                                     >  2:3880476/1‑91 (MQ=255)
                                                      aCATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTa                                     <  1:3596829/91‑1 (MQ=255)
                                                       cATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTaa                                    >  2:1442384/1‑91 (MQ=255)
                                                       cATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTaa                                    <  1:1106168/91‑1 (MQ=255)
                                                          atTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAAt                                 >  1:781889/1‑91 (MQ=255)
                                                           ttttCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATc                                >  2:2283031/1‑91 (MQ=255)
                                                            tttCCAGAAATATAGTGCACTAACATTGCTTGTTGATCATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCa                               <  2:3830217/91‑1 (MQ=255)
                                                             ttCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCaa                              <  2:3345913/91‑1 (MQ=255)
                                                                  gAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACcat                         >  1:1887128/1‑91 (MQ=255)
                                                                   aaaTATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACcatc                        <  2:2921498/91‑1 (MQ=255)
                                                                    aaTATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACcatca                       >  2:1081935/1‑91 (MQ=255)
                                                                    aaTATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACcatca                       >  2:2773914/1‑91 (MQ=255)
                                                                      tataGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCatcatc                     <  1:1285578/91‑1 (MQ=255)
                                                                       ataGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATcc                    <  2:15039/91‑1 (MQ=255)
                                                                        taGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCg                   >  2:3325872/1‑91 (MQ=255)
                                                                          gTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGAc                 <  2:2492051/91‑1 (MQ=255)
                                                                           tGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACt                <  2:2137832/91‑1 (MQ=255)
                                                                           tGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACt                <  1:2450704/91‑1 (MQ=255)
                                                                           tGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACt                <  1:1513318/91‑1 (MQ=255)
                                                                            gCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCatcatc                     >  1:1238332/1‑85 (MQ=255)
                                                                            gCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACtt               >  1:2107539/1‑91 (MQ=255)
                                                                            gCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACtt               >  1:236178/1‑91 (MQ=255)
                                                                            gCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACtt               >  1:3394959/1‑91 (MQ=255)
                                                                            gCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACtt               >  2:1546917/1‑91 (MQ=255)
                                                                            gCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGCTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACtt               >  1:2831041/1‑91 (MQ=255)
                                                                             cACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACttt              <  1:588850/91‑1 (MQ=255)
                                                                             cACGAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACttt              <  1:3007001/91‑1 (MQ=255)
                                                                               cTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTaa            <  1:726418/91‑1 (MQ=255)
                                                                                tAACATTGCTTGTTGAACATGCAATACACGTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAAc           <  1:3115768/91‑1 (MQ=255)
                                                                                gtgtaTTGCATGTTCAACAAGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAAc           >  1:3188382/5‑91 (MQ=38)
                                                                                  aCATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCa         <  1:1072533/91‑1 (MQ=255)
                                                                                    aTTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCAtt       >  1:400488/1‑91 (MQ=255)
                                                                                    aTTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCAtt       >  1:3772822/1‑91 (MQ=255)
                                                                                      tGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCa     >  1:2548764/1‑91 (MQ=255)
                                                                                         |                                                                                          
ACCCAACCTGGATCTGCTGTACACCAATAAACATCATCTTCTTGTAAATCTAATACATATTTTCCAGAAATATAGTGCACTAACATTGCTTGTTGAACATGCAATACACCTTTAGGTTGCCCAGTAGAACCTGATGTATAATGTAAAATCAAACCATCATCCGACTTTAACCATTCAATG  >  USA300TCH1516_ALE/1850833‑1851012

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