Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 299,457 G→T 6.6% *390L (TGA→TTA)  tarF_1 → Teichoic acid glycerol‑phosphate transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE299,4570GT 6.6% 951.1 / 36.9 280*390L (TGA→TTA) tarF_1Teichoic acid glycerol‑phosphate transferase
Reads supporting (aligned to +/- strand):  ref base G (156/105);  new base T (5/14);  total (161/119)
Fisher's exact test for biased strand distribution p-value = 6.83e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 283 total aligned reads displayed.

ATATTGAAGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTTGTG  >  USA300TCH1516_ALE/299367‑299546
                                                                                          |                                                                                         
atatTGAAGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTg                                                                                           >  1:2280925/1‑91 (MQ=255)
atatTGAAGCACGAACGGTATTTTATCAAAAATTCTGTTCAGTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTg                                                                                           >  2:2998338/1‑91 (MQ=255)
 tatTGAAGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGa                                                                                          >  1:2898664/1‑91 (MQ=255)
  atTGAAGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGAt                                                                                         >  2:2461492/1‑91 (MQ=255)
  atTGAAGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGAt                                                                                         <  1:3650461/91‑1 (MQ=255)
   ttGAAGCACGAACGGTATTTTATCAAAGATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTTAtc                                                                                        <  2:2378226/91‑2 (MQ=255)
   ttGAAGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTTAtc                                                                                        <  2:1071467/91‑2 (MQ=255)
   ttGAAGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATa                                                                                        <  1:1285598/91‑1 (MQ=255)
   ttGAAGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATa                                                                                        <  1:1397261/91‑1 (MQ=255)
      aaGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACt                                                                                     >  2:2155945/1‑91 (MQ=255)
       aGCACGAGCGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACta                                                                                    <  2:36053/91‑2 (MQ=255)
       aGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTTATCACtt                                                                                    <  1:3392234/91‑1 (MQ=255)
        gCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTTATCACTTa                                                                                   <  1:187982/91‑1 (MQ=255)
        gCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTa                                                                                   <  2:139895/91‑1 (MQ=255)
        gCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTa                                                                                   >  1:621108/1‑91 (MQ=255)
          aCGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTTATCACTTaaa                                                                                 <  1:2334238/91‑1 (MQ=255)
          aCGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTTATCACTTaaa                                                                                 <  1:3087284/91‑1 (MQ=255)
          aCGAACGGTATTGTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTaaa                                                                                 >  1:3030661/1‑91 (MQ=255)
           cGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTTATCACTTaaaa                                                                                <  1:2140129/91‑1 (MQ=255)
           cGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTaaaa                                                                                >  2:356033/1‑91 (MQ=255)
              aCGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACa                                                                             >  1:3569961/1‑91 (MQ=255)
              aCGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACa                                                                             >  1:220620/1‑91 (MQ=255)
              aCGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACa                                                                             <  2:1106754/91‑1 (MQ=255)
               cGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACaa                                                                            >  1:3626411/1‑91 (MQ=255)
               cGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACaa                                                                            >  1:889256/1‑91 (MQ=255)
               cGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACaa                                                                            >  1:2607689/1‑91 (MQ=255)
                    ttttATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATaaaa                                                                       >  1:1013403/1‑91 (MQ=255)
                        aTCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTCAAGTGATAACTTAAAAACAATAAAAAAtt                                                                   >  1:3176243/1‑91 (MQ=255)
                         tCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATta                                                                  >  2:2364477/1‑91 (MQ=255)
                         tCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATta                                                                  <  1:1229491/91‑1 (MQ=255)
                           aaaaaTTCTGTTCATTAGAAGATGGATAAGCGTCACAACGAATTTGCCAAACGATTTTGAAGTGATAACTTAAAAACAATAAAAAATtata                                                                <  2:3694393/91‑1 (MQ=255)
                           aaaaaTTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATtata                                                                <  1:560390/91‑1 (MQ=255)
                           aaaaaTTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATtata                                                                <  2:478759/91‑1 (MQ=255)
                                ttCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAaatta                                                           >  1:175863/1‑91 (MQ=255)
                                ttCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAaatta                                                           >  1:448566/1‑91 (MQ=255)
                                 tCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAaattaa                                                          >  1:3501090/1‑91 (MQ=255)
                                 tCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAaattaa                                                          >  1:694952/1‑91 (MQ=255)
                                   tGTTCATTAGAAGATGGACAAGCGTCACAACGAACTTGCCAAACGATTTTTAAGAGATAACTTAAAAACAATAAAAAAATATAaaaaaatt                                                        >  2:1582278/1‑91 (MQ=255)
                                    gTTCATTAGAAGATGGACAAGCGTCAGAACGTATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAattaatta                                                       >  2:563671/1‑91 (MQ=255)
                                    gTTCATCAGAAGCTGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAACAAATTATCAattaatta                                                       >  2:3409316/1‑91 (MQ=255)
                                     ttCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAg                                                      >  1:1306329/1‑91 (MQ=255)
                                     ttCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAg                                                      >  2:2886949/1‑91 (MQ=255)
                                     ttCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAg                                                      >  1:995305/1‑91 (MQ=255)
                                      ttaTTTGAATATGGATAATTGTCACAATGAATTTGTTAAATGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGt                                                     <  2:2550939/89‑1 (MQ=255)
                                      tCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGt                                                     <  1:3264549/91‑1 (MQ=255)
                                       tcatAGAAGATAGACAAGCGTCACAACAAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGtt                                                    <  2:3801177/88‑1 (MQ=255)
                                        aTTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTa                                                   >  2:2706077/1‑91 (MQ=255)
                                          tAGAAGGTGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAg                                                 <  1:1690520/91‑1 (MQ=255)
                                              agaTGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGat                                             >  1:3889253/1‑91 (MQ=255)
                                               gaTGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGata                                            >  1:2944572/1‑91 (MQ=255)
                                                aTGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGatat                                           >  1:269243/1‑91 (MQ=255)
                                                 tGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGAtata                                          <  1:1474270/91‑1 (MQ=255)
                                                 tGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGAtata                                          <  1:1080136/91‑1 (MQ=255)
                                                 tGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGAtata                                          >  2:2928287/1‑91 (MQ=255)
                                                 tGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGAtata                                          <  1:2807232/91‑1 (MQ=255)
                                                  ggACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATataa                                         <  1:2454504/91‑1 (MQ=255)
                                                  ggACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATataa                                         >  1:947287/1‑91 (MQ=255)
                                                  ggACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATataa                                         <  2:1494151/91‑1 (MQ=255)
                                                     cAAGCGTCACAACGAATTTGCCAAACGATTGTTAAGTGCTAACTTAAAAACAATAAAAAATGATAAATTAATTAGTCAAGTGATatcaaca                                      >  2:1264720/1‑89 (MQ=255)
                                                       aGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAaataat                                    >  1:3941675/1‑91 (MQ=255)
                                                        gCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAataata                                   >  1:1945466/1‑91 (MQ=255)
                                                        gCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAataata                                   >  1:2972369/1‑91 (MQ=255)
                                                          gTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATaaa                                 <  2:2015154/91‑1 (MQ=255)
                                                            cacaACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACg                               <  1:1038720/91‑1 (MQ=255)
                                                             acaACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGa                              >  2:1851295/1‑91 (MQ=255)
                                                             acaACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGa                              >  1:3711503/1‑91 (MQ=255)
                                                              caACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAAGAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGaa                             <  2:3485430/91‑1 (MQ=255)
                                                              caACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGaa                             >  1:325245/1‑91 (MQ=255)
                                                              caACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGaa                             >  1:216670/1‑91 (MQ=255)
                                                               aaCGAATTTGCCAAACGATTTTTTTTTTTTTTCTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGaaa                            >  2:3751155/1‑91 (MQ=255)
                                                                aCGAATTTGCCGAACGATTTTTCGGTGATAACTTAAAAACAATAAAAAAGTATAAATTATTTAGTTAAGTGATATAAATAATAAACGAAAt                           <  2:3009609/91‑1 (MQ=255)
                                                                 cGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATg                          >  1:505765/1‑91 (MQ=255)
                                                                   aagttGCCAAACGATTTTTAAGTGATAACCGAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGtt                        <  2:2222899/88‑1 (MQ=255)
                                                                    aTTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGttt                       >  1:2650864/1‑91 (MQ=255)
                                                                    aTTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGttt                       >  1:1324521/1‑91 (MQ=255)
                                                                     tttGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTg                      >  2:2184055/1‑91 (MQ=255)
                                                                      ttGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGc                     <  1:1544023/91‑1 (MQ=255)
                                                                      ttGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGc                     >  2:2645748/1‑91 (MQ=255)
                                                                        gCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCtt                   >  1:1814894/1‑91 (MQ=255)
                                                                        gCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCtt                   <  2:1513460/91‑1 (MQ=255)
                                                                        gCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCtt                   >  2:2984161/1‑91 (MQ=255)
                                                                         ccAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTg                  >  1:2537811/1‑91 (MQ=255)
                                                                         ccAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTg                  <  2:188830/91‑1 (MQ=255)
                                                                         ccAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTg                  >  2:3305142/1‑91 (MQ=255)
                                                                               gATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGtt            <  2:3366403/91‑1 (MQ=255)
                                                                                 tttttAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGttat          <  1:704562/91‑1 (MQ=255)
                                                                                 tttttAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGttat          >  1:758746/1‑91 (MQ=255)
                                                                                   tttAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGttatta        >  1:1074490/1‑91 (MQ=255)
                                                                                   tttAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGttatta        <  2:1763955/91‑1 (MQ=255)
                                                                                   tttAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGttatta        >  2:2428618/1‑91 (MQ=255)
                                                                                    ttAAGTTATCACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTtattat       >  2:70387/1‑91 (MQ=255)
                                                                                    ttAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTtattat       <  2:1698167/91‑1 (MQ=255)
                                                                                     tAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTattatt      >  1:2523287/1‑91 (MQ=255)
                                                                                      aaGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTAttt     <  2:3098736/91‑1 (MQ=255)
                                                                                       aGTTATCACTTCAAAACAATAAAAAATTATCTATTAATTAGTGCAGTGATATAAATACTACTCTAAATGTTTGTTTGTTAGTttttttttt    >  2:3523049/1‑90 (MQ=255)
                                                                                       aGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTtg    <  1:376399/91‑1 (MQ=255)
                                                                                       aGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTtg    >  1:1634095/1‑91 (MQ=255)
                                                                                       aGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTtg    <  2:2259697/91‑1 (MQ=255)
                                                                                         tGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTtgtg  <  1:3335539/91‑1 (MQ=255)
                                                                                         tGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTtgtg  <  2:2058901/91‑1 (MQ=255)
                                                                                          |                                                                                         
ATATTGAAGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTTGTG  >  USA300TCH1516_ALE/299367‑299546

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: