Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,653,690 G→T 10.7% Q97K (CAA→AAA)  aroK ← Shikimate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,653,6900GT10.7% 207.0 / 8.6 74Q97K (CAA→AAA) aroKShikimate kinase
Reads supporting (aligned to +/- strand):  ref base G (38/28);  new base T (6/2);  total (44/30)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

ATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAATATTGTTCACCATGTTTACAAAATA  >  USA300TCH1516_ALE/1653550‑1653827
                                                                                                                                            |                                                                                                                                         
aTTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTg                                                                                                                                           <  2:239578/141‑1 (MQ=255)
aTTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTg                                                                                                                                           <  2:239586/141‑1 (MQ=255)
  tataTCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGGGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGAt                                                                                                                                         >  1:1089398/1‑141 (MQ=255)
             ttCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTAATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATCTATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCaaaaaa                                                                                                                              <  2:115641/141‑1 (MQ=255)
                 gcgAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTaa                                                                                                                          <  1:89467/141‑1 (MQ=255)
                        aTACAATTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCctc                                                                                                                   <  2:985726/141‑1 (MQ=255)
                        aTACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCctc                                                                                                                   <  2:712690/141‑1 (MQ=255)
                        aTACAAGTCATTTAACTGCTTGATTGTCTTATCATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAATCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCctc                                                                                                                   <  1:541944/141‑1 (MQ=255)
                            aaGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCa                                                                                                               <  2:395233/141‑1 (MQ=255)
                               tCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCACATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCActc                                                                                                            <  1:1051001/141‑1 (MQ=255)
                                 aTTTAACTGCTTGATTGTCTTATTTTTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACtctc                                                                                                          <  2:258330/141‑1 (MQ=255)
                                    tAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCaat                                                                                                       >  1:16997/1‑141 (MQ=255)
                                      aCTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCaataa                                                                                                     >  2:1080556/1‑141 (MQ=255)
                                          cTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATacc                                                                                                 >  2:978826/1‑141 (MQ=255)
                                             taTTGTCTTATTATTTGCATTAGGTCGATGTGGTTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATagcacc                                                                                              <  1:751530/140‑1 (MQ=255)
                                             gATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGCTTTTTCCAAACATTAAATGCCTCTTCATTCTCAATAATaccacc                                                                                              <  2:645066/141‑1 (MQ=255)
                                             gATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATaccacc                                                                                              <  1:902530/141‑1 (MQ=255)
                                             gATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTGTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATaccacc                                                                                              <  1:758531/141‑1 (MQ=255)
                                             gATTGTCTGATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATaccacc                                                                                              <  2:179509/141‑1 (MQ=255)
                                              gttGTCTTTTTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATaccacca                                                                                             <  1:371839/140‑1 (MQ=255)
                                              aTTGTCTTATTATTTGCATTAGGGCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATCATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATaccacca                                                                                             <  2:8930/141‑1 (MQ=255)
                                               ttGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATACTaccacccc                                                                                            >  2:339978/1‑139 (MQ=255)
                                               ttGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATaccaccac                                                                                            >  1:950362/1‑141 (MQ=255)
                                               ttGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATaccaccac                                                                                            >  1:627663/1‑141 (MQ=255)
                                                tGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAGTGTTTTTTTGATTTTTCAAAAAATTAAATTCCTCTTCACTCTCAATAATaccaccacc                                                                                           >  2:728684/1‑141 (MQ=255)
                                                 ttCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATAccaccacca                                                                                          <  1:714660/140‑1 (MQ=255)
                                                 gTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATAccaccacca                                                                                          <  1:4929/141‑1 (MQ=255)
                                                 gTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATAccaccacca                                                                                          <  1:1080190/141‑1 (MQ=255)
                                                 gTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTGTCATATAATTAAATGCCTCTTCACTATCAATAATAcccccacca                                                                                          >  2:403105/1‑141 (MQ=255)
                                                 gTCTTATTATTTGCAGTAGTTGCATGTGGGTCATCATTCCTTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATAccaccacca                                                                                          <  2:114224/141‑1 (MQ=255)
                                                               tttAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATc                                                                            <  2:672955/140‑1 (MQ=255)
                                                                 tAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCtc                                                                          >  1:55463/1‑140 (MQ=255)
                                                                   ggTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCGTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCtccc                                                                        >  2:359490/1‑138 (MQ=255)
                                                                     tCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCACATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAgg                                                                      >  1:317540/1‑140 (MQ=255)
                                                                        aTGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTTTTTTTTTTTTTAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTGGCAATTATATCTGccttttt                                                                   >  2:272104/1‑135 (MQ=255)
                                                                        aTGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGtt                                                                   >  1:1013810/1‑141 (MQ=255)
                                                                           tgGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAAt                                                                >  2:985748/1‑141 (MQ=255)
                                                                            gggTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCCCTAGTAGAATTTATATCTGTATTTTTAATa                                                               >  2:213049/1‑141 (MQ=255)
                                                                            gggTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATa                                                               >  2:580025/1‑141 (MQ=255)
                                                                            gggTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATa                                                               >  1:362726/1‑141 (MQ=255)
                                                                                  tcatTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCt                                                         >  1:942870/1‑141 (MQ=255)
                                                                                  tcatTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCt                                                         >  1:308182/1‑141 (MQ=255)
                                                                                           tgACTATATTGATCATCAAGCTTCCAATCTAACCAAATAATGTTTTTTTGTTTTTTCAAAAAATTAAATGCCTCTACACTCTCAATAGTACCCCCCCCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACAt                                                <  2:281148/140‑1 (MQ=255)
                                                                                             aCTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTTTTTTTTTTTAAAATTAAATACCTCTTCACTCTCAATAATACCAGCACCAGTAGCAATTATTTATGCATTATTATTTCATTCTTTCACAAAtgt                                              >  2:878168/1‑141 (MQ=255)
                                                                                             aCTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTTTTTTTTTTAAAAATTAAATGCCTCTTCACTCTCAATAATACTACCACCAGTTGCAATTATTTCTGCATTCTTAATACATTCTTGCAAACAttt                                              >  2:84528/1‑139 (MQ=255)
                                                                                                    ataataTCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAActct                                       >  1:28164/1‑141 (MQ=255)
                                                                                                     taataTCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTa                                      <  1:396456/141‑1 (MQ=255)
                                                                                                     taataTCAATATTACAATCTAACCAAATAATGTTTTTTTGATTCTTCAAAAAATTAAATGCCTTTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTa                                      <  2:301170/141‑1 (MQ=255)
                                                                                                        tatCAATATTACAATCTAACCAAATAATGTTTTTTTTATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAAt                                   <  1:594730/141‑1 (MQ=255)
                                                                                                        tatCAATATTACAATCTAACCAAATAATGTTTTTTTTATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAAt                                   <  1:380449/141‑1 (MQ=255)
                                                                                                         atCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAAtt                                  <  1:1095441/141‑1 (MQ=255)
                                                                                                           cAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAGCATGTGAACTCTAAATTcc                                <  1:662799/141‑1 (MQ=255)
                                                                                                             atatTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTg                              <  2:16997/141‑1 (MQ=255)
                                                                                                              tatTACAATCTAACCAAATAATGTTTTTTTTTTTTTTTTTAAAATAAAATGCATCATCAGACTCAATAATACCTCCATCAGTATCAATTATATCTGCATTGTTAATACAATCTTCCAAATATGAGATATCTAAATACCTgt                             >  2:39346/1‑140 (MQ=255)
                                                                                                                 tACAATCTAACCAAATAATGTTTTTTTGAACTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCCCCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAat                          >  2:285664/1‑141 (MQ=255)
                                                                                                                  aCAATCTAACCAAATAATGTTTTTTTTTTTTTTTTTAAAATTAAATGCCTCTACACTCTAAATAATACCACCACGAGTAGCAATTATTTCTGCAGCGTTAATACATTCTTTCAAACATGAGAAATATAAATTCCttatttt                         >  2:441259/1‑135 (MQ=255)
                                                                                                                  aCAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCCGAAAGGATGATAACTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAata                         >  1:255136/1‑141 (MQ=255)
                                                                                                                  aCAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAatt                         >  1:122478/1‑140 (MQ=255)
                                                                                                                  aCAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAata                         >  1:122459/1‑141 (MQ=255)
                                                                                                                    aaTCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAATTCTAAATTCTTGAAAtttt                       >  2:870672/1‑141 (MQ=255)
                                                                                                                    aaTCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAATAtt                       >  1:998142/1‑141 (MQ=255)
                                                                                                                         aaCCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTGGCAAACATGTGAACTCTAAATTCCTGAAATATTGTTCa                  >  1:639473/1‑141 (MQ=255)
                                                                                                                         aaCCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTATACACTCTCAATAATACCAACACCGGTAGCAATTATATCTGCAGTGTTAATACATTTTTGAAAACATGTGACTTCTAAATTCCTGAAATATTGTTcc                  >  2:984457/1‑140 (MQ=255)
                                                                                                                          aCCAAATAATGTTTTTTTTTTTTTTTAAAAAATTAAATGCCTCTTCACTCTCAATAATATCACCACCAGTAGCAATTATATCTGCATTGTTATTACATTCTTGCAAAAATGTGAACTCTATATTCTTTTAATATTGTTTAc                 >  2:418119/1‑141 (MQ=255)
                                                                                                                          aCCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTCAATTCCTGAAATATTGTTCAc                 <  1:571148/141‑1 (MQ=255)
                                                                                                                          aCCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAACTTCCTGAAATATTGTTCAc                 >  1:417569/1‑141 (MQ=255)
                                                                                                                          aCCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCGTGCAAACATGTGAACTCTAAATTCCTGAAATATTGTTCAc                 >  2:563612/1‑141 (MQ=255)
                                                                                                                              aataatGTTTTTTTGCTTTTTAAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAATATTGTTCACCATg             >  1:221117/1‑141 (MQ=255)
                                                                                                                              aataatGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCTAACATGTGCACTCTAAATTCCTGACATATTGTTAACCAtt             >  2:193603/1‑140 (MQ=255)
                                                                                                                              aataatGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAATATTGTTCACCATg             >  1:457613/1‑141 (MQ=255)
                                                                                                                                  atGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAATATTGTTCACCATGTTTa         >  1:429654/1‑141 (MQ=255)
                                                                                                                                   tGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAATATTGTTCACCATGTTTAc        >  1:35557/1‑141 (MQ=255)
                                                                                                                                    gTTTTTTTGATTTTTCACAAAAGTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAATATTGTTCACCATGTTTACa       >  1:1099466/1‑141 (MQ=255)
                                                                                                                                    gTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAATATTGTTCACCATGTTTACa       >  1:326839/1‑141 (MQ=255)
                                                                                                                                        ttttGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAATATTGTTCACCATGTTTACAAAat   <  1:562541/141‑1 (MQ=255)
                                                                                                                                         tttGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAATATTGTTCACCATGTTTACAAAata  >  1:826554/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                         
ATTATATCTTAAATTCCGCGAGCAATACAAGTCATTTAACTGCTTGATTGTCTTATTATTTGCATTAGGTCGATGTGGGTCATCATTGATTCGACTATATATAATATCAATATTACAATCTAACCAAATAATGTTTTTTTGATTTTTCAAAAAATTAAATGCCTCTTCACTCTCAATAATACCACCACCAGTAGCAATTATATCTGCAGTGTTAATACATTCTTGCAAACATGTGAACTCTAAATTCCTGAAATATTGTTCACCATGTTTACAAAATA  >  USA300TCH1516_ALE/1653550‑1653827

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