Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 76,969 G→A 31.5% intergenic (+509/‑107) USA300HOU_RS00140 → / → USA300HOU_RS00350 hypothetical protein/Monoacylglycerol lipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE76,9690GA31.5% 110.3 / 52.4 80intergenic (+509/‑107)USA300HOU_RS00140/USA300HOU_RS00350hypothetical protein/Monoacylglycerol lipase
Reads supporting (aligned to +/- strand):  ref base G (28/26);  new base A (11/14);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCAAA  >  USA300TCH1516_ALE/76831‑77104
                                                                                                                                          |                                                                                                                                       
tctTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAggg                                                                                                                                       <  2:681043/141‑1 (MQ=255)
 ctTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGt                                                                                                                                      >  1:1132939/1‑141 (MQ=255)
   tCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTa                                                                                                                                    <  1:675293/141‑1 (MQ=255)
     aGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCATAAAGGGTTAAt                                                                                                                                  >  2:640648/1‑141 (MQ=255)
      gATCCACGGGTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAAtt                                                                                                                                 <  2:179916/141‑1 (MQ=255)
           aCGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAGTTtatat                                                                                                                            >  1:137345/1‑141 (MQ=255)
             ggATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAg                                                                                                                          <  1:101584/141‑1 (MQ=255)
                  ttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTa                                                                                                                     <  1:406450/141‑1 (MQ=255)
                      ccATGCCACGAAATTAGCATCATGCTAGCACGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTACAGGGTTAATTTATCTAGTATTACTGc                                                                                                                 >  1:44399/1‑141 (MQ=255)
                      ccATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACTGAACCTATAAAAAACATGGGCTAAGAAGGGGAATTTATATAGTATTACTGc                                                                                                                 <  2:849143/141‑1 (MQ=17)
                           ccACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACa                                                                                                            <  2:170193/141‑1 (MQ=255)
                           ccACGAAATTAGCACCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACa                                                                                                            <  2:322442/141‑1 (MQ=255)
                            cACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACAt                                                                                                           >  2:1890/1‑141 (MQ=255)
                               gAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGa                                                                                                        >  2:456123/1‑141 (MQ=255)
                               gAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTAATGCAAACATTGa                                                                                                        >  1:308954/1‑141 (MQ=255)
                                     aGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTggg                                                                                                  >  2:742168/1‑141 (MQ=255)
                                     aGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAAAATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTggg                                                                                                  >  1:1118985/1‑141 (MQ=255)
                                      gCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGt                                                                                                 >  1:521323/1‑141 (MQ=255)
                                           atGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCaaaa                                                                                            <  1:129073/141‑1 (MQ=255)
                                                 gCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAgg                                                                                      <  1:722797/141‑1 (MQ=255)
                                                  cAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATAGTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAggg                                                                                     >  2:638439/1‑141 (MQ=255)
                                                   aaGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAgggg                                                                                    <  1:435626/141‑1 (MQ=255)
                                                   aaGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAgggg                                                                                    >  1:967819/1‑141 (MQ=255)
                                                    aGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGt                                                                                   <  1:493765/141‑1 (MQ=255)
                                                    aGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGt                                                                                   <  1:719774/141‑1 (MQ=255)
                                                         aGCGCACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTAAATCTAGATCGTGACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaa                                                                              <  1:865460/141‑1 (MQ=37)
                                                         aGCGAACACTGACAGGATAAATTCGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTCCTGTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaa                                                                              <  1:169454/141‑1 (MQ=37)
                                                          gccaaCAATGCCACGATACATTATGGGTAAGCTATATTTTTTTCCTTTGCAACAGAACCCATAAAAAACATGGGATTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaa                                                                             <  1:10818/138‑1 (MQ=255)
                                                          gCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaa                                                                             >  1:357915/1‑141 (MQ=255)
                                                                 cTGACGTGCTAAATTAGTAGTTAGCTATATTTTCTTACTTTGCAACAGCATCTATAAAAAACATGGGCTTAAAGGTTTAATTTATATAGTATTACTGCTAACATTGATTTGGGTCTAAAATTAGGGGTATTAAAAAactga                                                                      <  2:506531/141‑1 (MQ=255)
                                                                  tGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaatgaa                                                                     >  2:468751/1‑141 (MQ=255)
                                                                    aCATGATAAATTAGTGGTTAGCGTTATTTTTTTACTTTGTAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaatgaatg                                                                   <  2:308185/141‑1 (MQ=255)
                                                                      aggATAAATTAGTGGGTAGCTATATGTTTTTACGTTGCAAAAGCAGCTATGAAAAAGGTGGGATGAAGAGGTTAATTTATATATTGTTACTGCTAACATTGAGTTGGGTTAAAAATTAGGGGTATTAAAAAAAtgaatgaa                                                                 <  2:1065413/139‑1 (MQ=255)
                                                                      aTGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTAATTTGGGTCAAAAATTAGGGGTATTAAAAAAAtgaatgaa                                                                 <  1:1014505/141‑1 (MQ=255)
                                                                      aTGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAAtgaatgaa                                                                 <  2:321974/141‑1 (MQ=255)
                                                                      aTGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAAtgaatgaa                                                                 <  1:1100603/141‑1 (MQ=255)
                                                                           aaaTTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATa                                                            >  1:695291/1‑141 (MQ=255)
                                                                            aaTTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTACATAGTATTACTGCAAACATTGCTTTGCATCAAGAATTAGGGGTATTAAAAAAATGAATGAATCATaa                                                           >  1:1117860/1‑141 (MQ=38)
                                                                                        gCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCt                                               <  1:1890/141‑1 (MQ=255)
                                                                                         ccatatTTTTTTACTTTGCAACAGAACCTATAAAAGACATGGGCTTTAACGGTTACTTTAAATAGTATAACTGCAAACATTGATTTGGGTCAAAAATTAAGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTg                                              >  2:126561/3‑141 (MQ=255)
                                                                                         cTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGAATTAAAAAAATGAATGAATCATAATTTTATCAAGCTg                                              >  2:154990/1‑141 (MQ=255)
                                                                                          tatataTTTTTACTTTGCAACAGAAAATATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGa                                             <  2:769259/141‑1 (MQ=255)
                                                                                             atTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTg                                          >  1:138835/1‑141 (MQ=255)
                                                                                             atTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGATGATTg                                          >  2:783033/1‑141 (MQ=255)
                                                                                                  tttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGagag                                     >  2:176541/1‑141 (MQ=255)
                                                                                                  tttACCGTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGagag                                     <  1:952825/141‑1 (MQ=255)
                                                                                                   ttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAgg                                    >  1:245634/1‑141 (MQ=255)
                                                                                                   ttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAgg                                    <  1:742168/141‑1 (MQ=255)
                                                                                                   ttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAgg                                    <  2:1079991/141‑1 (MQ=255)
                                                                                                    tACTTTGCAACACAACCTATAAAAAACATGGGCATAAAGGGTTAATTTAGATAGTATTACTGCAAACATTGATTTGGTTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGt                                   <  2:670969/141‑1 (MQ=255)
                                                                                                     tcTTTGCAACAGAATCTCTAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGtt                                  <  2:1115145/140‑1 (MQ=255)
                                                                                                     acattGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGtt                                  <  2:643406/138‑1 (MQ=255)
                                                                                                     aCTTTGCAACAGAACCTATAACAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGtt                                  <  2:945375/141‑1 (MQ=255)
                                                                                                     aCTTTGCAACAGAACCTATAAAAAACATGGGCTGAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGtt                                  <  2:1115287/141‑1 (MQ=255)
                                                                                                      ccttGCAACAGAACCTATCACAAACATGGGCTTAAAGGGTTAATTTAGATAGTATAAATACAAACATGTATTTGGGTCAAAAACTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTa                                 <  1:276109/139‑1 (MQ=255)
                                                                                                      cTTTGCAACAGAACCTATAAAAAACATGGGCGTAAAGGGTTAATTTATATAGTATTATTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTa                                 >  1:1123785/1‑141 (MQ=255)
                                                                                                        ttGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAACTGAGTGGTGAGGtcaga                               >  1:1027909/1‑137 (MQ=255)
                                                                                                        ttGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATGTTATCAAGCAGATTGGAGAGGATaaa                               >  2:795695/1‑141 (MQ=255)
                                                                                                        ttGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTaaa                               >  2:754661/1‑141 (MQ=255)
                                                                                                        ttGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTaaa                               >  1:672501/1‑141 (MQ=255)
                                                                                                         tGCAGCAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGGCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATGTTATCAAGCTGATTGGAGAGGTTaaaa                              >  2:326091/1‑141 (MQ=255)
                                                                                                            aaCAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGc                           <  2:280233/141‑1 (MQ=255)
                                                                                                               aGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCAtt                        >  2:139066/1‑141 (MQ=255)
                                                                                                                gAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATta                       <  1:1097785/141‑1 (MQ=255)
                                                                                                                gAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATta                       <  2:1118733/141‑1 (MQ=255)
                                                                                                                gAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATta                       >  2:39368/1‑141 (MQ=255)
                                                                                                                  aCCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATtata                     >  1:1102714/1‑141 (MQ=255)
                                                                                                                    cTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATtatata                   >  2:395481/1‑141 (MQ=255)
                                                                                                                     tataAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATaa                  >  1:784558/1‑141 (MQ=255)
                                                                                                                          aaaaCATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAAttt             <  1:106018/141‑1 (MQ=255)
                                                                                                                            aaCATAGACTAAAAGTGATAATTTATATAGTATTACTGCAAACATTGATTTGAATAAAAAAAAAAAAAAAAAAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAAtataa           <  2:78052/141‑2 (MQ=255)
                                                                                                                             aCATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATAAAAAAAATTAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTAtt          >  1:276692/1‑141 (MQ=255)
                                                                                                                                tgggtgtAAAGGGTTAATTTATATCGTATTACTGCAAACAGTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAg       <  2:182115/135‑1 (MQ=255)
                                                                                                                                 gggCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGt      >  2:1027029/1‑141 (MQ=255)
                                                                                                                                 gggCTTAAAAGGTTAATTTATATAGTATTACGGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGt      <  1:139066/141‑1 (MQ=255)
                                                                                                                                  ggCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGtg     >  1:694183/1‑140 (MQ=255)
                                                                                                                                  ggCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTc     >  1:32823/1‑141 (MQ=255)
                                                                                                                                  ggCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTc     <  2:967819/141‑1 (MQ=255)
                                                                                                                                  ggCTTAAAAGGTTAATTTATATAGAATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGGTTGGAGAGGTTAAAATGCATTATAGAAAAGTTATTAAGta     >  2:1083935/1‑140 (MQ=255)
                                                                                                                                   gCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCa    >  1:184407/1‑141 (MQ=255)
                                                                                                                                     ttAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCaaa  >  1:442960/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                       
TCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCAAA  >  USA300TCH1516_ALE/76831‑77104

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: