Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,519,951 T→G 8.2% 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,9510TG 8.2% 431.8 / 30.5 142P227P (CCA→CCCnarGRespiratory nitrate reductase 1 alpha chain
Reads supporting (aligned to +/- strand):  ref base T (55/75);  new base G (10/2);  total (65/77)
Fisher's exact test for biased strand distribution p-value = 1.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 157 total aligned reads displayed.

TAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAA  >  USA300TCH1516_ALE/2519861‑2520041
                                                                                          |                                                                                          
  aaGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTgg                                                                                          <  1:1696669/91‑1 (MQ=255)
   aGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCACAAATTTGGGga                                                                                         <  2:1747660/91‑1 (MQ=255)
    ggTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGgag                                                                                        >  1:2468862/1‑91 (MQ=255)
     gTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGgaga                                                                                       <  2:1084312/91‑1 (MQ=255)
      tACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGagag                                                                                      <  2:2799566/91‑1 (MQ=255)
      tACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGagag                                                                                      <  2:2819831/91‑1 (MQ=255)
       aCATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAgg                                                                                     >  1:1787332/1‑91 (MQ=255)
       aCATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAgg                                                                                     >  1:3254060/1‑91 (MQ=255)
        cATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGc                                                                                    >  1:3436278/1‑91 (MQ=255)
        cATTAGAGCCCCACATAATAATGTATGATGCGTGATACCAGTCACGTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGc                                                                                    <  2:3210685/91‑1 (MQ=255)
         aTTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCa                                                                                   >  2:2075236/1‑91 (MQ=255)
         aTTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCa                                                                                   >  2:2075252/1‑91 (MQ=255)
          ttAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAg                                                                                  <  2:661217/91‑1 (MQ=255)
          ttAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAg                                                                                  >  1:890547/1‑91 (MQ=255)
          ttAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAg                                                                                  <  2:661097/91‑1 (MQ=255)
          ttAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAg                                                                                  >  1:890183/1‑91 (MQ=255)
           tAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAgg                                                                                 <  1:2282168/91‑1 (MQ=255)
            agagCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAggt                                                                                <  2:1684095/91‑1 (MQ=255)
              agCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAggtgg                                                                              <  1:709556/91‑1 (MQ=255)
              agCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAggtgg                                                                              <  2:1989534/91‑1 (MQ=255)
                     caTAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCt                                                                       <  2:1910924/91‑1 (MQ=255)
                      ataataATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTg                                                                      >  2:369901/1‑91 (MQ=255)
                       taataaTGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGc                                                                     >  1:3507989/1‑91 (MQ=255)
                         ataatGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCCGGTGGTAACTCTGCAt                                                                   >  2:3778225/1‑91 (MQ=255)
                           aatGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATAc                                                                 >  2:628145/1‑91 (MQ=255)
                           aatGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATAc                                                                 >  1:1664357/1‑91 (MQ=255)
                              gTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAg                                                              <  2:219433/91‑1 (MQ=255)
                                 tgatgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCa                                                           >  2:3340131/1‑91 (MQ=255)
                                  gatgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCAt                                                          >  1:2838275/1‑91 (MQ=255)
                                  gatgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCAt                                                          >  2:1812596/1‑91 (MQ=255)
                                  gatgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCAt                                                          >  1:3718319/1‑91 (MQ=255)
                                  gatgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCAt                                                          <  2:798954/91‑1 (MQ=255)
                                   atgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATa                                                         <  1:2170745/91‑1 (MQ=255)
                                    tgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTTGTAAATCTGCATACCAGTCATaa                                                        >  1:1407250/1‑91 (MQ=255)
                                     gCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATaaa                                                       <  1:1141558/91‑1 (MQ=255)
                                        ttATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAAc                                                    <  2:1760988/91‑1 (MQ=255)
                                         tataCGAGTCACTTGATTCAGGCACATCTGTTTGCTCTCGCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACt                                                   <  1:3672583/91‑1 (MQ=255)
                                         tataCCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCGTAAAAACt                                                   >  1:3457298/1‑91 (MQ=255)
                                          ataCCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTa                                                  <  1:3154406/91‑1 (MQ=255)
                                           taCCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTaa                                                 <  1:111994/91‑1 (MQ=255)
                                           taCCAGTCACGTGAGTCGGGCACATCTGTTGGCGCGCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTaa                                                 <  2:3347311/91‑1 (MQ=255)
                                              cAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCa                                              >  2:3851377/1‑91 (MQ=255)
                                              cAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCa                                              <  1:2966023/91‑1 (MQ=255)
                                               aGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCAt                                             <  1:2508250/91‑1 (MQ=255)
                                               aGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCAt                                             <  1:151441/91‑1 (MQ=255)
                                               aGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCAt                                             <  2:1492012/91‑1 (MQ=255)
                                                 tCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCAttt                                           <  1:1267322/91‑1 (MQ=255)
                                                     ttGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCAGCATACCAGTCATAGAAACTAAGCATTTcacc                                       >  1:3381216/1‑91 (MQ=255)
                                                         ttCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCAccacca                                   >  2:1007350/1‑91 (MQ=255)
                                                          tgaGGCACATCTGTTTGCGCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCaa                                  <  2:3517667/89‑1 (MQ=255)
                                                          tCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCaa                                  <  1:355087/91‑1 (MQ=255)
                                                           cAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAg                                 <  2:3933608/91‑1 (MQ=255)
                                                           cAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAg                                 <  1:2905075/91‑1 (MQ=255)
                                                           cAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAg                                 <  1:1885953/91‑1 (MQ=255)
                                                           cAGGCACATCTGTTTGCTCTCCCCAAATTTGTCGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAACATTTCACCACCAAg                                 <  1:3364849/91‑1 (MQ=255)
                                                            aGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATCAAAACTAAGCATTTCAGCACCAAGc                                >  2:199912/1‑91 (MQ=255)
                                                            aGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGc                                >  2:1877837/1‑91 (MQ=255)
                                                            aGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGc                                <  2:1496544/91‑1 (MQ=255)
                                                             ggCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAgaa                               >  1:2702064/1‑89 (MQ=255)
                                                             ggCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCa                               <  2:3575206/91‑1 (MQ=255)
                                                              gCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCaa                              <  2:3058005/91‑1 (MQ=255)
                                                                agagcTGTTTGCTCTCCCCAGATGGGTGGAGAGCCAGGTGGTAAAACTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAAt                            <  2:1807917/87‑1 (MQ=255)
                                                                acaTCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAAt                            >  1:3316001/1‑91 (MQ=255)
                                                                 caTCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAAtt                           >  1:2532148/1‑91 (MQ=255)
                                                                 caTCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAAtt                           >  2:2479869/1‑91 (MQ=255)
                                                                  aTCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATtg                          <  2:1830606/91‑1 (MQ=255)
                                                                   tCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATtga                         <  2:1240972/91‑1 (MQ=255)
                                                                   tCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATtga                         >  2:2730555/1‑91 (MQ=255)
                                                                     tgattGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCACACCAGTCATAAACCCTAAGCATTTCACCACCAAGCCAATcgatg                       >  2:1240380/4‑91 (MQ=255)
                                                                     tGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATtgatg                       <  1:2333830/91‑1 (MQ=255)
                                                                      gTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATtgatga                      <  1:2668047/91‑1 (MQ=255)
                                                                      gTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATtgatga                      <  2:2270214/91‑1 (MQ=255)
                                                                       tttGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGaa                     >  2:475137/1‑91 (MQ=255)
                                                                        ttGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAAt                    >  2:962117/1‑91 (MQ=255)
                                                                         tGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATc                   <  1:587833/91‑1 (MQ=255)
                                                                          gCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCg                  <  2:3710217/91‑1 (MQ=255)
                                                                          gCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCg                  <  1:1280853/91‑1 (MQ=255)
                                                                          gCTCTCCCCAAATTTCTGGAGAGGCTTGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCg                  <  2:3431200/91‑1 (MQ=255)
                                                                           ctctCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGa                 <  2:3099927/91‑1 (MQ=255)
                                                                           ctctCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGa                 <  2:3441604/91‑1 (MQ=255)
                                                                            tctcCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAg                <  1:2897422/91‑1 (MQ=255)
                                                                             ctcCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGc               <  1:3801146/91‑1 (MQ=255)
                                                                              tcCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCa              >  2:2804732/1‑91 (MQ=255)
                                                                              tcCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCa              >  2:2648949/1‑91 (MQ=255)
                                                                               ccccAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCAc             >  1:161418/1‑91 (MQ=255)
                                                                               ccacAAATTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCAc             >  1:2812133/4‑91 (MQ=255)
                                                                                  cAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctg          >  1:1764384/1‑91 (MQ=255)
                                                                                  cAAATTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctg          >  2:23338/1‑91 (MQ=255)
                                                                                  cAAATTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctg          >  1:1313634/1‑91 (MQ=255)
                                                                                      tttGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCtgctgc      <  2:1678098/91‑1 (MQ=255)
                                                                                      tttGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCtgctgc      >  2:1022967/1‑91 (MQ=255)
                                                                                       ttGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCa     >  1:3431965/1‑91 (MQ=255)
                                                                                       ttGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCa     >  1:135342/1‑91 (MQ=255)
                                                                                        tgtgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCAt    <  1:3110407/91‑1 (MQ=255)
                                                                                         gtgGAGAGGCAGGTGGTAAATCTGCATACCAGTCCTAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATa   >  2:3593517/1‑91 (MQ=255)
                                                                                         gtgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATa   <  1:1307312/91‑1 (MQ=255)
                                                                                         gtgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATa   <  2:1346008/91‑1 (MQ=255)
                                                                                         gtgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATa   >  1:3170401/1‑91 (MQ=255)
                                                                                          tgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATaa  >  2:771568/1‑91 (MQ=255)
                                                                                          tgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATaa  >  1:2382229/1‑91 (MQ=255)
                                                                                          |                                                                                          
TAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAA  >  USA300TCH1516_ALE/2519861‑2520041

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