Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 984,359 A→G 5.9% K101E (AAA→GAA)  USA300HOU_RS04865 → 65 kDa membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE984,3590AG 5.9% 255.6 / 2.8 82K101E (AAA→GAA) USA300HOU_RS0486565 kDa membrane protein
Reads supporting (aligned to +/- strand):  ref base A (42/34);  new base G (2/3);  total (45/37)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATGGCATTATGGCTTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGATATTAATGTAAAAACTAAATAAACACT  >  USA300TCH1516_ALE/984220‑984498
                                                                                                                                           |                                                                                                                                           
aGATGGCATTATGGCTTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCaaa                                                                                                                                            >  2:473377/1‑141 (MQ=255)
aGATGGCATTATGGCTTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCaaa                                                                                                                                            >  2:570096/1‑141 (MQ=255)
  aTGGCATTATGGCTTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAAt                                                                                                                                          >  1:155964/1‑141 (MQ=255)
     gaaTTATGGCTTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTg                                                                                                                                       <  1:341563/139‑1 (MQ=255)
     gCATTATGGCTTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTg                                                                                                                                       >  2:68993/1‑141 (MQ=255)
     gCATTATGGCTTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTg                                                                                                                                       >  1:649217/1‑141 (MQ=255)
     gCATTATGGCCTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTg                                                                                                                                       <  1:49005/141‑1 (MQ=255)
         tATGGTTTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGt                                                                                                                                   <  2:722962/141‑1 (MQ=255)
         tATGGCTTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGt                                                                                                                                   <  1:762894/141‑1 (MQ=255)
           tGGCTTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGttt                                                                                                                                 <  1:433521/141‑1 (MQ=255)
            ggCTTTCAATCAATCTTACTTAAACTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTa                                                                                                                                <  2:341610/141‑1 (MQ=255)
             gCTTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAAGCCAAAATCTGCCGTTTac                                                                                                                               >  2:719178/1‑141 (MQ=255)
              cTTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTaca                                                                                                                              <  2:437428/141‑1 (MQ=255)
               tttCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTacac                                                                                                                             <  2:567694/141‑1 (MQ=255)
                 tcaatcaaTCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGGAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGa                                                                                                                           >  2:394172/1‑141 (MQ=255)
                   aatcaatcTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAGAATCTGCCGTTTACACGAtt                                                                                                                         <  1:467654/141‑1 (MQ=255)
                      caatcTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACt                                                                                                                      >  1:650808/1‑141 (MQ=255)
                      caatcTTACTTAAATTTACCAAAAGATAGCCAATCATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACt                                                                                                                      >  1:527788/1‑141 (MQ=255)
                       aatcTTACTTAAATTTACCACAGGTTAAGCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACtt                                                                                                                     >  1:758019/1‑141 (MQ=255)
                        atcTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTg                                                                                                                    <  2:20566/141‑1 (MQ=255)
                        atcTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTg                                                                                                                    >  2:582347/1‑141 (MQ=255)
                            tacttaAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGaaa                                                                                                                <  1:128486/141‑1 (MQ=255)
                             acttaAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAg                                                                                                               >  2:710390/1‑141 (MQ=255)
                               ttaAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGAt                                                                                                             <  1:570096/141‑1 (MQ=255)
                                taAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGaag                                                                                                            >  2:3259/1‑139 (MQ=255)
                                  aagtgagcAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGgt                                                                                                          <  2:569757/134‑1 (MQ=255)
                                  aaTTTACCAAAAGATAGCCAATTCTCATTAGTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATgga                                                                                                          >  2:543269/1‑140 (MQ=255)
                                       aCCAAAAGATAGCCAATTTTCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTaa                                                                                                     <  2:42851/141‑1 (MQ=255)
                                        ccAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTaaa                                                                                                    >  1:379289/1‑141 (MQ=255)
                                        ccAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTaaa                                                                                                    >  1:765557/1‑141 (MQ=255)
                                          aaaaGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTaaaaa                                                                                                  <  1:624748/141‑1 (MQ=255)
                                           aaaGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGGAAAAGCTGCCGTTTACACGCTTACTTGGAAAGATGGTAGTaaaaaa                                                                                                 >  2:42767/1‑141 (MQ=255)
                                           aaaGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTaaaaaa                                                                                                 <  1:134054/141‑1 (MQ=255)
                                              gATAGCCAAGTATCATATTCAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAaagaa                                                                                              <  1:220922/141‑1 (MQ=255)
                                                  gCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTgg                                                                                          >  2:641951/1‑141 (MQ=255)
                                                            atatTTAGATTTAGGACATAAAGTTAAAGCTTTGGCATATGATGAAAGCGGCGTAACACCTGAGAAGATTCGAGATGCACAATCTTCAGTGTACACGATTACTTGGCAAGACGGTGGTAAAAAATAAGTGGATCTTacgaa                                                                                >  1:220008/1‑141 (MQ=255)
                                                               tttagatttagGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTaagaaaga                                                                             >  2:403362/1‑141 (MQ=255)
                                                                 tagatttagGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATa                                                                           >  2:613226/1‑141 (MQ=255)
                                                                     tttagGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTa                                                                       >  2:280626/1‑141 (MQ=255)
                                                                         ggAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTAcaca                                                                   >  1:436024/1‑141 (MQ=255)
                                                                          gAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAg                                                                  <  1:124816/141‑1 (MQ=255)
                                                                           aaaTAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAAAAATTTTTAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGc                                                                 >  1:149144/1‑141 (MQ=255)
                                                                            aaTAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAGAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCa                                                                >  2:287359/1‑141 (MQ=255)
                                                                            aaTAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAGAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTAACCAGCa                                                                >  2:717545/1‑141 (MQ=255)
                                                                             aTAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAAGCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGTAGTGGATATTAAGAAAGATAGCTACACAGCaa                                                               <  2:657085/141‑1 (MQ=255)
                                                                              tAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCaaa                                                              >  2:276010/1‑141 (MQ=255)
                                                                                aaGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGGTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATATTAAGAAAGATAGCTCAACAGCAAACt                                                            >  2:142973/1‑141 (MQ=255)
                                                                                  gTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTg                                                          >  1:613495/1‑141 (MQ=255)
                                                                                  gTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTg                                                          >  2:715578/1‑141 (MQ=255)
                                                                                  gTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTg                                                          >  1:118264/1‑141 (MQ=255)
                                                                                  gTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACGTGGAAAGAGGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTg                                                          >  2:370628/1‑141 (MQ=255)
                                                                                   ttAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGt                                                         <  2:194299/141‑1 (MQ=255)
                                                                                     aaaGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGttt                                                       <  1:132800/141‑1 (MQ=255)
                                                                                     aaaGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGttt                                                       <  1:511186/141‑1 (MQ=255)
                                                                                     aaaGCTTTGCTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGttt                                                       >  2:9843/1‑141 (MQ=255)
                                                                                      aaGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTg                                                      >  1:473831/1‑141 (MQ=255)
                                                                                        gCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGat                                                    <  1:132292/141‑1 (MQ=255)
                                                                                        gCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGat                                                    >  2:300834/1‑141 (MQ=255)
                                                                                        gCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGat                                                    <  1:589639/141‑1 (MQ=255)
                                                                                         cTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGatt                                                   <  1:474472/141‑1 (MQ=255)
                                                                                          tttGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGattc                                                  >  1:209561/1‑141 (MQ=255)
                                                                                          tttGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAAATAGCTACACAGCAAACTTGTCTCattc                                                  <  2:758019/141‑1 (MQ=255)
                                                                                                   agatgaACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAATAAGTGGATCTTAAGAAAGACAGCTACACAGCAAACTTGTTTGATTCAAATTCAAt                                         <  1:594714/140‑1 (MQ=255)
                                                                                                             ggTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAATCAAAtt                               >  1:243566/1‑141 (MQ=255)
                                                                                                             ggTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAAtt                               >  1:122973/1‑141 (MQ=255)
                                                                                                             ggTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAAtt                               >  1:419983/1‑141 (MQ=255)
                                                                                                             ggTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAAtt                               >  2:443594/1‑141 (MQ=255)
                                                                                                             ggTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAAtt                               <  1:285615/141‑1 (MQ=255)
                                                                                                             ggTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAAtt                               >  1:558812/1‑141 (MQ=255)
                                                                                                               tgtAACACCTGAGAAGATTCGAAATGCAGAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGTTAATAGTAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGa                             <  1:717545/141‑1 (MQ=255)
                                                                                                               tgtAACACCTGAGAAGATTCGAAATGCAGAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGa                             <  1:287359/141‑1 (MQ=255)
                                                                                                               tgtAACACCTGAGAAGATTCGAAATGCAAAATCTTCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGa                             <  2:485254/141‑1 (MQ=255)
                                                                                                               tgtAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGa                             <  2:379289/141‑1 (MQ=255)
                                                                                                               tgtAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGa                             <  2:765557/141‑1 (MQ=255)
                                                                                                               tgtAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGa                             <  1:710390/141‑1 (MQ=255)
                                                                                                                 tAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGata                           >  2:433323/1‑141 (MQ=255)
                                                                                                                 tAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTAGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGata                           >  2:564521/1‑141 (MQ=255)
                                                                                                                   acacCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGATAtt                         <  1:300834/141‑1 (MQ=255)
                                                                                                                    aacCCGAGAAGATTCGAAATGCAAAATCTGCCGGTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGATATTa                        <  2:195854/140‑1 (MQ=255)
                                                                                                                               aTTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGATATTAATGTAAAAACt             <  1:582347/141‑1 (MQ=255)
                                                                                                                                ttCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGATATTAATGTAAAAACta            >  2:356308/1‑141 (MQ=255)
                                                                                                                                      aTGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGATATTAATGTAAAAACtaaataaa      >  2:604933/1‑141 (MQ=255)
                                                                                                                                        gCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGATATTAATGTAAAAACTAAATAAaca    >  2:597049/1‑141 (MQ=255)
                                                                                                                                        gCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGATATTAATGTAAAAACTAAATAAaca    >  2:597048/1‑141 (MQ=255)
                                                                                                                                          aaaaTCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGATATTAATGTAAAAACTAAATAAACACt  <  2:515561/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
AGATGGCATTATGGCTTTCAATCAATCTTACTTAAATTTACCAAAAGATAGCCAATTATCATATTTAGATTTAGGAAATAAAGTTAAAGCTTTGTTATATGATGAACGCGGTGTAACACCTGAGAAGATTCGAAATGCAAAATCTGCCGTTTACACGATTACTTGGAAAGATGGTAGTAAAAAAGAAGTGGATCTTAAGAAAGATAGCTACACAGCAAACTTGTTTGATTCAAATTCAATTAAACAAATTGATATTAATGTAAAAACTAAATAAACACT  >  USA300TCH1516_ALE/984220‑984498

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