Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 110,651 C→T 10.2% D170D (GAC→GAT)  USA300HOU_RS00515 → putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE110,6510CT10.2% 197.5 / 11.0 78D170D (GAC→GAT) USA300HOU_RS00515putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (34/36);  new base T (4/4);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

AGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGT  >  USA300TCH1516_ALE/110514‑110789
                                                                                                                                         |                                                                                                                                          
aGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGga                                                                                                                                         >  2:749230/1‑141 (MQ=255)
aGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTAAATATGGAGAATTTAAAAATATGAAAAATTATAGAGTCGga                                                                                                                                         >  2:280838/1‑141 (MQ=17)
  ttatCCGGTTATGTTGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCTAAAACTTTAAGTTCTTTGTTCAATAGGGCGACTTTAAAAATTTGAAAAATTATAAAGACGgaga                                                                                                                                       <  1:166445/140‑1 (MQ=17)
  atatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAATAATTATAAAGACGgaga                                                                                                                                       >  2:620896/1‑141 (MQ=255)
          ttAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCat                                                                                                                               >  1:375507/1‑141 (MQ=255)
           tAAAATGGTAGATAATAAAATCATTACAACAAAAGAAATAAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCata                                                                                                                              <  1:301543/141‑1 (MQ=255)
             aaaTGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCAtata                                                                                                                            >  2:48730/1‑141 (MQ=255)
                tGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATc                                                                                                                         >  1:219383/1‑141 (MQ=255)
                   tcaataataAAATCATTCCAACAAAATAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                      <  1:681502/138‑1 (MQ=255)
                   aagATAATAAAATCATTACAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                      <  1:277463/139‑1 (MQ=255)
                      ataataAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAgg                                                                                                                   <  1:147807/141‑1 (MQ=255)
                       taataaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATctaaagtt                                                                                                                  >  2:728062/1‑134 (MQ=255)
                        aataaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTg                                                                                                                 >  1:505193/1‑141 (MQ=255)
                          taaAATCAATCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                                                               <  1:426212/141‑1 (MQ=255)
                               tCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGtt                                                                                                          >  2:640855/1‑141 (MQ=255)
                               tCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGtt                                                                                                          >  1:301922/1‑141 (MQ=255)
                                  ttCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAACTTTGAAAAATTATAAAGACGAAGATATTACATATAATACAAAAGTACCGAGTTAtt                                                                                                       >  1:192670/1‑141 (MQ=17)
                                     aaacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGCAAATTAAATCTGTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                                                                    <  1:710790/140‑1 (MQ=18)
                                             gAAATTAAAGATGAAATAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGTGACTTTAAAAATTTGAAAAACTAGAAAGACTGAGATATTTCATATAATCCAGATGTGTCTAGTTACTAGGCTAAataa                                                                                            >  2:328527/1‑140 (MQ=14)
                                               aaTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTAAATATTCAAAAAGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCAAATAATAAAGAGGTGGGAAAGTATGCGGGTAAatataa                                                                                          >  2:636885/1‑139 (MQ=17)
                                                 ttAAAGATGAAAAAATACAAAAAGAAATCGAAATCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                        <  2:2769/141‑1 (MQ=255)
                                                    aaaaagAAAAAAAAAAAAAATGAATCGAAAACTTTAAGTTCTTTGTTCAAAATGGGAACTTTAAAAATTTAAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                     <  2:458402/136‑1 (MQ=14)
                                                    aaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                     <  1:275709/141‑1 (MQ=255)
                                                    aaGATGAAAAAATAAAAAAAGAAATAGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                     <  1:61785/141‑1 (MQ=255)
                                                    aaGACGAAAAAATAAAAAAAGAAATCGAAAACTCTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                     <  1:183868/141‑1 (MQ=255)
                                                       aggAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                  <  2:337364/139‑1 (MQ=255)
                                                           aaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                              >  2:669031/1‑141 (MQ=255)
                                                            aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                             <  2:392013/141‑1 (MQ=255)
                                                             aaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                                                            >  1:642143/1‑141 (MQ=255)
                                                               aTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                          >  1:100031/1‑141 (MQ=255)
                                                                  aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGCTATTTCACAAAGCTCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                       >  1:589109/1‑141 (MQ=17)
                                                                    aaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTAGGCTAAATATCAATTAACTAATGATGATtatgc                                                                     >  1:41057/1‑139 (MQ=255)
                                                                       aGAAATCAAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATGTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGt                                                                  <  2:199090/141‑1 (MQ=255)
                                                                           aTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                              >  1:154594/1‑141 (MQ=255)
                                                                           aTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                              >  2:150672/1‑141 (MQ=255)
                                                                             cGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAcaca                                                            >  1:545785/1‑141 (MQ=255)
                                                                               aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                          <  1:352604/141‑1 (MQ=255)
                                                                                aaaCTGTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                         <  2:713680/141‑1 (MQ=255)
                                                                                 aaCTTCAAGGTCTGTTTTCAACATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTa                                                        <  1:560925/141‑1 (MQ=17)
                                                                                  aCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                                       <  2:512104/141‑1 (MQ=255)
                                                                                   cTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACg                                                      <  1:48730/141‑1 (MQ=255)
                                                                                    tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                     <  1:612143/141‑1 (MQ=255)
                                                                                      tAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa                                                   <  1:577223/141‑1 (MQ=255)
                                                                                      tAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAAACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa                                                   <  1:652168/141‑1 (MQ=255)
                                                                                          ttCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAga                                               >  2:334045/1‑141 (MQ=255)
                                                                                             tttGTTCAATATGGCGACTTTAACAATTTGAAAAATTATAAAGACGGAGATCTTTCATATAATCCCGAGGTGCCGAGTTATTCGGCTAAATATCTATTCGCTTATGCTGTGTTTTTTGTAAAACACTTACGCACAAcatct                                            >  1:185143/1‑139 (MQ=17)
                                                                                                  tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgatat                                       >  1:628937/1‑141 (MQ=255)
                                                                                                       atGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTGATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcga                                  <  2:179828/141‑1 (MQ=255)
                                                                                                       atGGCGACTTTAAAAATTTGAAAAATTAGAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcga                                  <  2:530576/141‑1 (MQ=255)
                                                                                                        tGGTGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgac                                 <  1:78560/141‑1 (MQ=255)
                                                                                                         ggCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATCTAATACAAGTACGCAAAAGATATGATATACcgacc                                >  2:690304/1‑140 (MQ=255)
                                                                                                          gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACggacga                               >  1:236354/1‑141 (MQ=255)
                                                                                                             aCTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTa                            >  2:406487/1‑141 (MQ=255)
                                                                                                              ttttAAAGATTTAATGAAGTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                           <  2:630671/140‑1 (MQ=11)
                                                                                                              cTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                           <  1:640855/141‑1 (MQ=255)
                                                                                                               tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaaa                          >  1:336622/1‑141 (MQ=255)
                                                                                                                 tAAAGATTTAAGGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                        <  2:19548/141‑1 (MQ=14)
                                                                                                                 tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                        <  1:113418/141‑1 (MQ=255)
                                                                                                                 tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACAAGTAAAGc                        >  1:21324/1‑141 (MQ=255)
                                                                                                                   aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                      >  2:256845/1‑141 (MQ=255)
                                                                                                                   aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                      <  1:611296/141‑1 (MQ=255)
                                                                                                                   aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                      >  1:44648/1‑141 (MQ=255)
                                                                                                                   aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGGTATGCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                      >  1:478395/1‑141 (MQ=255)
                                                                                                                   aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGATATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATAAGATATACCGACGAGTAACGCTc                      >  2:191991/1‑141 (MQ=255)
                                                                                                                   aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCAAGAGGTGCCGAGTTATTCAGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTATGTATAAAATGTTATATACCGCCGAGTAAAACTc                      >  2:5020/1‑141 (MQ=17)
                                                                                                                    aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                     <  2:608419/141‑1 (MQ=255)
                                                                                                                    aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                     <  2:609759/141‑1 (MQ=255)
                                                                                                                     aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCa                    <  2:301922/141‑1 (MQ=255)
                                                                                                                        ttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAg                 <  2:702951/135‑1 (MQ=255)
                                                                                                                           atgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa              >  1:294523/4‑141 (MQ=255)
                                                                                                                           atgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAAATATGATATACCGACGAGTAAAGCTCCAAAGTTa              >  2:748000/4‑141 (MQ=255)
                                                                                                                               aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACAGACGAGTAAAGCTCCAAAGTTATTGt          >  2:284712/1‑141 (MQ=255)
                                                                                                                               aTTATAAAGACGGAGATATGTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGt          <  1:489781/141‑1 (MQ=255)
                                                                                                                                ttATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGtt         <  1:258222/141‑1 (MQ=255)
                                                                                                                                   taAAGACTGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaa      <  1:643249/141‑1 (MQ=255)
                                                                                                                                   aagaGACGGAAATAGTTCATATAATCCAGAGGTGCCGGGTTATTCGGCTAAAGATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaa      <  1:215978/138‑1 (MQ=18)
                                                                                                                                    aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaaa     >  1:388131/1‑141 (MQ=255)
                                                                                                                                       gACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGt  <  2:200481/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
AGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGT  >  USA300TCH1516_ALE/110514‑110789

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: ‑