Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 111,512 A→G 40.3% S194S (TCA→TCG USA300HOU_RS00520 → putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE111,5120AG40.3% 46.6 / 74.3 75S194S (TCA→TCGUSA300HOU_RS00520putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (26/18);  new base G (10/20);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 3.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

TACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAA  >  USA300TCH1516_ALE/111374‑111650
                                                                                                                                          |                                                                                                                                          
tACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGTTTTTTTTAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATTATATAATCCAGAGGTGCCGATTTATTCAgt                                                                                                                                          >  2:398198/1‑140 (MQ=21)
 aCTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCTCTATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCTAGAGTTGCCGAATTATTCAGCt                                                                                                                                         >  2:586106/1‑141 (MQ=25)
 aCTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATTAAGGATAAGATATATCATATAATCCAGAGGTGCCGAGTTATTCATCt                                                                                                                                         >  2:159350/1‑141 (MQ=25)
 aCTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTCTAAAGATTTAATGAATTATTAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCt                                                                                                                                         >  2:108192/1‑141 (MQ=255)
  cgtaGGAAGTGAAAGAATAAAATATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGAAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCta                                                                                                                                        <  2:276886/138‑2 (MQ=18)
  atGAGGAATTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCta                                                                                                                                        <  2:402795/140‑2 (MQ=255)
   tGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCtaa                                                                                                                                       <  1:679637/141‑3 (MQ=255)
    gAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATGCAGAGGTGCCGAGTTATTCGGCTaaa                                                                                                                                      >  1:118749/1‑141 (MQ=255)
     aGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                                                                                                     <  1:591133/141‑1 (MQ=255)
      ggAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATTAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                                                                    <  2:29001/141‑1 (MQ=255)
      ggAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                                                                    <  2:712376/141‑1 (MQ=255)
      ggAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAata                                                                                                                                    <  2:351944/141‑1 (MQ=255)
              aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                                                            >  2:659023/1‑141 (MQ=255)
              aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTaa                                                                                                                            >  1:187572/1‑141 (MQ=255)
               agaagaAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAAc                                                                                                                           >  1:560552/1‑141 (MQ=255)
                gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                                                          <  2:720150/141‑1 (MQ=255)
                gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACt                                                                                                                          >  2:521695/1‑141 (MQ=255)
                gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTACCGAGTTATTCACCCCAATATCCATTAACt                                                                                                                          >  2:682167/1‑141 (MQ=255)
                   aaaaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                       <  1:159350/140‑1 (MQ=255)
                    aaaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATAAATTAACTAatt                                                                                                                      >  2:528309/1‑140 (MQ=255)
                     aaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                                                                     <  1:314505/141‑1 (MQ=255)
                        aaTAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAAtgatga                                                                                                                  <  2:620454/140‑1 (MQ=255)
                         aaaaaaaaaaaaATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                                                                 <  1:316442/139‑1 (MQ=255)
                         aTAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                                                                 >  1:405964/1‑141 (MQ=255)
                            aaaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                                                              >  1:646319/1‑141 (MQ=255)
                            aaaaaaGAAATCGAAAACGTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGGAGGAGAGAGATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                                                              <  1:357430/141‑1 (MQ=17)
                                   aaaTCGAAAACTTTAAGTTTTTTGCGAAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                       <  2:681755/141‑1 (MQ=255)
                                    aaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCCGCTCAATATCAATTAACTAATGATGATTCTAATGTaca                                                                                                      >  1:3256/1‑139 (MQ=255)
                                     aTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                                     >  1:29307/1‑141 (MQ=255)
                                     aTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                                     >  2:293899/1‑141 (MQ=255)
                                     aTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                                     >  2:207321/1‑141 (MQ=255)
                                         aaaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                 <  1:160276/141‑1 (MQ=255)
                                         aaaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                 <  1:648518/141‑1 (MQ=255)
                                              tttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGTGGTATGACAAACAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGt                                                                                            <  1:354718/141‑1 (MQ=18)
                                              tttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGt                                                                                            <  1:161230/141‑1 (MQ=255)
                                                tAAGTTTTTTGAGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGCGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaa                                                                                          <  1:534017/141‑1 (MQ=255)
                                                 aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                                                         >  1:216435/1‑141 (MQ=255)
                                                 aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCACTATCAATTAACTAATGCTGATTCTAATGTCAACCAAttcactcaa                                                                                         >  2:7900/1‑134 (MQ=21)
                                                  aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTACAGATAATTGTTATGTAAAACAATTACGCaaaa                                                                                        >  1:384685/1‑141 (MQ=18)
                                                  aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaaaa                                                                                        >  1:498819/1‑141 (MQ=255)
                                                  aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaaaa                                                                                        <  1:126931/141‑1 (MQ=255)
                                                  aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaaaa                                                                                        <  1:224108/141‑1 (MQ=255)
                                                  aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCAAGTTATTCAGATCAATATTGACTAACGAATGCTAATTATAACGTAAAACAATTAAGTaaaa                                                                                        >  2:729454/1‑141 (MQ=25)
                                                     tttttGCGCAATATAGCAACTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGat                                                                                     <  2:118749/141‑1 (MQ=255)
                                                        ttGCGCAATACAGTAATTTTAAAGATTTAGTGAATTATAAAGAGGGAGATATATCATATAATCCAGAGGTGCCGAGTGATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATg                                                                                  <  2:11588/141‑1 (MQ=14)
                                                         tGCGTAATATAGTAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGa                                                                                 <  1:254184/141‑1 (MQ=14)
                                                           cgaaaTATAGCAATTTTAAAGATTTAATGAATTATAAAGATAGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGata                                                                               <  2:505045/138‑3 (MQ=14)
                                                             cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATCTGAtata                                                                             >  1:599723/1‑141 (MQ=14)
                                                                 ataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAccga                                                                         <  1:80919/141‑3 (MQ=11)
                                                                   aGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAccgacg                                                                       <  1:126765/141‑5 (MQ=11)
                                                                   aGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAccgacg                                                                       <  2:481588/141‑5 (MQ=11)
                                                                    gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATACAGAGGTGTCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacga                                                                      >  2:522608/1‑139 (MQ=11)
                                                                         tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATaaa                                                                 >  1:314381/1‑141 (MQ=255)
                                                                         tttAAAGATATAATGAAATATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATaaa                                                                 <  1:412482/141‑1 (MQ=255)
                                                                           tAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAga                                                               >  1:132875/1‑140 (MQ=255)
                                                                              aGATTTAATGAATTATAAAGATGGAGATATATCTTATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGcc                                                            <  2:358221/141‑1 (MQ=255)
                                                                                     aggAATGATAAAGATGGAGATAGATCATATAATCAAGAGGTGCCTAGTTATTCAGATCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtg                                                     <  2:245357/139‑1 (MQ=17)
                                                                                       gAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgtt                                                   >  1:488795/1‑141 (MQ=255)
                                                                                         aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGTCAAGTTAGGCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttgt                                                 <  2:532495/141‑1 (MQ=255)
                                                                                         aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttgt                                                 <  2:169322/141‑1 (MQ=255)
                                                                                                  aTGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATAAGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAggaa                                        >  2:260009/1‑139 (MQ=255)
                                                                                                       gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAggg                                   <  2:479877/141‑1 (MQ=255)
                                                                                                        atatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGa                                  >  1:644435/1‑141 (MQ=255)
                                                                                                            atCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAttt                              <  2:498819/141‑1 (MQ=255)
                                                                                                            atCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAttt                              >  1:233536/1‑141 (MQ=255)
                                                                                                                   aaTACAGAGGTGACGAGTTATTCAGATCAATATCAATGAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                       <  2:713989/141‑1 (MQ=255)
                                                                                                                    aTCCCGAGTTGCCGAGTTATTCAGCTTAATAGCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACGAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGtt                      <  1:154619/141‑1 (MQ=14)
                                                                                                                      ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                    >  1:460816/1‑141 (MQ=255)
                                                                                                                      ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAAGGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                    >  1:716937/1‑141 (MQ=255)
                                                                                                                        agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatc                  >  1:118527/1‑141 (MQ=255)
                                                                                                                            gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGtt              >  1:642521/1‑141 (MQ=255)
                                                                                                                             tGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGCAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTg             >  1:54335/1‑141 (MQ=255)
                                                                                                                                     tATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGCATTTCAAAGGTTCATCAGTTggttctata     >  1:185559/1‑134 (MQ=255)
                                                                                                                                     tATTCAACTCAATATCAATTAACCAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAACTAAAGCGCCGAAGCCGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaa     <  1:640823/141‑1 (MQ=255)
                                                                                                                                        tCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaaa  >  1:148330/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAA  >  USA300TCH1516_ALE/111374‑111650

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: