Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 111,560 T→C 8.2% R210R (CGT→CGC)  USA300HOU_RS00520 → putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE111,5600TC 8.2% 202.6 / 4.3 75R210R (CGT→CGC) USA300HOU_RS00520putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (42/24);  new base C (4/3);  total (46/29)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATATTTACTTTAC  >  USA300TCH1516_ALE/111420‑111697
                                                                                                                                            |                                                                                                                                         
tttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACATTTACGt                                                                                                                                           >  2:11495/1‑141 (MQ=255)
 ttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAACACAATTACGTa                                                                                                                                          >  1:233021/1‑141 (MQ=255)
   aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                                                                                                        <  2:740748/141‑1 (MQ=255)
   aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaaa                                                                                                                                        >  1:95508/1‑141 (MQ=255)
   aaGTGTTTTGCGAAAAATAGCAATTTTAAAGATTTAACGAATTATAAAGAGGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                                                                                                        <  2:592748/141‑1 (MQ=14)
    aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaaaa                                                                                                                                       <  1:198212/141‑1 (MQ=255)
     gTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAg                                                                                                                                      >  2:255361/1‑141 (MQ=255)
      ttttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGa                                                                                                                                     >  1:317132/1‑141 (MQ=255)
      ttttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGa                                                                                                                                     <  1:129667/141‑1 (MQ=255)
      ttttttGCGAAATATAGTAATTTTAAAGATTTAATGAATTACAAAGTTGTTGATATATCATATAATCCAGAGGTGCCGATTTAGTCGGCTAAATATCAATTATCTAATGATGATTATAATGTAATACAATTACGCATAAGa                                                                                                                                     <  2:337879/141‑1 (MQ=17)
         tttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTACTGATGATTATAATGTAAAACAATTACGCAAAAGatat                                                                                                                                  >  1:206727/1‑141 (MQ=255)
              gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACAt                                                                                                                             >  1:73905/1‑141 (MQ=255)
               cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATa                                                                                                                            >  2:123982/1‑141 (MQ=255)
               cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATAGAGATATATGATATAATCAAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGCTTATAATGTAAAACAATTGCGCAAAAGATAGGatcta                                                                                                                            >  2:508911/1‑137 (MQ=14)
                  tataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACca                                                                                                                         <  1:438728/141‑1 (MQ=255)
                  aaatgCAATTTTAAATATTGAATGAACTATAAAGATGGAGATATATCAGATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACca                                                                                                                         <  1:198377/137‑1 (MQ=17)
                   ataGCAATTTTAAAGATTTAATGAATTTTCTATAATAAAATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATCATTATAATGTAAAACAACTACGTAAAAGACATGACATACcac                                                                                                                        <  2:731392/141‑2 (MQ=17)
                       cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATAATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACaaa                                                                                                                    >  2:272173/1‑141 (MQ=255)
                           tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAGTTAACTAATGATCATTATAATGTAAAACAATTACGTAAACGAAATGACATACCAACAAATaaa                                                                                                                >  2:525134/1‑141 (MQ=255)
                              aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATGAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcg                                                                                                             >  1:751954/1‑141 (MQ=255)
                                 gATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCg                                                                                                          >  1:247297/1‑141 (MQ=255)
                                 gATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCg                                                                                                          >  2:480718/1‑141 (MQ=255)
                                  aTTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGa                                                                                                         <  1:604907/141‑1 (MQ=255)
                                    ttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTGTTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAg                                                                                                       >  2:546142/1‑141 (MQ=255)
                                        tGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCAAAGCtgt                                                                                                   >  2:396386/1‑141 (MQ=255)
                                         gAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgtt                                                                                                  >  2:203415/1‑141 (MQ=255)
                                         gAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgtt                                                                                                  <  1:662098/141‑1 (MQ=255)
                                          aaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAACCAATTACGTAAAAGATATGAAATACCAACAAATAAAGCGCCGAAGCagttg                                                                                                 >  2:122129/1‑141 (MQ=255)
                                           aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttgt                                                                                                <  1:27577/141‑1 (MQ=255)
                                           aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttgt                                                                                                >  1:6155/1‑141 (MQ=255)
                                             tataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGACGCTGttgttg                                                                                              >  2:414828/1‑141 (MQ=255)
                                              ataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTGGTTGa                                                                                             >  1:600701/1‑141 (MQ=255)
                                                       gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAg                                                                                    >  2:516019/1‑141 (MQ=255)
                                                         gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAggg                                                                                  <  1:409492/141‑1 (MQ=255)
                                                             tatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAtt                                                                              >  1:76600/1‑141 (MQ=255)
                                                                     aaTCCAGAGGTGCCGAGTTATTCAACTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                                                                      <  2:720207/141‑1 (MQ=255)
                                                                      aTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGtt                                                                     >  1:21136/1‑141 (MQ=255)
                                                                        ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATAGCAACAAACAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                                                                   >  1:565955/1‑141 (MQ=255)
                                                                          agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGGATTTAAAGGGTtcatc                                                                 >  1:58295/1‑141 (MQ=255)
                                                                          agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACAAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTAAAGGGAATTTAAAAGGTtcatc                                                                 >  2:661092/1‑141 (MQ=255)
                                                                            agGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAg                                                               >  1:582258/1‑141 (MQ=255)
                                                                            agGTGCCGAGTTATTCAGCACAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCACCGAAGCTATTGTTGACAGGTACAGGGAATTTAAAAGGAtcatcat                                                               >  2:476703/1‑140 (MQ=14)
                                                                                gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTgg                                                           >  2:72852/1‑141 (MQ=255)
                                                                                 ccGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGa                                                          >  1:659718/1‑141 (MQ=255)
                                                                                  cGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGat                                                         >  1:24379/1‑141 (MQ=255)
                                                                                        atcaaGCTCAAAATCAACTAACTAAGGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaa                                                   <  2:515031/137‑1 (MQ=255)
                                                                                          tCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaaa                                                 >  1:274780/1‑141 (MQ=255)
                                                                                           cAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTGAAAGAGTCATCAGTTGGATATaaaaaaa                                                >  1:357798/1‑141 (MQ=255)
                                                                                           cAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGACGAAGAACAATGTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaaaa                                                >  1:594271/1‑141 (MQ=17)
                                                                                            aGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAAt                                               <  2:245834/141‑1 (MQ=255)
                                                                                            aGCTCAATAGGAATTAACTAACGATGATTATAATGTAAAACAATTACGTAAAAGGTATGACATAACAACAAATAAAGCGCCTAAGCTGTTGTTGTAAGGTACAGGGAATTTAAAGGGTTCATCAGTTGGATATAAAAAAAt                                               <  2:328394/141‑1 (MQ=17)
                                                                                                 aaTATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAAt                                          <  2:518228/141‑1 (MQ=255)
                                                                                                 aaTATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAAt                                          <  1:123982/141‑1 (MQ=255)
                                                                                                  aaatCAATTAACTAATGATGAATATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAAtt                                         <  2:703415/139‑1 (MQ=255)
                                                                                                     tCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTAc                                      <  2:32356/141‑1 (MQ=255)
                                                                                                     tCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTAc                                      >  1:227332/1‑141 (MQ=255)
                                                                                                           aaCTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACCGGGAATTGAAAAAGTTCATCAGTTGGATATAACAAAAATGAATCTACtttttt                                >  1:5272/1‑141 (MQ=17)
                                                                                                            aCTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTa                               >  1:74427/1‑141 (MQ=255)
                                                                                                              tAATGATGATTATAATGTAGAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTaga                             <  2:20646/141‑1 (MQ=255)
                                                                                                              tAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTaga                             >  2:632023/1‑141 (MQ=255)
                                                                                                                atgatgATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAgaga                           >  2:461668/1‑141 (MQ=255)
                                                                                                                  gatgatTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAAt                         <  1:661228/141‑1 (MQ=255)
                                                                                                                     gatTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATaaa                      >  1:382729/1‑141 (MQ=255)
                                                                                                                     gatTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTAAAAAGTTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATaaa                      >  1:291321/1‑141 (MQ=255)
                                                                                                                       ttATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATaaaaa                    >  2:441480/1‑141 (MQ=255)
                                                                                                                         ataATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAGGGTACAGGGAATTTAAAAGGTTCATCAATTGGATATAAAGAAATTGAATTTACTTTTTTAGGATATAAAGATg                  >  1:343121/1‑141 (MQ=14)
                                                                                                                         ataATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATg                  <  1:632023/141‑1 (MQ=255)
                                                                                                                         ataATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATg                  >  2:676810/1‑141 (MQ=255)
                                                                                                                          taATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGa                 <  2:33011/141‑1 (MQ=255)
                                                                                                                             tGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGaaaa              >  2:89272/1‑141 (MQ=255)
                                                                                                                               tAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAata            <  1:546142/141‑1 (MQ=255)
                                                                                                                                   aaaaatACGGAAAAGATATGACATACCAACAAATAAAGCGCCTAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTAGATATAAAAAAATTGAATCTACTTTTTTAGAGAATAAAAATGAAAATAttta        <  2:5272/136‑1 (MQ=18)
                                                                                                                                      aTTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAtttactt     >  1:497726/1‑141 (MQ=255)
                                                                                                                                       ttACGGAGAAGATATGACATACCAACAAAGAAAGCGCCGAAGCTGGTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAGTTTACGTTTTTAGAGAATAAAAATGAAAATAtttacttt    <  2:3767/141‑1 (MQ=18)
                                                                                                                                        tACGTAAAAGGTATGAAATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGCGAATTTAAAAGTTTAATCAGTTGCATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAtttacttta   <  2:184907/141‑1 (MQ=255)
                                                                                                                                         aCGTAAAAGATATGACATACCAACAAATAAAGCGCCAAATCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAtttactttac  <  1:480718/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                         
TTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATATTTACTTTAC  >  USA300TCH1516_ALE/111420‑111697

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