Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 297,176 C→T 6.1% S337L (TCA→TTA)  tarK → Teichoic acid ribitol‑phosphate polymerase TarK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE297,1760CT 6.1% 195.7 / 2.6 68S337L (TCA→TTA) tarKTeichoic acid ribitol‑phosphate polymerase TarK
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base T (2/2);  total (32/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 = 4.00e-01

CTTTAGACCATCACAAGAAATTATTCAAGTGTGGCATGCTGTTGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAAGAGATGGAAGATGAAT  >  USA300TCH1516_ALE/297036‑297316
                                                                                                                                            |                                                                                                                                            
cTTTAGACCATCACAAGAAATTATTCAAGTGTGGCATGCTGTGGGTGCCTTTAAAACAGTTGGCTTTAGTCGCACATGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTtc                                                                                                                                              <  1:105406/141‑1 (MQ=255)
  ttAGACCATCACAAGAACTTATTCAAGTGTGGCATGCTGGTGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTtcat                                                                                                                                            <  2:640789/141‑1 (MQ=255)
  ttAGACCATCACAAGAAATTATTCAAGTGTGGCATGCTGTTGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGGAGATGTTtcat                                                                                                                                            >  2:369387/1‑141 (MQ=255)
        cATCACAAGAAATTATTCAAGTGTGGCATGCTGTTGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGaaa                                                                                                                                      >  2:483400/1‑141 (MQ=255)
            accaGAAAGTATTCAAGTGTGGCATGCTGTTGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCTATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGa                                                                                                                                  <  2:20271/138‑1 (MQ=255)
               aGAAATTATTCAAGCGTGGCATGCTGTTGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGATACCGttat                                                                                                                               >  1:402667/1‑141 (MQ=255)
                         cAAGTGTGGCATGCTGTTGTTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTTATCAGAAACCGATATTCCATTTTAt                                                                                                                     >  2:199782/1‑141 (MQ=255)
                         cAAGTGTGGCATGCTGTTGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTTATCAGAAACCGATATTCCATTTTAt                                                                                                                     >  2:648918/1‑141 (MQ=255)
                         cAAGTGTGGCATGCTGTTGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTAt                                                                                                                     >  1:615029/1‑141 (MQ=255)
                          aaGTGTGGCATGCTGTTGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATg                                                                                                                    <  1:135829/141‑1 (MQ=255)
                          aaGTGTGGCATGCTGTTGGGCCTTTTAAAACAGTTGTCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTTATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATg                                                                                                                    <  1:391370/141‑1 (MQ=255)
                                ggCCTGCTGTTGGTGCCTTGAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAg                                                                                                              <  1:196596/141‑1 (MQ=255)
                                  cATGCTGTTGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCa                                                                                                            >  2:427104/1‑141 (MQ=255)
                                   aTGCTGTTGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCAt                                                                                                           <  1:693764/141‑1 (MQ=255)
                                   aTGCTGTTCGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAACCCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCAt                                                                                                           <  1:193476/141‑1 (MQ=255)
                                     gCTGTTGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCAttt                                                                                                         <  2:513724/141‑1 (MQ=255)
                                         ttGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGaa                                                                                                     <  2:120838/141‑1 (MQ=255)
                                             tGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCGTTAAACCATCGTAGTTACACGAAAACATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCTTTTGGAATTag                                                                                                 >  1:322409/1‑141 (MQ=255)
                                             tGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTag                                                                                                 <  2:661847/141‑1 (MQ=255)
                                               ccTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTGACATCAATCCTGAAGCATTTGGAATTagag                                                                                               <  1:768636/141‑1 (MQ=21)
                                                 tttAAAACAGTTGGCTTTAGTAGTACAGGTAAAAAAGGTGGTCCGTTAATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGACACCGACATTCCATTTTATGCTGAAGGATTTGGAATTatggaa                                                                                             >  2:564603/1‑136 (MQ=21)
                                                   tAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGagaaga                                                                                           <  2:255454/141‑1 (MQ=255)
                                                     aaaCAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGagaagaca                                                                                         >  2:237855/1‑139 (MQ=255)
                                                        cAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAAtg                                                                                      <  1:306794/141‑1 (MQ=255)
                                                        cAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAAtg                                                                                      <  2:122057/141‑1 (MQ=255)
                                                        cAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAAtg                                                                                      >  1:91894/1‑141 (MQ=255)
                                                           gtGGCTTTAGTCGTACAGGTAAAAAAGTTGGTCCGTTTATGGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAAtgttg                                                                                   <  2:184254/140‑1 (MQ=255)
                                                             ggCTTTAGTCGTACAGGTAACAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTa                                                                                 >  1:489079/1‑141 (MQ=255)
                                                               cTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTAcc                                                                               >  1:411827/1‑141 (MQ=255)
                                                                 ttAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTAATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCaa                                                                             >  2:62913/1‑141 (MQ=255)
                                                                         aCAGGTAAAAAAGGTGGTCCGTTTATCGCTTCATTAAACCATCGTAGTTACACGAAAGCAGATGTTTCATCAGAAACCGATGTTCCATTTTATGCTGACGCATTGGGAATAACAGACGAAAATGTTGTCCTAACACgtgtc                                                                     >  2:271705/1‑140 (MQ=255)
                                                                             gTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCAc                                                                 <  1:650161/141‑1 (MQ=255)
                                                                                aaaaaGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTa                                                              <  2:587523/141‑1 (MQ=255)
                                                                                 aaaaGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTGCCACGTAc                                                             >  1:405292/1‑141 (MQ=255)
                                                                                     ggtggtCCGTTTATCGATTCGTTAAACCATCGTAGTTACACGAACGCATATGATTCATGAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGAt                                                         >  1:471916/1‑141 (MQ=255)
                                                                                       tgggCCCTTTATCGATTCATTAAACCATCCTAGTTACACGAAAGCATATGTTACATCTGGCACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGt                                                       <  1:366742/141‑1 (MQ=255)
                                                                                          tCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACt                                                    <  2:667708/141‑1 (MQ=255)
                                                                                               ttATCTATTTATTAAACCATCTTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCACCATGTGTACCACGTACTGATGTACTATTtg                                               <  2:354387/141‑1 (MQ=25)
                                                                                                 aTCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACAAACAGGTGTACAACGTAATGATGTACTATTtgtt                                             >  1:544177/1‑139 (MQ=255)
                                                                                                     aTTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAg                                         <  1:262253/141‑1 (MQ=255)
                                                                                                            aaaCCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGCTGAAGGTTATGc                                  >  2:117599/1‑141 (MQ=255)
                                                                                                            aaaCCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGc                                  <  1:698699/141‑1 (MQ=255)
                                                                                                             aaCCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCa                                 <  2:489079/141‑1 (MQ=255)
                                                                                                              aCCATCGTAGTTACACGAAAGCATATGTTTTATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCaa                                <  1:199782/141‑1 (MQ=255)
                                                                                                              aCCATCGTAGTTACACGAAAGCATATGTTTTATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCaa                                <  1:648918/141‑1 (MQ=255)
                                                                                                              aCCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCaa                                <  2:592154/141‑1 (MQ=255)
                                                                                                              aCCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCaa                                >  2:484330/1‑141 (MQ=255)
                                                                                                              aCCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATAAAGCTTATGCaa                                >  1:147595/1‑141 (MQ=255)
                                                                                                                  tCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTAATGCAAcaca                            >  2:441633/1‑141 (MQ=255)
                                                                                                                   cGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACaa                           >  1:250635/1‑141 (MQ=255)
                                                                                                                    gTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGCTGTACTATTTGATGACTGTTATGCCTCACaaa                          >  2:726681/1‑141 (MQ=255)
                                                                                                                     tAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGGTGAAGCTTATGCAACACAAAt                         >  1:197129/1‑141 (MQ=255)
                                                                                                                        ttACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTaa                      >  1:588504/1‑141 (MQ=255)
                                                                                                                        ttACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTaa                      >  2:374974/1‑141 (MQ=255)
                                                                                                                        tgacaacaaaGCATATGTTTCATCAGAAACGGATATTCCATTTTATGCTGAAGCATTTTGAATTAGAGAAGAAAATGTTGTACCAACAGGCGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTaa                      <  2:93500/134‑1 (MQ=21)
                                                                                                                         tACACGAAAGTATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTaaa                     <  2:456816/141‑1 (MQ=255)
                                                                                                                         tACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGACGCATCTGGAATTAGAGACGAAACGCTTGTACCAACAGGTGTACCACGTACTGATGTACTAGTTGATGAAGCTTATGCCACACAAATCaaa                     >  1:763914/1‑141 (MQ=255)
                                                                                                                         tACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTaaa                     <  2:541405/141‑1 (MQ=255)
                                                                                                                              gAAAGCGTATGTTTCATCAGAAACCGATATTGCATTTTACGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAaga                <  2:21978/141‑1 (MQ=255)
                                                                                                                              gAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAaga                <  2:649851/141‑1 (MQ=255)
                                                                                                                              gAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAaga                <  2:662369/141‑1 (MQ=255)
                                                                                                                                 aGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAAgagag             >  2:253478/1‑140 (MQ=255)
                                                                                                                                 aGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAAGAGAt             <  1:89649/141‑1 (MQ=255)
                                                                                                                                   gatatGTTTCATCAGGAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAAGAGATgg           <  1:600595/140‑1 (MQ=255)
                                                                                                                                   cATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAAGAGATgg           >  1:777426/1‑141 (MQ=255)
                                                                                                                                     tatGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAAGAGATGGaa         >  2:485935/1‑141 (MQ=255)
                                                                                                                                     tatGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATAAAACAAGAGATGGaa         >  2:698677/1‑141 (MQ=255)
                                                                                                                                     tatGGTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAAGAGATGGaa         >  2:91590/1‑141 (MQ=255)
                                                                                                                                      atGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAAGAGATGGAAg        <  1:240204/141‑1 (MQ=255)
                                                                                                                                       tGTTTGATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAACATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAAGAGATGGAAGa       <  1:62913/141‑1 (MQ=255)
                                                                                                                                         tttCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAAGAGATGGAAGATg     <  2:692064/141‑1 (MQ=255)
                                                                                                                                            catcaGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACAAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGATTATGCAGCAAAAATTAAACAAGAGAGGGAAGATGAAt  >  2:214022/1‑141 (MQ=21)
                                                                                                                                            |                                                                                                                                            
CTTTAGACCATCACAAGAAATTATTCAAGTGTGGCATGCTGTTGGTGCCTTTAAAACAGTTGGCTTTAGTCGTACAGGTAAAAAAGGTGGTCCGTTTATCGATTCATTAAACCATCGTAGTTACACGAAAGCATATGTTTCATCAGAAACCGATATTCCATTTTATGCTGAAGCATTTGGAATTAGAGAAGAAAATGTTGTACCAACAGGTGTACCACGTACTGATGTACTATTTGATGAAGCTTATGCAACACAAATTAAACAAGAGATGGAAGATGAAT  >  USA300TCH1516_ALE/297036‑297316

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