Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 111,516 C→A 28.9% 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,5160CA28.9% 135.0 / 58.5 90Q196K (CAA→AAA) USA300HOU_RS00520putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (28/36);  new base A (10/16);  total (38/52)
Fisher's exact test for biased strand distribution p-value = 8.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

TGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAA  >  USA300TCH1516_ALE/111377‑111648
                                                                                                                                           |                                                                                                                                    
tGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGGGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCa                                                                                                                                     <  2:159541/141‑1 (MQ=255)
 gAGGCAGTTAAAGAAGAAAACATAAAAAAAGACATCGAAAACTTTACGTTTTGTGCGCAATATAGCAATTTTAAAGATTTAATGAGTTATAAAGATGGAGCTATATCATATAATCCAGAGGTGCCGCGTTATTCGGCTaaa                                                                                                                                    >  2:196822/1‑141 (MQ=18)
 gAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTCTGCGCAATATAGCAACTTTAAAGATTTAATGAATTATAAAGATGCAGATATATCATATCATCCAAAGGTCACGAGTTATTCACCTcca                                                                                                                                    >  2:239525/1‑139 (MQ=21)
  aGGAAGTTAAAGAAGCAAACATAAAAAAAGAAAACGATGACTTTAACTTGTTTGCGCAATATAGCAATTTTAAAGATTTAATGATTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAat                                                                                                                                   <  1:144058/141‑1 (MQ=25)
  aGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAat                                                                                                                                   >  2:303711/1‑141 (MQ=255)
        ttAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAt                                                                                                                             >  1:402807/1‑141 (MQ=255)
          aaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                                                                           <  1:21080/141‑1 (MQ=255)
          aaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTa                                                                                                                           >  1:453307/1‑141 (MQ=255)
          aaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCAGAGATATTCAGCTCAATATCAATTa                                                                                                                           >  1:16436/1‑141 (MQ=255)
          aaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGACAAGTTATTCGGCTAAATATAAAGTa                                                                                                                           >  2:547401/1‑141 (MQ=21)
           aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                                                          <  2:648046/141‑1 (MQ=255)
             aaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                                                        <  1:153914/140‑1 (MQ=255)
             aaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACt                                                                                                                        <  2:448038/140‑1 (MQ=255)
              aatacaaCATAAAAAAAGAAATCTAAAACTTTAAGTTTTTTGCCCAATATATCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTa                                                                                                                       <  2:674200/136‑1 (MQ=255)
              aagaaAACATAAACAAAGAAATCGAAATCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGCTGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTa                                                                                                                       <  2:29414/141‑1 (MQ=255)
              aagaaAACATAAAAAAATAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGCGCTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTa                                                                                                                       <  2:107327/141‑1 (MQ=255)
                 aaaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGGATTAGAAATATGGAGATATATCCTATAATCCAGAGGTGCCGAGGTATTCGGCTAAATATCAATTAACTAatg                                                                                                                    >  1:32986/1‑141 (MQ=17)
                 aaaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatg                                                                                                                    >  1:128853/1‑141 (MQ=255)
                 aaaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatg                                                                                                                    <  1:404205/141‑1 (MQ=255)
                 aaaaCATAAAAAAAGAAATCGAAAACTGTAAGTTTTTTGCGCAATATAGCAAGTTTAAAGATTTAATGAATTATAAAGATGGAGATATAGCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatg                                                                                                                    >  1:148818/1‑141 (MQ=17)
                   aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                                                  <  1:194864/141‑1 (MQ=255)
                   aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                                                  >  2:181803/1‑141 (MQ=255)
                   aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCAAGAGGTGCCGAGTTATTCGGCTCAATATCAATTAACTAataaa                                                                                                                  >  2:296422/1‑138 (MQ=255)
                   aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGAAATTTTAAAGATTTAATGAATTATAAAGATGGTGGTCGCGTGAACGATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                                                  >  1:193719/1‑141 (MQ=17)
                     cATCAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                                                                <  1:638338/141‑1 (MQ=255)
                     cATAAAAAAAGAAATCGAAAACTTGAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                                                                <  1:638344/141‑1 (MQ=255)
                     aaTAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                                                                <  1:294794/140‑1 (MQ=255)
                      aaaaaaaaaGAAATCGAAAACTTTAAGCTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATgatgat                                                                                                               <  2:250543/139‑1 (MQ=255)
                      aTAAAAAGATAAAGCGAAAACTTTAAGTTTTTGACGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATgatgat                                                                                                               <  2:663374/141‑1 (MQ=21)
                      aTAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                                                               >  2:269163/1‑141 (MQ=255)
                       tAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATGATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTCACTAatgagaatt                                                                                                              >  1:71893/1‑136 (MQ=255)
                         aaaaaaaaaaTCGAAAACTTTAAGTTTTTTGCACAATATAGTAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAACCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtat                                                                                                            <  1:81583/141‑1 (MQ=255)
                          aaaaaGAAATAGAAAACTTTAAGTTTTTGGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                                                                                                           <  2:172771/141‑1 (MQ=255)
                            aaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAAt                                                                                                         <  2:56329/141‑1 (MQ=255)
                            aaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAAt                                                                                                         >  1:269807/1‑141 (MQ=255)
                             aaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATg                                                                                                        <  1:416624/141‑1 (MQ=255)
                              aGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGt                                                                                                       <  1:560716/141‑1 (MQ=255)
                              aGAAATCGAAAAATTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGt                                                                                                       <  1:239344/141‑1 (MQ=255)
                               gAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTa                                                                                                      <  2:300591/141‑1 (MQ=255)
                               gAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTa                                                                                                      <  2:56135/141‑1 (MQ=255)
                                aaaTCTAAAAATTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGAAGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                     <  1:524303/141‑1 (MQ=17)
                                aaaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                     <  2:276124/141‑1 (MQ=255)
                                aaaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAAAATTAAATCGTGAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                     >  1:306678/1‑141 (MQ=17)
                                aaaCCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGGCGATTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                     <  2:322826/141‑1 (MQ=255)
                                 aaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                                    <  2:186367/141‑1 (MQ=255)
                                 aaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                                    <  2:178940/141‑1 (MQ=255)
                                 aaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                                    >  1:641360/1‑141 (MQ=255)
                                   tCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAAc                                                                                                  >  1:464316/1‑141 (MQ=255)
                                      aaaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                               <  2:128853/141‑1 (MQ=255)
                                        aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTa                                                                                             >  1:145693/1‑141 (MQ=255)
                                         aCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                                                                            <  1:303711/141‑1 (MQ=255)
                                          cTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACg                                                                                           <  2:301773/141‑1 (MQ=255)
                                          cTTTAAGTTGTTTGCGCAATCTAGCAATTTTAAATATTTAATGAATTATACCGACGGAGATATCTCATATAATCCATATGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACg                                                                                           <  2:113365/141‑1 (MQ=14)
                                           tttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACgc                                                                                          <  1:540559/141‑2 (MQ=255)
                                           ttgaaGTTTTGAGCGCAAACAAGCAATTTTAAAGATTTAACGAATTATAAAGATGGAGATATATCATATAACCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACgc                                                                                          <  2:544596/138‑2 (MQ=14)
                                                     ttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATg                                                                                <  1:219731/141‑1 (MQ=255)
                                                     ttGCGCAATATAGCAATTTTAAAGATTTAATGAAATATAAAAATGGAAATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATg                                                                                <  1:6284/141‑1 (MQ=255)
                                                      tacaaaaTAAAGCAATTTTAAATATTTAATAAAATATAAAGAGGGAGACACATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAAATAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGa                                                                               <  1:93472/136‑1 (MQ=17)
                                                      tGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGa                                                                               <  1:22734/141‑1 (MQ=255)
                                                        cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTATAACAATTACGTAAAAGATATGACa                                                                             >  1:439343/1‑141 (MQ=255)
                                                        cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGTTATGACa                                                                             >  2:661394/1‑141 (MQ=255)
                                                         gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAGTATCAATTAACTAATGATGATTATAATGTAAAACAATTAAGTAAAAGATATGACAt                                                                            >  2:404715/1‑141 (MQ=255)
                                                          cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAAGGAATGAGGATTATAATGTAAAACAATTACGTAAAAGATATGACATa                                                                           >  2:72392/1‑141 (MQ=255)
                                                            atatAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAcc                                                                         >  1:341366/1‑141 (MQ=255)
                                                               taGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaac                                                                      <  1:323065/141‑1 (MQ=255)
                                                                aGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaaca                                                                     <  2:230325/141‑1 (MQ=255)
                                                                 gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaacaa                                                                    <  2:269807/141‑1 (MQ=255)
                                                                  cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacgag                                                                   >  2:194192/1‑138 (MQ=11)
                                                                         aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATCTGACATACACACAAATAAAgcg                                                            >  2:543104/1‑141 (MQ=255)
                                                                          aaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                                                           >  2:70677/1‑141 (MQ=255)
                                                                          aaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                                                           >  2:530474/1‑141 (MQ=255)
                                                                                tAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGc                                                     <  1:20562/141‑1 (MQ=255)
                                                                                aaaTGAATTATAAAGATGGAGATCTATCTTATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGc                                                     <  1:113266/140‑1 (MQ=255)
                                                                                      aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttgt                                               <  1:271950/141‑1 (MQ=255)
                                                                                               aTGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAAACCAAGAGAAGTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTa                                      <  1:160819/141‑1 (MQ=17)
                                                                                                        tatCATATAATCCAGAGGTGCCGCGCTCTCCACCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTCCGTACAAGATATGACATACCAACAACTAAAACGCCGAAGCTGTCGTTGAAAGGTACAGGGAAtt                             >  2:634891/1‑141 (MQ=17)
                                                                                                        tatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAGGGTACAGGGAAtt                             >  1:278598/1‑141 (MQ=255)
                                                                                                                   ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                  >  1:196551/1‑141 (MQ=255)
                                                                                                                    cAGAGGTGCCGAGTTATTCAGATCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTacat                 >  1:244722/1‑141 (MQ=255)
                                                                                                                     agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatc                >  2:548858/1‑141 (MQ=255)
                                                                                                                      gagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatca               <  2:644701/141‑1 (MQ=255)
                                                                                                                      gagGTGCCGAGTTATTCAGCACAATATCAATTAACTAAAGATGATTAAAATGTAAAACAATTACGTACAAGATATGACATAGCAGCAAATGAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatca               <  2:111511/141‑1 (MQ=14)
                                                                                                                       agGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAg              >  2:27320/1‑141 (MQ=255)
                                                                                                                          tGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTg           >  1:584880/1‑141 (MQ=255)
                                                                                                                          tGCCGAGTTATTCAGCTCAAGATCAATTAACTAATGATGATTATAATGTAAAACAATTTCGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTAGAGGGAATTTACAAGGTTCATCAGTTg           >  1:421897/1‑141 (MQ=255)
                                                                                                                           gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTgg          >  2:498653/1‑141 (MQ=255)
                                                                                                                            ccGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGa         >  1:276233/1‑141 (MQ=255)
                                                                                                                              gAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGGTCATCAGTTGGata       >  2:10244/1‑141 (MQ=255)
                                                                                                                               aGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTAAGatat      <  2:673190/141‑1 (MQ=255)
                                                                                                                                   aTTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaa  <  1:658363/141‑1 (MQ=255)
                                                                                                                                   aTTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaa  <  1:245920/141‑1 (MQ=255)
                                                                                                                                   aTTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaa  <  1:524446/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                    
TGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAA  >  USA300TCH1516_ALE/111377‑111648

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: