Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,985,405 T→A 7.0% P314P (CCA→CCT fumC ← Fumarate hydratase class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,985,4050TA 7.0% 674.1 / 28.5 202P314P (CCA→CCTfumCFumarate hydratase class II
Reads supporting (aligned to +/- strand):  ref base T (98/89);  new base A (4/10);  total (102/100)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

Only 100 of 204 total aligned reads displayed.

AACTGTATCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAATTTTCATTAAGTCT  >  USA300TCH1516_ALE/1985315‑1985495
                                                                                          |                                                                                          
 acTGTATCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTg                                                                                           <  2:2738361/91‑1 (MQ=255)
 acTGTATCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTg                                                                                           <  2:3518063/91‑1 (MQ=255)
  cTGTATCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTgg                                                                                          >  2:1558254/1‑91 (MQ=255)
   tGTATCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGt                                                                                         >  1:2469967/1‑91 (MQ=255)
    ggaTCATGACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGtt                                                                                        <  1:277114/89‑1 (MQ=255)
     tATCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTc                                                                                       <  2:3498693/91‑1 (MQ=255)
     tATCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCAGGTTc                                                                                       <  2:1968702/91‑1 (MQ=255)
     tATCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCAGGTTc                                                                                       <  2:2075839/91‑1 (MQ=255)
      aTCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCa                                                                                      <  2:2935251/91‑1 (MQ=255)
      aTCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCAGGTTCa                                                                                      <  2:1360122/91‑1 (MQ=255)
      aTCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCAGGTTCa                                                                                      <  1:1542818/91‑1 (MQ=255)
       tCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCAt                                                                                     <  2:3028569/91‑1 (MQ=255)
       tCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCAt                                                                                     >  2:2524620/1‑91 (MQ=255)
       tCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCAt                                                                                     <  2:3723397/91‑1 (MQ=255)
       tCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCAGGTTCAt                                                                                     <  2:1436144/91‑1 (MQ=255)
        cATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCAGGTTCAtc                                                                                    <  1:2585399/91‑2 (MQ=255)
            aCCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCa                                                                                >  2:1175428/1‑91 (MQ=255)
                aTTACTTGGACTGCAACCATTGTTAGCATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGa                                                                            <  1:3612555/91‑1 (MQ=255)
                aTTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGa                                                                            >  1:337825/1‑91 (MQ=255)
                  tACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATa                                                                          <  1:912826/91‑1 (MQ=255)
                  tACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATa                                                                          >  2:1071263/1‑91 (MQ=255)
                   aCTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAg                                                                         >  1:231146/1‑91 (MQ=255)
                    cTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGa                                                                        <  2:2710179/91‑1 (MQ=255)
                       ggACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAAt                                                                     >  2:3574411/1‑91 (MQ=255)
                       ggACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAAt                                                                     >  2:2284153/1‑91 (MQ=255)
                       ggACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAAt                                                                     >  1:1643817/1‑91 (MQ=255)
                       ggACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAAt                                                                     <  1:1063086/91‑1 (MQ=255)
                        gACTGCAACCATCGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAGAtt                                                                    >  1:3005084/1‑91 (MQ=255)
                         aCTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAAttt                                                                   >  2:1974398/1‑91 (MQ=255)
                             cAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGc                                                               >  1:2527999/1‑91 (MQ=255)
                             cAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGc                                                               >  1:1317894/1‑91 (MQ=255)
                             cAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGc                                                               >  1:2121166/1‑91 (MQ=255)
                              aaCCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGcc                                                              <  1:2433987/91‑1 (MQ=255)
                              aaCCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGcc                                                              >  1:1791784/1‑91 (MQ=255)
                                ccATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCaa                                                            >  2:3203546/1‑91 (MQ=255)
                                ccATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCaa                                                            >  1:2085911/1‑91 (MQ=255)
                                   ttGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGCCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAACTTCTGCCAAAcc                                                         >  2:1053303/1‑91 (MQ=255)
                                     gTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAg                                                       <  2:2505432/91‑1 (MQ=255)
                                       tAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCt                                                     <  2:2110652/91‑1 (MQ=255)
                                       tAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCt                                                     >  1:1417706/1‑91 (MQ=255)
                                       tAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCt                                                     >  2:2156189/1‑91 (MQ=255)
                                       tAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCt                                                     <  1:2262922/91‑1 (MQ=255)
                                         acacttCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAAc                                                          >  2:1722740/5‑84 (MQ=255)
                                         acacTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAAc                                                          <  1:1722740/84‑1 (MQ=255)
                                         aCATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCg                                                   >  1:3526295/1‑91 (MQ=255)
                                          cATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGt                                                  >  1:2102680/1‑91 (MQ=255)
                                          cATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGt                                                  <  2:3777168/91‑1 (MQ=255)
                                           aTTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTg                                                 <  2:3875889/91‑1 (MQ=255)
                                             ttCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGc                                               >  1:134870/1‑91 (MQ=255)
                                              tCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGCTAGAAATTTCTGCCAAACCAGCTCGTGGcc                                              >  2:2374679/1‑91 (MQ=255)
                                              tCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGcc                                              <  2:337825/91‑1 (MQ=255)
                                               cacaTTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGccc                                             >  1:1930801/1‑91 (MQ=255)
                                               cacaTTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGccc                                             >  2:744000/1‑91 (MQ=255)
                                               cacaTTGTGTAGGATTAACTTTACCAGGCATAATTGAGGAACCGGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGccc                                             <  2:3828443/91‑1 (MQ=255)
                                                 caTTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTg                                           >  2:115504/1‑91 (MQ=255)
                                                 caTTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTg                                           >  2:3356896/1‑91 (MQ=255)
                                                  aTTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGa                                          <  1:1229236/91‑1 (MQ=255)
                                                    tgtgTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAg                                        >  1:3049772/1‑91 (MQ=255)
                                                     gtgtAGGTTTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGACCCTGAAGc                                       <  1:1951174/91‑1 (MQ=255)
                                                      tgtAGGATTAGCTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGcc                                      <  2:515881/91‑1 (MQ=255)
                                                      tgtAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGcc                                      <  2:1485338/91‑1 (MQ=255)
                                                      tgtAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGcc                                      >  1:719122/1‑91 (MQ=255)
                                                      tgtAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAACTCGTGGCCCTGAAGcc                                      <  2:2085911/91‑1 (MQ=255)
                                                       gtAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCa                                     >  2:2423524/1‑91 (MQ=255)
                                                        tAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCaa                                    <  1:3824617/91‑1 (MQ=255)
                                                         aGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGTTCGTGGCCCTGAAGCCAAc                                   >  1:2459341/1‑91 (MQ=255)
                                                         aGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAAc                                   >  2:1594145/1‑91 (MQ=255)
                                                         aGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAAc                                   >  2:2692603/1‑91 (MQ=255)
                                                          ggATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAAcc                                  >  2:1445696/1‑91 (MQ=255)
                                                          ggATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAAcc                                  <  1:2743232/91‑1 (MQ=255)
                                                          ggATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAAcc                                  <  2:1329379/91‑1 (MQ=255)
                                                             ttAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCAtc                               <  2:2345334/91‑1 (MQ=255)
                                                               aaCTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCAtctc                             <  1:3708255/91‑1 (MQ=255)
                                                                aCTTTACCAGGCATTATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTca                            <  1:2207025/91‑1 (MQ=255)
                                                                aCTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTca                            >  1:1098137/1‑91 (MQ=255)
                                                                 cTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTcac                           >  2:1396950/1‑91 (MQ=255)
                                                                    tACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCAcatc                        >  1:84011/1‑91 (MQ=255)
                                                                    tACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCAcatc                        >  2:573439/1‑91 (MQ=255)
                                                                     aCCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCAcatca                       <  2:376961/91‑1 (MQ=255)
                                                                      ccAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCAcatcat                      >  1:1905187/1‑91 (MQ=255)
                                                                        cggCATAATTGATGAACCTGGTTCATTTTCAGGGAAAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTa                    <  1:612054/90‑1 (MQ=255)
                                                                        aGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTa                    >  1:3547206/1‑91 (MQ=255)
                                                                        aGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTa                    <  1:1411975/91‑1 (MQ=255)
                                                                         ggCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAg                   >  2:2804183/1‑91 (MQ=255)
                                                                             tAATTGATGAACCTGGTTCATTTTCAGGGGTAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAAt               <  1:3115985/91‑1 (MQ=255)
                                                                             tAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAAt               <  2:2668203/91‑1 (MQ=255)
                                                                                ttGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAAtttt            >  2:613771/1‑91 (MQ=255)
                                                                                ttGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAAtttt            <  2:397321/91‑1 (MQ=255)
                                                                                 tgatgaACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAATTTTc           >  1:2836723/1‑91 (MQ=255)
                                                                                  gatgaACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAATTTTCa          >  2:1915632/1‑91 (MQ=255)
                                                                                  gatgaACCTGGTTCATGTTCAGGGAGAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAATTTTCa          <  1:1855138/91‑1 (MQ=255)
                                                                                   atgaACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGc                  >  2:1853069/1‑82 (MQ=255)
                                                                                     gaACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAATTTTCATTa       >  1:2084068/1‑91 (MQ=255)
                                                                                     gaACCAGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAATTTTCATTa       >  1:470895/1‑91 (MQ=255)
                                                                                      aaCCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAATTTTCATTaa      <  2:123535/91‑1 (MQ=255)
                                                                                       aCCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAATTTTCATTAAg     <  2:3005084/91‑1 (MQ=255)
                                                                                         cTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAATTTTCATTAAGtc   <  1:1737061/91‑1 (MQ=255)
                                                                                         cTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAATTTTCATTAAGtc   >  1:1962723/1‑91 (MQ=255)
                                                                                         cTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAATTTTCATTAAGtc   <  1:774945/91‑1 (MQ=255)
                                                                                          tGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAATTTTCATTAAGtct  <  1:1366381/91‑1 (MQ=255)
                                                                                          |                                                                                          
AACTGTATCATTACCCATTACTTGGACTGCAACCATTGTTAACATTTCACATTGTGTAGGATTAACTTTACCAGGCATAATTGATGAACCTGGTTCATTTTCAGGGATAGAAATTTCTGCCAAACCAGCTCGTGGCCCTGAAGCCAACCATCTCACATCATTAGCAATTTTCATTAAGTCT  >  USA300TCH1516_ALE/1985315‑1985495

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: