Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,742,872 A→G 5.0% V292A (GTG→GCG)  relA ← GTP pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,742,8720AG 5.0% 230.5 / 2.4 80V292A (GTG→GCG) relAGTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base G (2/2);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

CGTTCGTATTTGGATTTCGAGCGGGTCTCCATTTGGGCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATATGTTTAGGTCTACCATTT  >  USA300TCH1516_ALE/1742739‑1743012
                                                                                                                                     |                                                                                                                                            
ctttCGTATTTGTGTTTCAGGCGGGTCTCCATTTGGGCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGCTTTCCATAACGTATGCACCAAccc                                                                                                                                       <  2:623938/139‑1 (MQ=255)
         ttGGATTTCGAGCGGGTCTCCATTTGGGCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGGTTCCATAACGTATGCACCAACTCAAGTATCGc                                                                                                                              <  2:458717/141‑1 (MQ=255)
         ttGGATTCCGCGCGGGTCTCCATTTCAGCCTACTACTGTAGTATGCAATGCGTGACACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGc                                                                                                                              <  1:241750/141‑1 (MQ=255)
         ttGGATGTCGTGCGGGTCCCCATTTGGGCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGc                                                                                                                              <  2:323572/141‑1 (MQ=255)
              tttCGAGCGGGTCTCCATTTGGGCCTACTACTGTAGTATGCAATGACTGATACAAATTGTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCCCAAAAc                                                                                                                         >  1:261911/1‑141 (MQ=255)
               ttCGAGCGGGTCTCCATTTGGGCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACa                                                                                                                        >  1:23996/1‑141 (MQ=255)
                  gAGCGGGTCTCCATTTGGGCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTGCATAACGTATGCCCCAACCCGAGTATCGCATAACAATc                                                                                                                     >  2:334526/1‑141 (MQ=255)
                  gAGCGGGTCTCCATTTGGGCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACaaca                                                                                                                     >  2:980243/1‑139 (MQ=255)
                    gCGGGTCTCCATTTGGGCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCACGTATCGCATAAAAATCAt                                                                                                                   >  1:235031/1‑141 (MQ=255)
                    gCGGGTCTCCATTTGGGCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCAt                                                                                                                   >  1:819166/1‑141 (MQ=255)
                       ggTCTCCATTTGGGCCTACTACTGTAGTATGCAATGACTGATACGAATTTTGTTTAGGCATTGCACTATAATCTTTACAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTaa                                                                                                                <  2:320681/141‑1 (MQ=255)
                        gTCTCCATTTGGGCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTaaa                                                                                                               >  1:924289/1‑140 (MQ=255)
                             cATTTGGGCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAAt                                                                                                          <  2:1102558/141‑1 (MQ=255)
                                  gggCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACa                                                                                                     >  1:996545/1‑141 (MQ=255)
                                  gggCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACa                                                                                                     >  1:439692/1‑141 (MQ=255)
                                   ggCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGTCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACaa                                                                                                    >  1:425216/1‑141 (MQ=255)
                                   ggCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACaa                                                                                                    >  1:693123/1‑141 (MQ=255)
                                           actTTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGt                                                                                            <  2:1078441/141‑1 (MQ=255)
                                           actGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAACATAATCCTCAAAACGTCCTGGCATCGGTTTCCAGAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGt                                                                                            <  2:749462/141‑1 (MQ=255)
                                             tGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATTCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTAt                                                                                          <  2:617208/141‑1 (MQ=255)
                                             tGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTAt                                                                                          <  2:1010970/141‑1 (MQ=255)
                                              gtagtaTGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATc                                                                                         >  1:351462/1‑141 (MQ=255)
                                                     gCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCCCATAACAATAATTAATAGAATTGAAAATAACCGGTATCGCcaaca                                                                                  >  2:474241/1‑141 (MQ=255)
                                                       aaTGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACaaa                                                                                <  2:829625/141‑1 (MQ=255)
                                                       aaTGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACaaa                                                                                <  2:124085/141‑1 (MQ=255)
                                                           aCTGATACAAATTTTGTTTAGGCATTGAAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCaa                                                                            <  2:380661/141‑1 (MQ=255)
                                                            cTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATCACACGTATCGCCAACAAATCaaa                                                                           >  2:104151/1‑141 (MQ=255)
                                                               aTACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAAAAATAATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAAt                                                                        <  2:281013/141‑1 (MQ=255)
                                                                  cAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCGTAACAATGATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTg                                                                     >  1:816317/1‑141 (MQ=255)
                                                                  cAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTg                                                                     >  1:575844/1‑141 (MQ=255)
                                                                   aaaTTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGa                                                                    >  1:948449/1‑141 (MQ=255)
                                                                         tttttCGGGATTGAAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGTACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAAt                                                              <  2:325145/139‑1 (MQ=255)
                                                                         tttttAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAAt                                                              <  2:971602/139‑1 (MQ=255)
                                                                         tGTTTAGGCATTGCAATATAATCTTTAAACCGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAAt                                                              <  2:974356/141‑1 (MQ=255)
                                                                         tGTTTAGGCATTGCAATATAATCCTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAAt                                                              <  2:184563/141‑1 (MQ=255)
                                                                               ggCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCGCCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGtttt                                                        >  1:104821/1‑141 (MQ=255)
                                                                               ggCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCGCCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGtttt                                                        >  2:850405/1‑141 (MQ=255)
                                                                               ggCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGtttt                                                        >  1:988914/1‑141 (MQ=255)
                                                                                gCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTACCATAACCTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGttttt                                                       >  2:690689/1‑141 (MQ=255)
                                                                                    ggCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCGAAAATTTGATCAAATTGttttttttt                                                   >  2:459589/2‑138 (MQ=255)
                                                                                       aaTATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCAGCGACCCAAGTATCGCGTAACAATCACAAATAGAATTGACAATAACACGTATGGCCAACAAATGAAAAATGTGATCAAATGGttatttgggcct                                                >  2:384576/1‑135 (MQ=255)
                                                                                       aaTATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGAAAAAAAATAAAAAATTTGATCAAATTGTTTTTTCTGCtt                                                <  2:23996/141‑1 (MQ=255)
                                                                                         ttaaaTCTTTAAAACCGCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTtca                                              <  2:97937/138‑1 (MQ=255)
                                                                                         tataGTCTTTCAAAGGTATTTGCATCGGTTTCCTTAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTtca                                              <  1:541810/141‑1 (MQ=255)
                                                                                         tataATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTtca                                              >  1:283339/1‑141 (MQ=255)
                                                                                         tataATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTtca                                              >  1:1066935/1‑141 (MQ=255)
                                                                                          ataATCTGTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTtcat                                             <  2:832598/141‑1 (MQ=255)
                                                                                           taATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACGAATCAAAAATTTGTTCAAATTGTTTTTTCTGCTtcatc                                            <  1:708567/141‑1 (MQ=255)
                                                                                           taATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTtcatc                                            <  2:231052/141‑1 (MQ=255)
                                                                                           taATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTtcatc                                            <  2:347022/141‑1 (MQ=255)
                                                                                           taATCTTGTACACGTCCTGGCATCGGTTTCCATAACGGATGCACCATTCCTAGTATCGCAGTACAATCAGTAATAGATTTGACAATTACTCGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTtcatc                                            <  2:485412/141‑1 (MQ=255)
                                                                                             aTCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGAATAACAATCATTTATAGAATTGACGATTACTCTTATCGCCAACAGATCAAAAATTTGATCAAATTGTTTTTTCTGCTTcatcat                                          <  2:369662/141‑1 (MQ=255)
                                                                                             aTCTTTAAAACCTCCTTGCATCGGTTTCCATAACGTATGCACAAACCCAAGTATCGCATAACAATCATTAATAGAATTGATAATAACACGTATCGCCAACAACTCAAAAATTTGATCAAATTGTTTTTTCTGCTTcatcat                                          <  2:451121/141‑1 (MQ=255)
                                                                                                      aCGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATa                                 >  2:1031570/1‑141 (MQ=255)
                                                                                                      aCGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATa                                 >  2:948980/1‑141 (MQ=255)
                                                                                                      aCGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATa                                 >  1:956270/1‑141 (MQ=255)
                                                                                                        gTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTCTAGCATAACAATCATTAATAGAATTGACAATAACACGGATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTACTTCATCATTTTCCCATaaa                               >  1:532993/1‑141 (MQ=255)
                                                                                                         tCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAAACATAAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTACATCAAATTGTTTTTTCTGCCTCATCATTTTCCGATaccc                              >  1:614321/1‑138 (MQ=255)
                                                                                                           cTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATaa                            >  1:792899/1‑140 (MQ=255)
                                                                                                           cTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATAc                            <  2:239375/141‑1 (MQ=255)
                                                                                                            tGGCATCTGTTTCCATAACGTATGCACCAACCCCAGTATCGCATAACAATCATTAATATAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACt                           <  2:395386/141‑1 (MQ=255)
                                                                                                            tGGCATCGGTTTCCATAACGTATGCACCAACCCCAGTAGCGCACTACAATTATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACt                           <  2:312189/141‑1 (MQ=255)
                                                                                                            tGGCATCGGTTTCCATAACGTATGCACCAACCAAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACt                           >  1:156088/1‑141 (MQ=255)
                                                                                                             ggCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTg                          >  1:977603/1‑141 (MQ=255)
                                                                                                             ggCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTg                          >  1:967514/1‑141 (MQ=255)
                                                                                                             ggCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAAAACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTg                          >  1:43975/1‑141 (MQ=255)
                                                                                                                 tCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAGCAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTaaa                      <  2:561887/141‑1 (MQ=255)
                                                                                                                 ccGGTTTCCATTACGGACGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTAGCGCCAACAAATCAAATATTTGATAAAATTTTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTaaa                      <  1:1033380/140‑1 (MQ=255)
                                                                                                                     tttCCATAACGTATGCGCCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATATg                  <  1:850405/141‑1 (MQ=255)
                                                                                                                     tttCCATAACGTATGCGCCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATATg                  <  2:104821/141‑1 (MQ=255)
                                                                                                                         cATAACTTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAAGTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATATGTTTa              <  2:622916/141‑1 (MQ=255)
                                                                                                                         cATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATATGTTTa              >  1:37548/1‑141 (MQ=255)
                                                                                                                             aCGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATATGTTTAGGTc          >  1:1010664/1‑141 (MQ=255)
                                                                                                                              cGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATATGTTTATGTCt         >  1:581236/1‑141 (MQ=255)
                                                                                                                              cGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATATGTTTAGGTCt         >  2:738097/1‑141 (MQ=255)
                                                                                                                              cGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATATGTTTAGGTCt         >  1:357437/1‑141 (MQ=255)
                                                                                                                                 aTGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATCTTTTTAGGTTTAca      >  1:515076/1‑140 (MQ=255)
                                                                                                                                  tGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGGCAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATATGTTTAGGTCTAcct     >  1:863594/1‑140 (MQ=255)
                                                                                                                                  tGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCGTCATCATTTTCCGATAAATACTGTAAATATGTTTAGGTCTACCa     >  1:144657/1‑141 (MQ=255)
                                                                                                                                   gCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATATGTTTAGGTCTACCAt    >  1:408577/1‑141 (MQ=255)
                                                                                                                                     aCCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATATGTTTAGGTCTACCAttt  <  1:308873/141‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                            
CGTTCGTATTTGGATTTCGAGCGGGTCTCCATTTGGGCCTACTACTGTAGTATGCAATGACTGATACAAATTTTGTTTAGGCATTGCAATATAATCTTTAAAACGTCCTGGCATCGGTTTCCATAACGTATGCACCAACCCAAGTATCGCATAACAATCATTAATAGAATTGACAATAACACGTATCGCCAACAAATCAAAAATTTGATCAAATTGTTTTTTCTGCTTCATCATTTTCCGATAAATACTGTAAATATGTTTAGGTCTACCATTT  >  USA300TCH1516_ALE/1742739‑1743012

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: ‑