Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,744,454 C→T G73D (GGC→GAC)  apt ← Adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,744,4540CT97.6% 567.3 / ‑0.5 167G73D (GGC→GAC) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (85/78);  total (88/79)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

GTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTT  >  USA300TCH1516_ALE/1744367‑1744540
                                                                                       |                                                                                      
gTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCaa                                                                                     <  1:2662420/91‑1 (MQ=255)
   aaaaCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATcc                                                                                  >  2:883756/1‑91 (MQ=255)
       cATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCAt                                                                              <  1:318484/91‑1 (MQ=255)
         tttGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTg                                                                            <  2:538134/91‑1 (MQ=255)
         tttGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTg                                                                            >  1:2765302/1‑91 (MQ=255)
         tttGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTg                                                                            <  1:930605/91‑1 (MQ=255)
          ttGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGa                                                                           >  2:1891698/1‑91 (MQ=255)
          ttGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGa                                                                           >  2:745898/1‑91 (MQ=255)
           tGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGaa                                                                          <  1:1169790/91‑1 (MQ=255)
           tGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGaa                                                                          <  1:1170953/91‑1 (MQ=255)
            gTACCATATTCTAGGTCATACTCATAACGAATGACTTCGCGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAAt                                                                         >  1:1016273/1‑91 (MQ=255)
            gTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAAt                                                                         <  1:910758/91‑1 (MQ=255)
            gTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAAt                                                                         >  2:2327483/1‑91 (MQ=255)
              aCCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATaa                                                                       >  2:3590425/1‑91 (MQ=255)
              aCCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATaa                                                                       >  2:3592505/1‑91 (MQ=255)
               ccATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAg                                                                      >  2:1446209/1‑91 (MQ=255)
                  tatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTa                                                                   <  2:3546599/91‑1 (MQ=255)
                  tatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTa                                                                   >  1:2476488/1‑91 (MQ=255)
                  tatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTa                                                                   <  1:442239/91‑1 (MQ=255)
                  tatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTa                                                                   <  1:1882710/91‑1 (MQ=255)
                  tatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTa                                                                   >  1:1405212/1‑91 (MQ=255)
                  tatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTa                                                                   <  2:2909728/91‑1 (MQ=255)
                   atTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTAc                                                                  <  2:1066887/91‑1 (MQ=255)
                    ttCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACa                                                                 <  2:2876965/91‑1 (MQ=255)
                     tCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAg                                                                <  1:2138145/91‑1 (MQ=255)
                     tCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAg                                                                >  1:315201/1‑91 (MQ=255)
                     tCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAg                                                                >  1:1588100/1‑91 (MQ=255)
                      cTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTATCAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAgg                                                               <  2:2157640/91‑1 (MQ=255)
                      cTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAgg                                                               <  1:1547734/91‑1 (MQ=255)
                      cTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAgg                                                               >  1:341195/1‑91 (MQ=255)
                        aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGAc                                                             >  1:3299675/1‑91 (MQ=255)
                        aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGAc                                                             <  1:2238165/91‑1 (MQ=255)
                           tCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGc                                                          >  1:2531894/1‑91 (MQ=255)
                           tCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGc                                                          <  2:2181732/91‑1 (MQ=255)
                            cATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGcc                                                         >  2:1848595/1‑91 (MQ=255)
                            cATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGcc                                                         <  2:3103084/91‑1 (MQ=255)
                               aCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAAt                                                      <  2:2163620/91‑1 (MQ=255)
                               aCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAAt                                                      <  2:177752/91‑1 (MQ=255)
                               aCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGGGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAAt                                                      <  2:2163377/91‑1 (MQ=255)
                                cTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATg                                                     >  2:2557373/1‑91 (MQ=255)
                                  cATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGAt                                                   >  1:728918/1‑91 (MQ=255)
                                     aaCGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaaa                                                <  1:888586/91‑1 (MQ=255)
                                        gAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGcc                                             >  1:1033409/1‑91 (MQ=255)
                                        gAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGcc                                             <  2:2629699/91‑1 (MQ=255)
                                        gAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGcc                                             <  1:2053098/91‑1 (MQ=255)
                                          aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAc                                           >  2:75158/1‑91 (MQ=255)
                                          aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAc                                           >  2:586297/1‑91 (MQ=255)
                                          aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAc                                           >  2:608553/1‑91 (MQ=255)
                                          aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAc                                           >  1:3392824/1‑91 (MQ=255)
                                           tGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcg                                          <  1:841241/91‑1 (MQ=255)
                                            gACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATACGCTACAGGACAGCCAATGATAAAGCCAcgc                                         >  2:2715446/1‑91 (MQ=255)
                                            gACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcgc                                         >  2:2715418/1‑91 (MQ=255)
                                            gACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcgc                                         >  2:2306600/1‑91 (MQ=255)
                                              cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACgcgc                                       >  1:2480458/1‑91 (MQ=255)
                                              cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACgcgc                                       >  1:1293122/1‑91 (MQ=255)
                                              cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACgcgc                                       <  2:1320593/91‑1 (MQ=255)
                                              cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGAAAAGCCACGCGCt                                      >  1:3577051/1‑91 (MQ=255)
                                                tCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCtt                                     <  1:2033864/91‑1 (MQ=255)
                                                 cACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTc                                    >  2:3502814/1‑91 (MQ=255)
                                                 cACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTc                                    >  2:1130123/1‑91 (MQ=255)
                                                   cGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAg                                  >  1:3052803/1‑91 (MQ=255)
                                                     aGGTAATTTCCGTTCTTTTCTAAGAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGt                                <  2:2774411/91‑1 (MQ=255)
                                                        tAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCt                             <  1:2357145/91‑1 (MQ=255)
                                                          aTTTCCCTTTTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTAc                           <  1:3468459/91‑1 (MQ=255)
                                                          aTTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTAc                           <  1:2872692/91‑1 (MQ=255)
                                                             tCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAAc                        <  1:1700283/91‑1 (MQ=255)
                                                             tCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAAc                        >  1:110511/1‑91 (MQ=255)
                                                             tCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAAc                        >  2:1259772/1‑91 (MQ=255)
                                                              cccTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACg                       >  2:877168/1‑91 (MQ=255)
                                                                 ttcttTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGata                    <  2:3211282/91‑1 (MQ=255)
                                                                   cttTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                  >  2:2619644/1‑91 (MQ=255)
                                                                   cttTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                  <  1:2644899/91‑1 (MQ=255)
                                                                   cttTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                  <  1:2702435/91‑1 (MQ=255)
                                                                   cttTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGAGAAAGCCACGCGCTTCAGGTCCTACAACGATATc                  >  1:2978307/1‑91 (MQ=255)
                                                                   cttTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGAGAAAGCCACGCGCTTCAGGTCCTACAACGATATc                  >  1:2978709/1‑91 (MQ=255)
                                                                    ttttCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCa                 <  2:2011171/91‑1 (MQ=255)
                                                                    ttttCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCa                 >  2:1390594/1‑91 (MQ=255)
                                                                    ttttCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCa                 <  2:343562/91‑1 (MQ=255)
                                                                    ttttCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCa                 >  2:1390367/1‑91 (MQ=255)
                                                                      ttCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAAc               <  2:2997296/91‑1 (MQ=255)
                                                                       tCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACa              >  2:2346939/1‑91 (MQ=255)
                                                                        cTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAACGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCaaaaa             >  2:2097971/1‑88 (MQ=255)
                                                                        cTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAt             >  2:505775/1‑91 (MQ=255)
                                                                         tAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtc            >  1:1417340/1‑91 (MQ=255)
                                                                         tAACAGGTGCAAAGTCAATCCCCATAGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtc            >  1:224646/1‑91 (MQ=255)
                                                                          aaCAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct           >  2:781516/1‑91 (MQ=255)
                                                                          aaCAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct           <  2:93605/91‑1 (MQ=255)
                                                                          aaCAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct           >  2:183174/1‑91 (MQ=255)
                                                                          aaCAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtcc           >  2:928002/1‑90 (MQ=255)
                                                                            cAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct         >  2:3395104/1‑91 (MQ=255)
                                                                            cAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct         >  2:196351/1‑91 (MQ=255)
                                                                            cAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct         <  1:1438699/91‑1 (MQ=255)
                                                                               gTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTc      >  2:1153895/1‑91 (MQ=255)
                                                                               gTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTc      >  1:973485/1‑91 (MQ=255)
                                                                               gTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTc      >  1:1514276/1‑91 (MQ=255)
                                                                                 gCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCATGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtt    >  2:2230821/1‑91 (MQ=255)
                                                                                 gCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtt    >  1:2768940/1‑91 (MQ=255)
                                                                                 gCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtt    >  2:1593265/1‑91 (MQ=255)
                                                                                  cAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCttt   <  1:51173/91‑1 (MQ=255)
                                                                                   aaaGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtttt  >  1:2441891/1‑91 (MQ=255)
                                                                                       |                                                                                      
GTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTT  >  USA300TCH1516_ALE/1744367‑1744540

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: