Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,519,951 T→G 13.5% P227P (CCA→CCC narG ← Respiratory nitrate reductase 1 alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,519,9510TG13.5% 469.9 / 70.4 169P227P (CCA→CCCnarGRespiratory nitrate reductase 1 alpha chain
Reads supporting (aligned to +/- strand):  ref base T (85/60);  new base G (21/2);  total (107/62)
Fisher's exact test for biased strand distribution p-value = 2.08e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 178 total aligned reads displayed.

AAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAAC  >  USA300TCH1516_ALE/2519863‑2520042
                                                                                        |                                                                                           
aaGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTgg                                                                                           <  2:2302338/91‑1 (MQ=255)
 aGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATTCCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGga                                                                                          <  1:1059443/91‑1 (MQ=255)
   gTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGgaga                                                                                        <  2:1841858/91‑1 (MQ=255)
   gTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGgaga                                                                                        <  2:1323455/91‑1 (MQ=255)
       aTTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCa                                                                                    <  2:2944677/91‑1 (MQ=255)
        ttAGAGCCCCACATACTAATGTATGATGCGTTATACCAGTCACTTGCTTCAGGCACATCTGATTGCTCTCCCCATGTTTGTGGAGTGGCAg                                                                                   >  2:2590884/1‑91 (MQ=255)
        ttAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCACAAATTTGGGGagagcaaa                                                                                   <  1:2054580/91‑5 (MQ=39)
          agagCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAggt                                                                                 >  1:2145024/1‑91 (MQ=255)
          agagCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAggt                                                                                 >  1:2385874/1‑91 (MQ=255)
          agagCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAggt                                                                                 >  1:1397973/1‑91 (MQ=255)
          agagCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAggt                                                                                 <  1:1625133/91‑1 (MQ=255)
            agCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCCGTTTGCTCTCCCCAAATTTGTGGAGAGGCAggtgg                                                                               >  1:964646/1‑91 (MQ=255)
            agCCCCACATAATAAAGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCCGTTTGCTCTCCCCAAATTTGTGGAGAGGCAggtgg                                                                               >  1:964645/1‑91 (MQ=255)
             gCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCACAAATTTGTGGAGAGGCAggtggt                                                                              <  2:2434992/91‑1 (MQ=255)
                ccACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTaaa                                                                           >  1:251783/1‑91 (MQ=255)
                ccACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTaaa                                                                           <  1:707316/91‑1 (MQ=255)
                ccACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTaaa                                                                           >  1:52833/1‑91 (MQ=255)
                ccACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTaaa                                                                           >  1:645493/1‑91 (MQ=255)
                ccACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTaaa                                                                           >  2:2087585/1‑91 (MQ=255)
                 cacaTAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAAt                                                                          >  2:1430687/1‑91 (MQ=255)
                    ataataATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTg                                                                       >  2:1593671/1‑91 (MQ=255)
                    ataataATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTg                                                                       >  1:1104684/1‑91 (MQ=255)
                         aatGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATAc                                                                  >  1:1728248/1‑91 (MQ=255)
                          atGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGAATAcc                                                                 >  2:3137575/1‑91 (MQ=255)
                          atGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTCTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATAcc                                                                 <  1:3087662/91‑1 (MQ=255)
                            gTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAg                                                               >  1:1819042/1‑91 (MQ=255)
                            gTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAg                                                               <  2:1419970/91‑1 (MQ=255)
                             tATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGt                                                              >  1:296364/1‑91 (MQ=255)
                             tATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGt                                                              >  2:617165/1‑91 (MQ=255)
                             tATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGt                                                              >  1:2600429/1‑91 (MQ=255)
                             tATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGt                                                              <  1:1719842/91‑1 (MQ=255)
                              atgatgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTc                                                             >  1:2095741/1‑91 (MQ=255)
                               tgatgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCa                                                            >  1:1975469/1‑91 (MQ=255)
                                gatgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCCTACCAGTCAt                                                           >  2:3348716/1‑91 (MQ=255)
                                gatgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCAt                                                           <  1:2662341/91‑1 (MQ=255)
                                gatgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCAt                                                           >  2:1623976/1‑91 (MQ=255)
                                gatgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCAt                                                           <  2:1790637/91‑1 (MQ=255)
                                 atgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATa                                                          >  2:1090486/1‑91 (MQ=255)
                                 atgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATa                                                          <  1:402617/91‑1 (MQ=255)
                                   gCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATaaa                                                        <  2:2551500/91‑1 (MQ=255)
                                     gTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATaaaaa                                                      >  1:1157423/1‑91 (MQ=255)
                                      ttATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAAc                                                     >  2:1537587/1‑91 (MQ=255)
                                      ttATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAAc                                                     >  1:368245/1‑91 (MQ=255)
                                          aCCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAg                                                 >  2:2143189/1‑91 (MQ=255)
                                            cAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAACAACTAAGCa                                               >  1:3711459/1‑91 (MQ=255)
                                            cAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCa                                               >  1:2044673/1‑91 (MQ=255)
                                            cAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCa                                               <  1:1651942/91‑1 (MQ=255)
                                              gTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCAtt                                             >  2:3493722/1‑91 (MQ=255)
                                              gTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCAtt                                             >  1:235866/1‑91 (MQ=255)
                                              gTCACTTAATTCAGGCACATCTGTTTGCTCTCCCCAGATTAGTGGAGAGGCAGGTGGTAAATCTGCCAACCAGTCAGAAAAACTAAGCAtt                                             >  2:3612593/1‑91 (MQ=255)
                                                 aCTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTca                                          <  2:1881596/91‑1 (MQ=255)
                                                   ttGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTcacc                                        >  2:2410492/1‑91 (MQ=255)
                                                    tGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTcacca                                       <  2:3665509/91‑1 (MQ=255)
                                                     gATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTcaccac                                      >  2:3760989/1‑91 (MQ=255)
                                                      aTTCAGGCACATCTGTTTGCTCTCCCCAAATTTATGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCaccacc                                     <  1:3183963/91‑1 (MQ=255)
                                                        tCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCaa                                   <  1:2270936/91‑1 (MQ=255)
                                                           ggCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCa                                >  1:1652258/1‑91 (MQ=255)
                                                            gCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCaa                               >  1:764336/1‑91 (MQ=255)
                                                            gCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCaa                               >  2:2472134/1‑91 (MQ=255)
                                                              acaTCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAAt                             <  2:2880691/91‑1 (MQ=255)
                                                              acaTCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTACACCACCAAGCAAAt                             <  2:2385874/91‑1 (MQ=255)
                                                               caTCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCCGTCATAAAAACTAAGCATTTCACCACCAAGCAAAtt                            >  2:3055054/1‑91 (MQ=255)
                                                               caTCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAAtt                            <  1:333247/91‑1 (MQ=255)
                                                                     tttGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGaa                      >  1:521936/1‑91 (MQ=255)
                                                                     tttGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGaa                      >  1:3863462/1‑91 (MQ=255)
                                                                     tttGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGaa                      >  1:1291326/1‑91 (MQ=255)
                                                                     tttGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGaa                      >  2:1476528/1‑91 (MQ=255)
                                                                       tGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATc                    <  1:3063598/91‑1 (MQ=255)
                                                                        gCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCg                   <  1:1035323/91‑1 (MQ=255)
                                                                         ctctCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCTGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGATTCGa                  <  1:3069741/91‑1 (MQ=255)
                                                                          tctcCTCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAg                 <  2:557429/91‑1 (MQ=255)
                                                                            tccacAAATTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCa               >  1:3777659/5‑91 (MQ=255)
                                                                             ccccAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCAc              <  2:3442533/91‑1 (MQ=255)
                                                                              cccAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCAcc             >  1:3691947/1‑91 (MQ=255)
                                                                              cacAAATTTGGGGAGAGGCAGGTGGTAAATCTGCCTACCAGTCATAAAAACTAAGTATTTCACCACCAAGCAAATTGATGAATCGAGCAcc             >  2:72204/3‑91 (MQ=255)
                                                                              cacAAATTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCAcc             >  2:2937908/3‑91 (MQ=255)
                                                                               ccAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACct            <  1:1016739/91‑1 (MQ=255)
                                                                               ccAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACct            <  1:1016826/91‑1 (MQ=255)
                                                                               acAAATTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACct            >  1:1293963/2‑91 (MQ=255)
                                                                               acAAATTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACct            >  1:1313215/2‑91 (MQ=255)
                                                                                 aaaTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgc          >  1:2295743/1‑91 (MQ=255)
                                                                                 aaaTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgc          >  1:2157901/1‑91 (MQ=255)
                                                                                 aaaTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgc          >  1:95905/1‑91 (MQ=255)
                                                                                 aaaTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgc          >  1:856953/1‑91 (MQ=255)
                                                                                  aaTTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgct         >  2:2496779/1‑91 (MQ=255)
                                                                                  aaTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgct         >  1:2928830/1‑91 (MQ=255)
                                                                                  aaTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgct         >  1:2255895/1‑91 (MQ=255)
                                                                                  aaTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgct         >  1:1703763/1‑91 (MQ=255)
                                                                                   aTTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgctg        >  1:241518/1‑91 (MQ=255)
                                                                                   aTTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgctg        <  1:711627/91‑1 (MQ=255)
                                                                                   aTTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgctg        <  1:1233939/91‑1 (MQ=255)
                                                                                   aTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgctg        >  2:2784787/1‑91 (MQ=255)
                                                                                   aTTTGGGGAGAGACAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgctg        >  1:2973516/1‑91 (MQ=255)
                                                                                    tttGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCtgctgc       <  2:2232571/91‑1 (MQ=255)
                                                                                    tttGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCtgctgc       >  1:694529/1‑91 (MQ=255)
                                                                                    tttGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCtgctgc       >  2:1039435/1‑91 (MQ=255)
                                                                                     ttGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCa      >  1:111452/1‑91 (MQ=255)
                                                                                       gtgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATa    >  1:1660109/1‑91 (MQ=255)
                                                                                        tgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATaa   >  1:1667400/1‑91 (MQ=255)
                                                                                        tgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAAc  >  2:3760584/1‑91 (MQ=255)
                                                                                        |                                                                                           
AAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAAC  >  USA300TCH1516_ALE/2519863‑2520042

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: