Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,552,031:1 +G 7.9% coding (1211/1734 nt) USA300HOU_RS13150 ← putative ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,552,0311.G 7.9% 870.9 / 72.2 215coding (1211/1734 nt)USA300HOU_RS13150putative ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (87/111);  new base G (12/5);  total (99/116)
Fisher's exact test for biased strand distribution p-value = 4.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 223 total aligned reads displayed.

CCTGGATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACATAATGTAG  >  USA300TCH1516_ALE/2551941‑2552121
                                                                                           |                                                                                          
ccTGGATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTAc                                                                                             >  1:1161844/1‑91 (MQ=255)
ccTGGATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTAc                                                                                             >  2:2171595/1‑91 (MQ=255)
 cTGGATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACt                                                                                           <  2:35103/91‑1 (MQ=255)
   ggATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTa                                                                                         <  2:1603648/91‑1 (MQ=255)
   ggATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTa                                                                                         >  1:766703/1‑91 (MQ=255)
    gATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTaa                                                                                        <  1:3533164/91‑1 (MQ=255)
      ttGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACAGCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATg                                                                                      >  1:3085712/1‑91 (MQ=255)
      ttGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATg                                                                                      <  1:3725088/91‑1 (MQ=255)
       tGCCAAACAATATGTTATTTTCAATGGTATCATTAAATAGATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGt                                                                                     <  1:3601180/91‑1 (MQ=255)
       tGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGt                                                                                     >  2:1263151/1‑91 (MQ=255)
       tGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGt                                                                                     <  1:1424652/91‑1 (MQ=255)
       tGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAAGTACGTTAATGt                                                                                     <  2:3421111/91‑1 (MQ=255)
        gCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTc                                                                                    >  2:1423214/1‑91 (MQ=255)
         ccAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCa                                                                                   <  2:1900549/91‑1 (MQ=255)
         ccAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCa                                                                                   >  2:13332/1‑91 (MQ=255)
          cAAACAATATGTTATCTTCAATCGTATCATTAAATAAATACACCTTTTGAAAGACTGCACTAATGTTCGCCATTAACGTCCTTAATGTCat                                                                                  <  2:2382120/91‑1 (MQ=255)
                 tatGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTa                                                                           >  2:1661687/1‑91 (MQ=255)
                  atGTTATTTTCAATCGTATCATTCAATAAATACACCTTTTGAAATACTGCACTAATTCTCGACACTAACGTACTTACTGTCATATCTTTac                                                                          >  2:2487583/1‑90 (MQ=255)
                  atGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTaa                                                                          >  2:1164886/1‑91 (MQ=255)
                   tGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAat                                                                         >  2:2113962/1‑91 (MQ=255)
                   tGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAat                                                                         >  2:554175/1‑91 (MQ=255)
                   tGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAat                                                                         >  1:2616100/1‑91 (MQ=255)
                    gTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAata                                                                        >  1:2891342/1‑91 (MQ=255)
                     ttATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAatat                                                                       >  2:3536185/1‑91 (MQ=255)
                       aTTTTCAATCGTACCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCa                                                                     >  1:2242280/1‑91 (MQ=255)
                        ttttCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCaa                                                                    >  2:203545/1‑91 (MQ=255)
                         tttCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAac                                                                   <  1:2845931/91‑1 (MQ=255)
                           tCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAacac                                                                 <  2:264025/91‑1 (MQ=255)
                             aaTCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCa                                                               >  2:824444/1‑91 (MQ=255)
                              aTCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCAt                                                              >  2:1919230/1‑91 (MQ=255)
                               tCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATc                                                             <  1:2091538/91‑1 (MQ=255)
                                cGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATATTTAATATCCACACCATCg                                                            >  2:2522703/1‑91 (MQ=255)
                                 gTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGa                                                           >  2:3672380/1‑91 (MQ=255)
                                   aTCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATg                                                         <  1:3095735/91‑1 (MQ=255)
                                    tCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGc                                                        <  2:3583939/91‑1 (MQ=255)
                                       ttAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGaa                                                     >  2:735954/1‑91 (MQ=255)
                                          aataaataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAatat                                                  <  1:658354/91‑1 (MQ=255)
                                           ataaataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATAtt                                                 <  2:483441/91‑1 (MQ=255)
                                             aaataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTc                                               >  2:1221092/1‑91 (MQ=255)
                                             aaataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTc                                               >  1:1560091/1‑91 (MQ=255)
                                             aaataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTc                                               >  1:1559183/1‑91 (MQ=255)
                                             aaataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTc                                               >  2:3137723/1‑91 (MQ=255)
                                              aataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTcc                                              >  1:1392238/1‑91 (MQ=255)
                                               ataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTTCTTAATGTCATATCTTTAGGATCAGCACCATCGATGCGAATATTTCca                                             <  2:3347679/91‑1 (MQ=255)
                                               ataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCca                                             >  1:253275/1‑91 (MQ=255)
                                               ataCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCca                                             >  1:2495225/1‑91 (MQ=255)
                                                taCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCcat                                            <  1:263109/91‑1 (MQ=255)
                                                taCACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCcat                                            >  1:320061/1‑91 (MQ=255)
                                                  cacCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCcatca                                          >  1:1872527/1‑91 (MQ=255)
                                                   acCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCcatcat                                         >  2:1103966/1‑91 (MQ=255)
                                                   acCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCcatcat                                         >  1:3565669/1‑91 (MQ=255)
                                                      ttttGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGa                                      >  1:108795/1‑91 (MQ=255)
                                                      ttttGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGa                                      >  1:1116435/1‑91 (MQ=255)
                                                       tttGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCGATCATCGat                                     <  1:1950603/91‑1 (MQ=255)
                                                         tGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGatat                                   <  1:2319764/91‑1 (MQ=255)
                                                          gAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATc                                  >  2:2401597/1‑91 (MQ=255)
                                                          gAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATc                                  <  1:3145834/91‑1 (MQ=255)
                                                           aaaTACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCa                                 <  1:1708717/91‑1 (MQ=255)
                                                            aaTACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCAt                                <  1:971511/91‑1 (MQ=255)
                                                             aTACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATa                               <  1:586621/91‑1 (MQ=255)
                                                             aTACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATa                               <  1:1726819/91‑1 (MQ=255)
                                                              tACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATaa                              >  2:1828558/1‑91 (MQ=255)
                                                               aCTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATaaa                             <  1:3017435/91‑1 (MQ=255)
                                                               aCTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATaaa                             <  1:1675551/91‑1 (MQ=255)
                                                                cTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACGATCGATGCGAATATTTCCATCATCGATATCATaaaa                            >  2:3844137/1‑91 (MQ=255)
                                                                  gCACTAATTTTCGACAGTAACGTACTTAATGTCATATCTTTAATATCAACACCAGCGATGCGAATATTTCCATCATCGATATCATAAAAgc                          <  2:1969289/91‑1 (MQ=255)
                                                                   cACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGCTATCATAAAAgcg                         >  2:1728999/1‑91 (MQ=255)
                                                                   cACTAATATTCGACATTAACGTACTTAATGCAATATCTTTAATATCAACACCATCGATACGAATAATTCCATCATCGATATCATAAAAgcg                         >  2:2427277/1‑91 (MQ=255)
                                                                    aCTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATGGATGCGAATATTTCCATCATCGATATCATAAAAgcgc                        <  1:3268636/91‑1 (MQ=255)
                                                                    aCTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAgcgc                        >  1:2471240/1‑91 (MQ=255)
                                                                    aCTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAgcgc                        >  1:3040509/1‑91 (MQ=255)
                                                                    aCTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAgcgc                        >  1:1053122/1‑91 (MQ=255)
                                                                     cTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCa                       <  2:923859/91‑1 (MQ=255)
                                                                      tAATTTTCGACATTAAGGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCaa                      <  2:3706562/91‑1 (MQ=255)
                                                                      tAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCaa                      >  1:3642862/1‑91 (MQ=255)
                                                                      tAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCaa                      >  2:3487922/1‑91 (MQ=255)
                                                                          tttCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGt                  <  2:3040509/91‑1 (MQ=255)
                                                                          tttCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGt                  <  2:3228290/91‑1 (MQ=255)
                                                                           tatgACATTAAGTACGTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATATAAGCGTAAGAGTa                 >  2:807564/4‑91 (MQ=38)
                                                                            tCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTaa                >  2:2350546/1‑91 (MQ=255)
                                                                            tCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTaa                <  1:1204765/91‑1 (MQ=255)
                                                                            atgACATTAAGTACGTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTaa                >  1:1466691/3‑91 (MQ=255)
                                                                             tgACATTAAGTACGTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAg               >  1:2793363/2‑91 (MQ=255)
                                                                              gACATTAAGTACGTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGt              >  1:469861/1‑91 (MQ=255)
                                                                               aCATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTg             <  2:720852/91‑1 (MQ=255)
                                                                               aCATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTg             <  2:721001/91‑1 (MQ=255)
                                                                               aCATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTg             >  2:440512/1‑91 (MQ=255)
                                                                               aCATTAACGTACTGAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTg             <  2:3763784/91‑1 (MQ=255)
                                                                               aCATTAAGTACGTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTg             >  1:1276409/1‑91 (MQ=255)
                                                                               aCATTAAGTACGTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTg             >  2:2564174/1‑91 (MQ=255)
                                                                                cATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGa            <  2:1354425/91‑1 (MQ=255)
                                                                                cATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGa            <  1:3239685/91‑1 (MQ=255)
                                                                                cATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGa            <  2:2891342/91‑1 (MQ=255)
                                                                                 aTTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTACGTGAc           >  2:3440155/1‑91 (MQ=255)
                                                                                   tAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACAt         <  2:2267651/91‑1 (MQ=255)
                                                                                     aCGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACATaa       <  1:3487922/91‑1 (MQ=255)
                                                                                       gTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACATAATg     <  2:578296/91‑1 (MQ=255)
                                                                                       gTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACATAATg     <  2:2200954/91‑1 (MQ=255)
                                                                                        tACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACATAATGt    >  2:3188651/1‑91 (MQ=255)
                                                                                         aCTTAATGTCATATCTGTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACATAATGTa   <  2:1515980/91‑1 (MQ=255)
                                                                                           |                                                                                          
CCTGGATTGCCAAACAATATGTTATTTTCAATCGTATCATTAAATAAATACACCTTTTGAAATACTGCACTAATTTTCGACATTAACGTACTTAATGTCATATCTTTAATATCAACACCATCGATGCGAATATTTCCATCATCGATATCATAAAAGCGCAAGAGTAAGTGACATAATGTAG  >  USA300TCH1516_ALE/2551941‑2552121

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: