Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 111,516 C→A 25.7% 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,5160CA25.7% 126.5 / 39.7 74Q196K (CAA→AAA) USA300HOU_RS00520putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (24/31);  new base A (11/8);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAA  >  USA300TCH1516_ALE/111380‑111649
                                                                                                                                        |                                                                                                                                     
ggAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAata                                                                                                                                   >  2:180065/1‑141 (MQ=255)
ggAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAata                                                                                                                                   >  2:313904/1‑141 (MQ=255)
      tAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGTGTGAATCGGCTAGATATCAAtt                                                                                                                             >  1:13096/1‑141 (MQ=21)
      tAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGACGAGTAATTGGGATAAATATCAAtt                                                                                                                             >  2:31809/1‑141 (MQ=18)
        aagaagAAAACATAAAAAAAGCAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGGTATTCAACACACTCTCACTTaa                                                                                                                           >  1:346839/1‑141 (MQ=21)
        aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                                                           >  2:559495/1‑141 (MQ=255)
        aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCACTTaa                                                                                                                           >  1:419299/1‑141 (MQ=255)
            agtagaCAGAAAAAAAGAAATCGAAAATTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTaa                                                                                                                       <  2:327575/136‑1 (MQ=18)
            agaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTaa                                                                                                                       <  2:569614/141‑1 (MQ=255)
                aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgat                                                                                                                   <  2:471593/141‑1 (MQ=255)
                 aCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgatg                                                                                                                  <  1:413376/141‑1 (MQ=255)
                   aTAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                                                                <  1:474936/141‑1 (MQ=255)
                    tAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGAtt                                                                                                               >  1:7140/1‑141 (MQ=255)
                     aaaaaaaGAAATCTAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATta                                                                                                              >  2:339714/1‑141 (MQ=255)
                      aaaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                                                             <  2:50110/141‑1 (MQ=255)
                         aaaGAAATCGAAAACTTTAAGTTTTTTTCTTTTTTTATCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCAAGAGGTGCCGAGTTATTCAGCTCAATATAAATTAACTATTGTTGATtttatt                                                                                                          >  2:68953/1‑139 (MQ=18)
                         aaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAAt                                                                                                          <  1:381520/141‑1 (MQ=255)
                         aaaGAAACCGAAAACTTTAAGGTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAAt                                                                                                          <  1:414212/141‑1 (MQ=255)
                           aGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAGTAATGATGATTATAATGt                                                                                                        >  2:12692/1‑141 (MQ=255)
                           aGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGt                                                                                                        <  2:263488/141‑1 (MQ=255)
                             aaaTCGAATATATTAAGTTTTTTTCGCAATATCGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCACAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                      <  1:444019/141‑1 (MQ=21)
                             aaaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                      <  1:592063/141‑1 (MQ=255)
                                 cGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                                                                  >  2:596418/1‑141 (MQ=255)
                                 cGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                                                                  >  1:381706/1‑141 (MQ=255)
                                 cGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                                                                  >  2:361153/1‑141 (MQ=255)
                                   aaaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                <  2:259776/141‑1 (MQ=255)
                                   aaaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                <  2:110156/141‑1 (MQ=255)
                                   aaaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATACCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                <  1:21495/141‑1 (MQ=255)
                                     aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATCAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTa                                                                                              >  1:354633/1‑141 (MQ=255)
                                     aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTa                                                                                              >  1:81613/1‑141 (MQ=255)
                                     aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTa                                                                                              <  1:297637/141‑1 (MQ=255)
                                           aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaaa                                                                                        <  1:14872/141‑1 (MQ=255)
                                           aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATGATAAAGATGGATATATATCATACAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAACTACGTaaa                                                                                        <  2:95781/141‑1 (MQ=255)
                                               tttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGat                                                                                    <  1:417315/141‑1 (MQ=255)
                                                   tGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGGTATTCAGCACCACAGCAAGTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGa                                                                                >  1:624404/1‑141 (MQ=21)
                                                    gcgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGAc                                                                               <  1:317732/141‑1 (MQ=255)
                                                     cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACa                                                                              >  1:208866/1‑141 (MQ=255)
                                                      gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTTATGATGATTAGAATGTAAAACAATTACGCAAAAGattttattt                                                                             >  2:539409/1‑134 (MQ=14)
                                                        aaTATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATAc                                                                           <  1:245025/141‑1 (MQ=255)
                                                           ataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAccga                                                                        <  2:110414/141‑3 (MQ=11)
                                                             aGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacg                                                                      >  2:185786/1‑140 (MQ=14)
                                                             aGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacg                                                                      >  2:653747/1‑140 (MQ=11)
                                                             aGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaaca                                                                      >  2:411212/1‑141 (MQ=255)
                                                              gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAccgaggc                                                                     >  2:326253/1‑136 (MQ=14)
                                                               taaTTGTAAAGATTTAATTAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAACGATATGACATACCAACaaa                                                                    <  2:538806/140‑1 (MQ=255)
                                                               cAATCTTAAAGATTTAATGAAATATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACaaa                                                                    <  2:232878/141‑1 (MQ=255)
                                                                 aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATa                                                                  >  2:483051/1‑141 (MQ=255)
                                                                       aaGATTTAATGAATTATAAAGATGGAAATATATCATATAATCCAGAGGTGCCGAGTAATTCAGCTCAATATCAATTAACTAATGATGATAATAATGTAAAACAATTACGTAAAAGATATGACAGACCAACAAATAAAgcgc                                                            >  2:352078/1‑141 (MQ=255)
                                                                         gATTTAAGGAATTAGAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCg                                                          <  2:184648/141‑1 (MQ=255)
                                                                          acttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGa                                                         <  1:209135/139‑1 (MQ=255)
                                                                          aTTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGa                                                         <  1:28139/141‑1 (MQ=255)
                                                                           tttAATGAATTATAAAGATGGCGCTAGTACAGATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAAAAATTACGTAAAAGATATGACATACCAAAAAATAAAacgaagaa                                                        >  2:305518/1‑136 (MQ=17)
                                                                           tttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGaa                                                        >  1:447381/1‑141 (MQ=255)
                                                                             tAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGc                                                      <  1:282181/141‑1 (MQ=255)
                                                                                          aggtGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAgg                                         <  2:615937/138‑1 (MQ=255)
                                                                                          aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAgg                                         <  2:482886/141‑1 (MQ=255)
                                                                                            aTGGAAATACATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTa                                       <  1:160889/141‑1 (MQ=255)
                                                                                               gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAg                                    >  2:262414/1‑141 (MQ=255)
                                                                                               gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCGGTTGTTGAAAGGTAAAg                                    >  2:239902/1‑141 (MQ=255)
                                                                                                   tatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGaa                                <  2:7140/141‑1 (MQ=255)
                                                                                                    atatCATATAATCCAGAGGTGCCGAGTTATTCAGATCAATATCAATTCACTATTGATGCTTCCAATGTAATCAAATTCTGCATGTGATATGACATACCAACAAAGAATTCGCCGAAGCTGTGGTTGAAAGGTACAGGGAAt                               >  2:487513/1‑141 (MQ=9)
                                                                                                       tCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAGCAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGGTGTTGAAAGGTACAGGGAATTTa                            >  1:531835/1‑141 (MQ=255)
                                                                                                        cATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTaa                           <  1:427391/141‑1 (MQ=255)
                                                                                                             aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                      <  2:310501/141‑1 (MQ=255)
                                                                                                             aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                      >  2:596411/1‑141 (MQ=255)
                                                                                                                   gagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatca                <  1:221264/141‑1 (MQ=255)
                                                                                                                      gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGtt             <  1:571796/141‑1 (MQ=255)
                                                                                                                      gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGtt             <  1:55912/141‑1 (MQ=255)
                                                                                                                         ccGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGGAGCTGTTGTTGAAAGGTACATGGAATTTAAAAGGTTCATCAGTTGGa          <  2:427202/141‑1 (MQ=255)
                                                                                                                         ccGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGa          <  1:349782/141‑1 (MQ=255)
                                                                                                                            aGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGatat       >  1:382752/1‑141 (MQ=255)
                                                                                                                             gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTCCATCAGGTGGAtata      >  2:9724/1‑141 (MQ=255)
                                                                                                                                 ttCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaa  <  1:180065/141‑1 (MQ=255)
                                                                                                                                 ttCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaa  >  1:23957/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                     
GGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAA  >  USA300TCH1516_ALE/111380‑111649

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