Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,513,930 A→T 100% L413F (TTA→TTT)  USA300HOU_RS07375 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,513,9300AT99.7% 1119.5 / ‑4.8 288L413F (TTA→TTT) USA300HOU_RS07375hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base T (151/135);  minor base C (1/0);  total (153/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 289 total aligned reads displayed.

CTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGT  >  USA300TCH1516_ALE/1513843‑1514017
                                                                                       |                                                                                       
ctctTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTAtt                                                                                      >  1:219637/1‑91 (MQ=255)
 tctTGGATAAAGCGAACAAAAAATACTGCAATAAAAGTACCACATCGAGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTg                                                                                     <  2:3811426/91‑1 (MQ=255)
  ctTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGt                                                                                    >  1:3743177/1‑91 (MQ=255)
  ctTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGt                                                                                    <  1:3687005/91‑1 (MQ=255)
  ctTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGt                                                                                    >  1:3425329/1‑91 (MQ=255)
    tGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTAt                                                                                  <  1:6015541/91‑1 (MQ=255)
     ggATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTAtt                                                                                 <  1:4884259/91‑1 (MQ=255)
      gATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTa                                                                                >  2:3949531/1‑91 (MQ=255)
      gATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTa                                                                                <  1:1557405/91‑1 (MQ=255)
        tAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAgg                                                                              <  2:1354972/91‑1 (MQ=255)
         aaaGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGa                                                                             >  2:577324/1‑91 (MQ=255)
         aaaGCATACAAAAAATACTGCAATAAAACTACCACATAGCGTAAAACAACCAGCATCTATATCAATCATAATATGTTATATTGTATTAggc                                                                             <  2:1794938/91‑2 (MQ=255)
          aaGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGat                                                                            >  2:3722605/1‑91 (MQ=255)
           aGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGatt                                                                           >  1:5571994/1‑91 (MQ=255)
            gCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGatta                                                                          >  1:386669/1‑91 (MQ=255)
            gCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGatta                                                                          >  1:3890326/1‑91 (MQ=255)
            gCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGatta                                                                          >  1:94437/1‑91 (MQ=255)
            gCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGatta                                                                          >  2:3340889/1‑91 (MQ=255)
            gCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGatta                                                                          >  2:4176128/1‑91 (MQ=255)
              gTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGattatt                                                                        <  1:2438665/91‑1 (MQ=255)
              gTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGattatt                                                                        >  1:5727487/1‑91 (MQ=255)
               tACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGAttatta                                                                       >  2:1493895/1‑91 (MQ=255)
                aaaaaaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTaa                                                                      <  2:1931974/89‑1 (MQ=255)
                 cAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAAt                                                                     >  2:1441193/1‑91 (MQ=255)
                  aaaaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATg                                                                    <  1:3312810/91‑1 (MQ=255)
                    aaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGcc                                                                  >  2:1286055/1‑91 (MQ=255)
                    aaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGcc                                                                  >  2:3655920/1‑91 (MQ=255)
                     aaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCa                                                                 >  1:201958/1‑91 (MQ=255)
                      aaTACTGCAATAAAAGTACCAAAAAAAATAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCat                                                                <  2:5508765/91‑1 (MQ=255)
                        tACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCatta                                                              >  2:1601653/1‑91 (MQ=255)
                        tACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCatta                                                              <  1:3199454/91‑1 (MQ=255)
                         aCTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                             >  2:5265928/1‑91 (MQ=255)
                         aCTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                             <  2:2636492/91‑1 (MQ=255)
                         aCTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                             <  2:5655307/91‑1 (MQ=255)
                          cTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattatt                                                            <  2:5581577/91‑1 (MQ=255)
                           tGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTAttt                                                           <  1:2580611/91‑1 (MQ=255)
                            gCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTg                                                          <  1:3106703/91‑1 (MQ=255)
                             cAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTgg                                                         >  1:5239643/1‑91 (MQ=255)
                             cAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTgg                                                         >  1:1733827/1‑91 (MQ=255)
                             cAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTgg                                                         >  2:1207786/1‑91 (MQ=255)
                               aTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGAt                                                       <  1:131425/91‑1 (MQ=255)
                               aTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATGATTAATGCCATTATTTGGAt                                                       <  1:5506055/91‑1 (MQ=255)
                                 aaaaGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTa                                                     <  2:3269681/91‑1 (MQ=255)
                                    aGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCa                                                  >  2:1665188/1‑91 (MQ=255)
                                    aGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCa                                                  <  1:5155786/91‑1 (MQ=255)
                                    aGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCa                                                  <  2:2526000/91‑1 (MQ=255)
                                    aGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCa                                                  <  1:3577890/91‑1 (MQ=255)
                                    aGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCa                                                  >  2:5806314/1‑91 (MQ=255)
                                    aGTACCACATCGCGTAAAAAAACCAGCATCTATATCACCCATAATATGTTTTATTGTATTCGGATTATTAATGCCATTATTTGGATTATCa                                                  >  1:4660919/1‑91 (MQ=255)
                                     gTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCAc                                                 <  2:2379646/91‑1 (MQ=255)
                                     gTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTATATTGTATTAGGATTATTAATGCCATTATTTGGATTATCAc                                                 <  1:5974865/91‑1 (MQ=255)
                                      tACCACATCGCGTAAAAAAACCAGCGTCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                                >  2:5032135/1‑91 (MQ=255)
                                      tACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                                >  2:710312/1‑91 (MQ=255)
                                      tACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                                >  1:2732739/1‑91 (MQ=255)
                                      tACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                                >  2:5627175/1‑91 (MQ=255)
                                      tACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                                >  2:4021579/1‑91 (MQ=255)
                                      tACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                                >  2:66790/1‑91 (MQ=255)
                                        ccACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTa                                              >  1:5344749/1‑91 (MQ=255)
                                        ccACATCGCGTAAAAAAACCAGCATCGATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTa                                              <  1:762605/91‑1 (MQ=255)
                                         cacaTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTAt                                             >  1:4736987/1‑91 (MQ=255)
                                           caTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATcc                                           <  1:1627153/91‑1 (MQ=255)
                                            aTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCt                                          >  1:2163121/1‑91 (MQ=255)
                                               gcgTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGt                                       >  1:1825212/1‑91 (MQ=255)
                                                 ttAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAt                                     <  2:3761855/90‑1 (MQ=255)
                                                   aaaaaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAttt                                   <  2:2298650/91‑1 (MQ=255)
                                                        aaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAt                              <  2:5239643/91‑1 (MQ=255)
                                                         aCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAtt                             >  1:1035970/1‑91 (MQ=255)
                                                         aCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAtt                             >  2:615269/1‑91 (MQ=255)
                                                         aCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAtt                             <  2:2989423/91‑1 (MQ=255)
                                                          ccAGCATTTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTg                            <  2:1040252/91‑1 (MQ=255)
                                                          ccAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTg                            >  2:1712653/1‑91 (MQ=255)
                                                           cAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGa                           <  1:4089510/91‑1 (MQ=255)
                                                           cAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGa                           >  2:4200564/1‑91 (MQ=255)
                                                           cAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGa                           <  1:3982935/91‑1 (MQ=255)
                                                            aGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGaa                          >  1:144460/1‑91 (MQ=255)
                                                             gCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAt                         <  1:4158381/91‑1 (MQ=255)
                                                              cATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAtt                        <  2:4987710/91‑1 (MQ=255)
                                                              cATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAtt                        <  2:4287194/91‑1 (MQ=255)
                                                              cATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAtt                        >  1:4957312/1‑91 (MQ=255)
                                                              cATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAtt                        <  2:5974981/91‑1 (MQ=255)
                                                              cATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAtt                        >  1:3795772/1‑91 (MQ=255)
                                                                tCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTaa                      >  1:3947116/1‑91 (MQ=255)
                                                                tCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTaa                      >  1:4238555/1‑91 (MQ=255)
                                                                 cTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAat                     >  2:782482/1‑91 (MQ=255)
                                                                   atatCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAatat                   <  2:3501748/91‑1 (MQ=255)
                                                                   atatCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAatat                   <  2:1960636/91‑1 (MQ=255)
                                                                    tatCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtatt                  >  1:1125967/1‑91 (MQ=255)
                                                                     atCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtatta                 <  1:4398589/91‑1 (MQ=255)
                                                                        aCTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATtatat              >  1:849920/1‑91 (MQ=255)
                                                                        aCTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATtatat              <  1:2357631/91‑1 (MQ=255)
                                                                         cTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTTTCCTTGTATTTATTATTGAATTAATATtatata             >  2:2730132/1‑91 (MQ=255)
                                                                         cTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATtatata             <  2:1841876/91‑1 (MQ=255)
                                                                          tCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTatatat            >  2:5422832/1‑91 (MQ=255)
                                                                           cATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATAtt           >  2:2063222/1‑91 (MQ=255)
                                                                             taataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTaa         <  2:4881702/91‑1 (MQ=255)
                                                                              aataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTaaa        >  2:4511089/1‑91 (MQ=255)
                                                                               atatGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAg       <  1:3530638/91‑1 (MQ=255)
                                                                               atatGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAg       <  2:939889/91‑1 (MQ=255)
                                                                                 atGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAt     >  2:2891310/1‑91 (MQ=255)
                                                                                 atGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAt     >  2:2237953/1‑91 (MQ=255)
                                                                                       |                                                                                       
CTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGT  >  USA300TCH1516_ALE/1513843‑1514017

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