Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,141,235 C→A 7.4% S285* (TCA→TAA)  USA300HOU_RS05655 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,141,2350CA 7.4% 620.0 / 8.5 247S285* (TCA→TAA) USA300HOU_RS05655hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (110/107);  new base A (7/19);  total (119/128)
Fisher's exact test for biased strand distribution p-value = 2.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

CAAGACAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATAAAAAGCA  >  USA300TCH1516_ALE/1141102‑1141367
                                                                                                                                     |                                                                                                                                    
caaGACAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACAGAAACATAAAGATGAAAAATc                                                                                                                                      >  1:2871720/1‑134 (MQ=255)
  aGACAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCaa                                                                                                                                    >  1:2335002/1‑134 (MQ=255)
   gACAGTAATAAAAAGATTGAAAATAAGCTACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCAAACAATAAAAAAGGTTAAAAATAATCACAAAAAACAACCGAAACATAAAGATGAAAAATCaaa                                                                                                                                   >  2:1760865/1‑134 (MQ=255)
       gTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAg                                                                                                                               >  1:2260684/1‑134 (MQ=255)
        taagaaAAAGATTGAAAATGAGCGACCTAAAGCATAAGGAACATTAAAAGTTGAAAATAGACCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAgg                                                                                                                              <  2:2841290/134‑1 (MQ=255)
         aataaAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGga                                                                                                                             >  2:114156/1‑134 (MQ=255)
          ataaAGAGATGGAAAATGAGCGACCTAAAGCATCAAGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaa                                                                                                                            <  2:573334/134‑1 (MQ=255)
          ataaAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaa                                                                                                                            <  1:1433403/134‑1 (MQ=255)
          ataaAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCATCCAACAAGAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaa                                                                                                                            <  1:1001842/134‑1 (MQ=255)
            aaaaaGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaaaa                                                                                                                          >  2:2635734/1‑134 (MQ=255)
             aaaaGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACAGAAACATAAAGATGAAAAATCAAAAAAGgaaaag                                                                                                                         >  1:1489586/1‑134 (MQ=255)
               aaGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaaaagaa                                                                                                                       >  2:939737/1‑134 (MQ=255)
                 gATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaaaagaaaa                                                                                                                     >  2:1650533/1‑134 (MQ=255)
                 gATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaaaagaaaa                                                                                                                     >  1:2537459/1‑134 (MQ=255)
                    tGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGTCCTGCAACAAAAAAAAAAGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGAAAAAGAAAAAag                                                                                                                  >  1:2571163/1‑134 (MQ=255)
                    tGAAAATAAGCAACTTATAGTATTATTTATATAAAAAGTTGAAAATAGCCCTCGAAAAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGAAAAAAAAAAAAAAAAGGAAAAGAAAAAag                                                                                                                  <  1:1017459/134‑1 (MQ=255)
                      aaaaTGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAACAAAAAagta                                                                                                                >  2:2741780/1‑134 (MQ=255)
                      aaaaAGAGCGAACTAAAGGATCAAGTACATTAAAAGTTGAAAATAAAAATACAAAAAAAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGaaaaaaaaa                                                                                                                <  2:762950/134‑4 (MQ=255)
                       aaaTGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAagtag                                                                                                               >  2:969934/1‑134 (MQ=255)
                       aaaTGAGCGAACTAAAGCATCAGGTAAATTAAAAGTTGAAAATAGCCCTCCAATAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAagtag                                                                                                               <  2:2208297/134‑1 (MQ=255)
                        aaTGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAagtagt                                                                                                              >  1:2807116/1‑134 (MQ=255)
                          tGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTg                                                                                                            >  1:1459088/1‑134 (MQ=255)
                          tGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATAACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAAGAAAAGAAAAAAGTAGTTg                                                                                                            >  2:1230574/1‑134 (MQ=255)
                           gAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTga                                                                                                           >  1:897537/1‑134 (MQ=255)
                            aGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaa                                                                                                          >  1:165734/1‑134 (MQ=255)
                              cGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACATTGAGCATTCTAAACTTAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaaaa                                                                                                        <  1:2650283/134‑1 (MQ=255)
                                aCCAAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaaaaag                                                                                                      <  1:1514264/134‑1 (MQ=255)
                                    aaaGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaaaaagaaaa                                                                                                  >  1:1150954/1‑134 (MQ=255)
                                     aaGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAATAATAATATAGTTTGAAAATAATCACAAAGAGCAAAAGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaaaaagaaaaa                                                                                                 <  1:192617/134‑1 (MQ=255)
                                        tattaGGTACATTAAAAGTTGAAAATAGCCCTCCAAAAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCg                                                                                              <  1:2062610/130‑1 (MQ=255)
                                           cAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAATATCTACTGAAACATTAAGATGAAAAATAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTa                                                                                           <  1:53333/134‑1 (MQ=255)
                                           aaGGAACATTAAAAATTAAAAAAAGCCCACCAACAAAAAAAAAAGTTAAAAATAATAAAAAAGAGCAACCGAAACATAAAAAAAAAAAAAAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTa                                                                                           <  2:230812/133‑1 (MQ=255)
                                            aGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTAc                                                                                          >  1:2642655/1‑134 (MQ=255)
                                             ggTACATTAAAAATTGAAAATAGCCGTCCAAAAAAAAAAAAGGTTAAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAATGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTAcc                                                                                         <  2:661631/134‑1 (MQ=255)
                                             agTACAATAAAAGTTAAAAATAGCCCTCCAACAATAAAAAAGGTTGAAATTAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTAcc                                                                                         <  1:2307238/133‑1 (MQ=255)
                                              gTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCa                                                                                        >  1:213885/1‑134 (MQ=255)
                                                 cATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCt                                                                                     <  1:1554700/134‑1 (MQ=255)
                                                   ttAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCttt                                                                                   <  1:1864890/134‑1 (MQ=255)
                                                      aaaGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTaa                                                                                >  1:1330412/1‑134 (MQ=255)
                                                       aaGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATGACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTaaa                                                                               <  2:2675640/134‑2 (MQ=255)
                                                       aaGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAATTGAAAAAGAAAAAGCGTTACCAGCTTTTAAt                                                                               >  2:550076/1‑134 (MQ=255)
                                                       aaGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATAAAAAAAGGAAAAGAAAAGAGAAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAAt                                                                               >  1:1637605/1‑134 (MQ=255)
                                                        aGTTGAAAATAGCCCTCCAACAACAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATa                                                                              >  2:234696/1‑134 (MQ=255)
                                                        aGTAGAAAATAGCCATACAAAAATAAAAAAAGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATa                                                                              <  2:76917/134‑1 (MQ=255)
                                                          ttGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAACTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATaga                                                                            <  1:1742876/134‑1 (MQ=255)
                                                            gAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAgaga                                                                          >  1:899805/1‑134 (MQ=255)
                                                            gAAAATAGCACTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAATAGATAGAGTATAAAAAGCGTTACCAGCTTTTAATAgaga                                                                          >  2:597778/1‑134 (MQ=255)
                                                              aaaTAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAAATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAgatg                                                                        >  1:223713/1‑134 (MQ=255)
                                                               aaTAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAATAGTAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAGAGCGTTACCAGCTTTTAATAGAgatga                                                                       <  1:401407/134‑1 (MQ=255)
                                                                aTAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAgaagat                                                                      >  2:2477438/1‑134 (MQ=255)
                                                                  aGCCCTCCAATAATAAAAAAGGTTGTTTTTAATATAAAAAAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAgatgatga                                                                    <  1:1221752/134‑1 (MQ=255)
                                                                  aGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAgatgatga                                                                    <  1:2171112/134‑1 (MQ=255)
                                                                   gcacTCCAAAAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAgatgatgat                                                                   <  2:324199/131‑1 (MQ=255)
                                                                     catacaacaaTAAAAAAGGTTGAAAATAATAACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAg                                                                 <  2:2537459/130‑1 (MQ=255)
                                                                       tCCAAAAATAAAAAAGGTTAAAAATAATCACAAAGAGCAATCGAAACATAAAGATTAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCa                                                               <  1:996733/134‑1 (MQ=255)
                                                                        ccAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGATGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCaa                                                              >  1:2820440/1‑134 (MQ=255)
                                                                            caaTAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAACCGGAAAAGACAAAAGCAGCAGAAAAAGAAAAAGCCTTACCAGCTTTTAATAGAGATGAGGATAGCAAGAAt                                                          >  2:1177182/1‑134 (MQ=255)
                                                                            caaTAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAAt                                                          <  1:2351247/134‑1 (MQ=255)
                                                                             aagaaaaaaGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAAta                                                         <  2:1733107/131‑1 (MQ=255)
                                                                               tAAAAAGGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCACAAAAGGCAAAGAAAAAAATAGTTGAAAAAGAAAAAGCGATAAAAGCTTATAATAGAGATGATGATAGCAAGAAtagt                                                       >  1:1022553/1‑134 (MQ=255)
                                                                                aaaaaaGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAAtagta                                                      >  2:1977467/1‑134 (MQ=255)
                                                                                 aaaaaaGTATAAAAGAATCACAAAAAGCAACCGAAACATAAAGATGAAAAATAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAAtagtag                                                     <  1:477118/134‑1 (MQ=255)
                                                                                 aaaaaaGATGAAAATAATCAAAAAGAGAAACAGAAAAATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAAtagtag                                                     <  1:345593/134‑1 (MQ=255)
                                                                                   aaaGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTc                                                   >  1:1573426/1‑134 (MQ=255)
                                                                                     aGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCaa                                                 <  1:753930/134‑1 (MQ=255)
                                                                                     aGGTTAAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCaa                                                 <  1:1895971/134‑1 (MQ=255)
                                                                                      ggTTGAAAATAATCACAAAGAGCAACCGAAACACAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAACAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATACTAGACAAt                                                >  2:536612/1‑134 (MQ=255)
                                                                                        ttGAAAATAATCACAAAAAGACACCAAAACATAAAGATTAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTa                                              <  2:1499594/134‑1 (MQ=255)
                                                                                          gAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTCGAAAAAGAAAAAGCGTTACAAGCTTTTAATAGAAATAATGATAACAAGAATAATAGTCAATTATc                                            >  1:2051472/1‑134 (MQ=255)
                                                                                           aaaaTAATCACAAAAAGCAACCGAAACATAAAAATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCt                                           <  1:1094186/134‑1 (MQ=255)
                                                                                           aaaaTAATAAGAAAGAGAGAACGAAATATAAAGATGGAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAGAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCt                                           <  2:1331422/134‑1 (MQ=255)
                                                                                             aataatCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAg                                         <  2:2494972/134‑1 (MQ=255)
                                                                                                  taacaAAGAGTAACCGAAACATAAAGATAAAAAATAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGata                                    <  1:596989/132‑1 (MQ=255)
                                                                                                    acaAAGTTCAAGCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTAACCAGCAATTAATAGAGATGATGATAGCAAGAATAGTAGTAAATTATCTAGTAatata                                  >  2:502083/1‑133 (MQ=255)
                                                                                                     tcaaGACCAACACAATCATTAAGAGAAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTa                                 <  1:2082557/132‑1 (MQ=255)
                                                                                                     gaaaGAGCAAACGAAACATAAAGATGAAAAATAAAAAAAGGAAAAGAAAAAAGTAGTTTAAAAAGAAAAATCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTa                                 <  1:488814/133‑1 (MQ=255)
                                                                                                      aaaGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAATCGTTATCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTaa                                <  1:18036/134‑1 (MQ=255)
                                                                                                        agagAAACCGAAACATAAAGAAGAAAAAAAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAg                              <  2:225827/134‑1 (MQ=255)
                                                                                                        aaagtaaCAGAAACATAAAGATGAAAAATAAAAAAAGGAAAAGAAAAAAGTAAAAAAAAAAAAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAg                              <  2:2097556/129‑1 (MQ=255)
                                                                                                         gagCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAAAGTTACCAGCTATAAATAGAGATGAAGATAGCGAGAATAGTAGTCAATTATCTAGTGATATTAAAGa                             >  1:96065/1‑134 (MQ=255)
                                                                                                          agCACCCGAAACATAAAGCTGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGaa                            <  1:520044/134‑1 (MQ=255)
                                                                                                          agCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGAGTTCCGAAAGTTTAAAAAAAATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGaa                            >  1:667045/1‑134 (MQ=255)
                                                                                                            aaaCCGAAAAATAAAGAAGAAAAAAAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACt                          <  2:1839163/133‑1 (MQ=255)
                                                                                                               ccGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTtga                       >  1:1132838/1‑134 (MQ=255)
                                                                                                                cGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTtgat                      >  2:1478318/1‑134 (MQ=255)
                                                                                                                 gAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTtgatg                     <  2:2162890/134‑1 (MQ=255)
                                                                                                                 aaaaGATAAAAAAAAAAAAAAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTtgatg                     <  2:1373794/133‑1 (MQ=255)
                                                                                                                  aaaCATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTtgataa                    >  2:823683/1‑132 (MQ=255)
                                                                                                                  aaaCATAAAGATGAAAAATAAAATATGTTAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTtgatga                    <  2:1395977/134‑1 (MQ=255)
                                                                                                                        aaaaATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCaaa              <  2:474029/134‑1 (MQ=255)
                                                                                                                        aaaGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGCTGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCaaa              <  1:1177182/134‑1 (MQ=255)
                                                                                                                          aGATGAAAAATCAAAAAATTAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATc            <  2:86182/134‑1 (MQ=255)
                                                                                                                           gtggAAAAAGCAAAAAAGTAAAAGAAGAAAGGAATTGAAAAAGAAAAAGTGTTACCAGCTTTTAATAGAGATGATGAGAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCa           <  2:279165/131‑1 (MQ=255)
                                                                                                                             tGAAAAATCAAAAAAAAAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATa         <  2:1689543/134‑1 (MQ=255)
                                                                                                                             tGAAAAATAAAAAAAAGAAAAAAAAAAAGTAGTTGAAAAAGAAAAAGCGTAAGCTGCGTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATa         <  2:1021163/134‑1 (MQ=255)
                                                                                                                               aaaaaTCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATaaa       >  2:1094700/1‑134 (MQ=255)
                                                                                                                               aaaaaTCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATaaa       >  1:1685950/1‑134 (MQ=255)
                                                                                                                                aaaaTCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATaaaa      >  1:2421975/1‑134 (MQ=255)
                                                                                                                                 aaaTCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATaaaaa     >  1:2747785/1‑134 (MQ=255)
                                                                                                                                    tCAAAAAAGGAAAAGAAACAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGCGATATTAAACAACTTGAGGAACAAAATCATAAAAAGcc  >  2:1917370/1‑133 (MQ=255)
                                                                                                                                     |                                                                                                                                    
CAAGACAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATAAAAAGCA  >  USA300TCH1516_ALE/1141102‑1141367

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