Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 111,516 C→A 31.6% Q196K (CAA→AAA)  USA300HOU_RS00520 → putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE111,5160CA31.6% 115.0 / 53.4 84Q196K (CAA→AAA) USA300HOU_RS00520putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base A (16/11);  total (47/37)
Fisher's exact test for biased strand distribution p-value = 8.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

TGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAA  >  USA300TCH1516_ALE/111377‑111651
                                                                                                                                           |                                                                                                                                       
tGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTGTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTc                                                                                                                                         <  1:742586/140‑1 (MQ=255)
  aGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAAATATAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTaaa                                                                                                                                       >  2:270182/1‑140 (MQ=255)
    gAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAata                                                                                                                                     <  1:352854/140‑1 (MQ=255)
     aaGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                                                                                                    <  2:305047/140‑1 (MQ=255)
      aGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGAGCAATATAGAAATTATAAAGATTTAATGAATTATAGAGATGGAGATATATCATATAATCAAGAGGTGCCGAGTTAGTCAGCTCCataac                                                                                                                                   >  2:307205/1‑138 (MQ=18)
       gTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCa                                                                                                                                  <  1:575409/140‑1 (MQ=255)
         tAAAGAAGAAAACATAAAAAAAGAAATCGCAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTCTCTATCAAt                                                                                                                                >  1:528196/1‑140 (MQ=255)
         tAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGCGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                                                                >  2:524038/1‑140 (MQ=255)
         tAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                                                                >  1:545558/1‑140 (MQ=255)
         tAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                                                                >  1:803311/1‑140 (MQ=255)
         tAAAGAAGAAAACATAAAAAAAGAAACCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                                                                >  2:607838/1‑140 (MQ=255)
          aaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATGATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                                               >  2:150698/1‑140 (MQ=255)
          aaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAtt                                                                                                                               >  1:197996/1‑140 (MQ=255)
          aaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTCGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAtt                                                                                                                               >  1:359879/1‑140 (MQ=255)
           aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATAGCAGTTa                                                                                                                              >  1:394689/1‑140 (MQ=255)
           aagaagAAAACATAAAAAAAAAAATCGAAAACCTTAAGCTTCATGCGAAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCAAGAGGTGCCGAGTTATTCAGCTCAATATCAATTa                                                                                                                              >  1:171075/1‑140 (MQ=255)
            agaagaAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                                                             <  1:271835/140‑1 (MQ=255)
            agaagaAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTaa                                                                                                                             <  1:718796/140‑1 (MQ=255)
            agaagaAAAAATAAAAAAAGAAATCGAAAACGTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                                                             <  2:756773/140‑1 (MQ=18)
             gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAAc                                                                                                                            <  1:688976/140‑1 (MQ=255)
              aagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGGGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACt                                                                                                                           <  2:64335/140‑1 (MQ=255)
                 aagaaaTAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                        <  2:276786/135‑1 (MQ=255)
                 aaaaCATAACAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                        >  1:638373/1‑140 (MQ=255)
                    aCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                                                     >  2:278561/1‑140 (MQ=255)
                    aCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                                                     >  1:791429/1‑140 (MQ=255)
                    aCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgat                                                                                                                     >  1:658983/1‑140 (MQ=255)
                     cATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGATAAATATCAATTAACTAatgatg                                                                                                                    >  2:190667/1‑140 (MQ=255)
                      aTAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAAtgatga                                                                                                                   <  2:337464/140‑1 (MQ=255)
                       tAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATgatgat                                                                                                                  <  1:8499/140‑1 (MQ=255)
                          aaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtat                                                                                                               <  1:434183/140‑1 (MQ=255)
                            aaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATaa                                                                                                             <  2:453468/140‑1 (MQ=255)
                            aaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATaa                                                                                                             <  1:163959/140‑1 (MQ=255)
                               gAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGt                                                                                                          <  1:435256/140‑1 (MQ=255)
                               gAAACCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGt                                                                                                          >  1:379905/1‑140 (MQ=255)
                                 aaTCGAAATCTTCAAGTTTTTTGCTCAATAGAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAAAATCCATTAACTATTGATGATTATAATGTaa                                                                                                        <  2:289965/140‑1 (MQ=17)
                                 aaTCGAAAACTTTAAGTTTTTTGCGTAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTATCTAATGATGATTATAATGTaa                                                                                                        <  2:491640/140‑1 (MQ=255)
                                 aaTCGAAAACTTTAAGTTTGTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                        <  1:740711/140‑1 (MQ=255)
                                   tCGAAAACTTTAAGTTTTTTGCGCGATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                                      <  2:600048/140‑1 (MQ=255)
                                   tCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                                      >  1:438873/1‑140 (MQ=255)
                                    cGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGCTTATAATGTAAAAc                                                                                                     >  2:43460/1‑140 (MQ=255)
                                       aacattaaaGTTTTTAACCCAAGAAATCAATTTTAAAGATTAAATGACTTATAAAGATGTAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                  <  1:298632/133‑1 (MQ=17)
                                            ttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGt                                                                                             <  1:649673/140‑1 (MQ=255)
                                            ttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGt                                                                                             <  1:649190/140‑1 (MQ=255)
                                             tAAGTTTTTTCCACAATATATCAATTTTACGCATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACgca                                                                                            <  2:164165/140‑3 (MQ=14)
                                              aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaa                                                                                           <  1:728402/140‑1 (MQ=255)
                                               aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGGGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAACGATGATTATAATGTAAAACAATTACGCaaa                                                                                          >  1:349147/1‑140 (MQ=18)
                                                gTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaaa                                                                                         <  2:752481/140‑1 (MQ=255)
                                                    tttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGata                                                                                     >  1:710363/1‑140 (MQ=255)
                                                        cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATAATATAATCCAGAGGTGCCGAGTTATTGGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAACCatatact                                                                                 >  2:283005/1‑137 (MQ=14)
                                                         gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGata                                                                                >  2:274537/1‑138 (MQ=255)
                                                            agatcgCAATTGTAAAGATTTAATGAATGATAAAGATGGAGATATATGATATAATCTAGAAGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATAc                                                                             <  1:183714/135‑1 (MQ=18)
                                                              ataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACca                                                                           <  2:587739/140‑1 (MQ=255)
                                                                 gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaaca                                                                        <  2:576484/140‑1 (MQ=255)
                                                                 gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaaca                                                                        >  2:621180/1‑140 (MQ=255)
                                                                  cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaacaa                                                                       <  2:578929/140‑1 (MQ=255)
                                                                   aaTTTTAAAGATTGAATGAATTATAAAGCTGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTATAATTTATGACATACCAACaaa                                                                      <  2:98670/140‑1 (MQ=18)
                                                                    aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacgagt                                                                     >  1:280412/1‑136 (MQ=11)
                                                                        tAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAg                                                                 >  1:9592/1‑140 (MQ=255)
                                                                        tAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAg                                                                 >  1:246960/1‑140 (MQ=255)
                                                                         aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGACGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTCATACAATTACGTAAGAGATATGACATACGAAGAAATAAAgc                                                                >  1:45455/1‑140 (MQ=18)
                                                                          aaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcg                                                               >  1:803817/1‑140 (MQ=255)
                                                                               ttAATGAATTATAAAGATGGAGATATATCATATAATCAAGAGGTGCAGAGTTATAAAGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGaa                                                          >  1:656218/1‑140 (MQ=17)
                                                                                tAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAg                                                         <  2:684171/140‑1 (MQ=255)
                                                                                      aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCAGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttg                                                   <  1:547281/140‑1 (MQ=255)
                                                                                             aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAg                                            >  1:683603/1‑140 (MQ=255)
                                                                                             aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAg                                            >  1:682435/1‑140 (MQ=255)
                                                                                                 ggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGGAAGGTAc                                        >  2:234236/1‑140 (MQ=255)
                                                                                                      tatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGCGGAATATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGa                                   >  1:507998/1‑140 (MQ=255)
                                                                                                          tCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTCGTTGAAAGGTACAGGGAAttt                               >  2:203020/1‑140 (MQ=255)
                                                                                                            atatAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTaa                             >  1:388592/1‑140 (MQ=255)
                                                                                                              ataATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACACAACAGATATGACATCCTATGAGGGAAGTCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTaaaa                           <  2:572550/140‑1 (MQ=14)
                                                                                                                  tCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGtt                       >  1:642582/1‑140 (MQ=255)
                                                                                                                  tCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGAAATTTAAAAGGtt                       >  2:694303/1‑140 (MQ=255)
                                                                                                                   ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtc                      >  1:109064/1‑140 (MQ=255)
                                                                                                                        gatGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCACCGAAAATGTTATTGAAAAATACAGAGAATATAAAAAATtcataaa                 >  2:501577/3‑137 (MQ=14)
                                                                                                                         gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGt                >  1:469074/1‑140 (MQ=255)
                                                                                                                         gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGt                >  1:349818/1‑140 (MQ=255)
                                                                                                                           gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTg              >  1:156874/1‑140 (MQ=255)
                                                                                                                            ccGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCCATAACGGGCGGGGTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTgg             >  2:778408/1‑140 (MQ=14)
                                                                                                                              gAGTTATTCAGCTCAACATCAATTAACCAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGat           <  1:316299/140‑1 (MQ=255)
                                                                                                                                   aTTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaa      <  1:400394/140‑1 (MQ=255)
                                                                                                                                    ttCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCAGTTGTTGAAAGGTCCAGCGTAGTAAAAAGGTTCATCAGTTGGATATaaaa     >  1:715293/1‑140 (MQ=17)
                                                                                                                                     tCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACACTCACGTAAGGGAGTTGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGATACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaa    <  1:182331/140‑1 (MQ=17)
                                                                                                                                     tCAACTCAATATCAATTAACGAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaa    <  2:586125/140‑1 (MQ=255)
                                                                                                                                      cAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAATGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaaa   >  1:651026/1‑140 (MQ=255)
                                                                                                                                       aGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaaaa  >  1:411572/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                       
TGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAA  >  USA300TCH1516_ALE/111377‑111651

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