Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 110,660 T→A 10.5% I173I (ATT→ATA)  USA300HOU_RS00515 → putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE110,6600TA10.5% 234.2 / 13.4 84I173I (ATT→ATA) USA300HOU_RS00515putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base A (4/5);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAATTT  >  USA300TCH1516_ALE/110523‑110800
                                                                                                                                         |                                                                                                                                            
gTTCAAATGGTTGATAATCCTATTATTGCAGTAATAGAAATTAAAGATGAAAAAATAAAAAAATAAAGCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCa                                                                                                                                           <  1:399603/141‑1 (MQ=17)
        ggTCGATAATAAAATCATTCCCAGCAAATAAACTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGATCTTTGTTCAATATGGCGACTTTAAAAATTCGAAAAATTATAAAGACGCAAATATTTCATATAATcc                                                                                                                                   >  2:246579/1‑141 (MQ=17)
                  aaaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAAAGGTGCCg                                                                                                                         >  1:244570/1‑141 (MQ=255)
                    aaTCATTCCAACAAAAGAATTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAg                                                                                                                       <  1:438634/141‑1 (MQ=255)
                    aaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAg                                                                                                                       <  1:39762/141‑1 (MQ=255)
                     aTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGt                                                                                                                      >  1:471839/1‑141 (MQ=255)
                     aTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCAGAGt                                                                                                                      >  1:432499/1‑141 (MQ=255)
                     aTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCGTTGTTAAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGAAGATATTTCATATAATAAAAAAGTGCCGaca                                                                                                                      >  2:707315/1‑139 (MQ=17)
                     aTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACATTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGACATATTTCATATAAACCAAAGGTACCCAGt                                                                                                                      >  2:304257/1‑141 (MQ=17)
                       cATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTa                                                                                                                    >  2:566398/1‑141 (MQ=255)
                        aTTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTAt                                                                                                                   <  1:38200/141‑1 (MQ=255)
                         ttCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCAAAAACTTTAAGTTCTTTGTTCAATATGGCAACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGATGCCAAGTAAtt                                                                                                                  >  2:112086/1‑141 (MQ=17)
                            caacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                                                                               >  1:617078/1‑141 (MQ=255)
                            caacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                                                                               >  1:244921/1‑141 (MQ=255)
                            caacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                                                                               <  1:77139/141‑1 (MQ=255)
                               aaaaaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGTTGGCTAGTTATTCGGCTa                                                                                                            <  2:684783/140‑1 (MQ=255)
                                   aGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                                        <  1:628775/141‑1 (MQ=255)
                                   aGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTAAGTTCAATATAGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCAAGAAGTACCGAGTTATTCAGCTAAatt                                                                                                        >  2:640264/1‑140 (MQ=17)
                                   aGAAAGTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCTAGTTATTCGGCTAAata                                                                                                        <  1:779660/141‑1 (MQ=255)
                                     aaaTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATc                                                                                                      >  1:737430/1‑141 (MQ=255)
                                      aaTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGGTATTCGGCTAAATATCa                                                                                                     >  1:310560/1‑141 (MQ=255)
                                      aaTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCATTTAACAAAATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCa                                                                                                     >  1:218125/1‑141 (MQ=18)
                                       agtAAAGATGAAAAAATAAAAAAAGAAATCGAAAAGTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                                                                    <  2:178257/139‑1 (MQ=255)
                                         tAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                  >  1:318458/1‑141 (MQ=255)
                                           aaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                                <  1:56933/141‑1 (MQ=255)
                                           aaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCAGAGTTATTCGGCTAAATATCAATTaa                                                                                                >  2:53675/1‑141 (MQ=255)
                                            aGATTAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                               <  2:121106/141‑1 (MQ=255)
                                            aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                               <  2:652386/141‑1 (MQ=255)
                                            aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                               <  1:747507/141‑1 (MQ=255)
                                            aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                               <  1:747503/141‑1 (MQ=255)
                                             gATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                              <  1:676116/141‑1 (MQ=255)
                                                  aaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                                         <  2:413812/141‑1 (MQ=255)
                                                      aTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                                     >  1:133636/1‑141 (MQ=255)
                                                         aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATTTCAATTAACTAATGATGATtat                                                                                  <  1:411888/141‑1 (MQ=255)
                                                         aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                                  <  2:458742/141‑1 (MQ=255)
                                                         aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                                  >  2:740523/1‑141 (MQ=255)
                                                         aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                                  <  1:161246/141‑1 (MQ=255)
                                                         aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCAGCTAAATATCAATTAACTAATGATGATtat                                                                                  <  2:435733/141‑1 (MQ=255)
                                                             aaGAAATCTAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATg                                                                              <  2:154470/141‑1 (MQ=255)
                                                               gAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                                                                            >  1:562167/1‑141 (MQ=255)
                                                               gAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                                                                            <  1:630681/141‑1 (MQ=255)
                                                                aaaTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                           >  1:510662/1‑141 (MQ=255)
                                                                 aaTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                                                                          >  2:679055/1‑141 (MQ=255)
                                                                      aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                     <  1:5389/141‑1 (MQ=255)
                                                                      aaaaCTTTAAGTTCTTTGTTCAATATGGCCACTTTAAATATTTTAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                     <  1:25532/141‑1 (MQ=255)
                                                                      aaaaCTTTAAGGTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                     <  1:4899/141‑1 (MQ=255)
                                                                         aCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                                                  <  1:520032/141‑1 (MQ=255)
                                                                           tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                                <  1:706481/141‑1 (MQ=255)
                                                                           tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                                <  1:266949/141‑1 (MQ=255)
                                                                            ttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCa                                                               <  2:57055/141‑1 (MQ=255)
                                                                            ttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCa                                                               >  1:612513/1‑141 (MQ=255)
                                                                              aaGTTCTTTGTTGAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                             <  2:217706/141‑1 (MQ=255)
                                                                              aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCgaa                                                             >  1:662283/1‑141 (MQ=255)
                                                                               aGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaaa                                                            >  1:487178/1‑141 (MQ=255)
                                                                                       gTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgat                                                    >  1:86142/1‑141 (MQ=255)
                                                                                       gTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAAAAAATACGCAAAAgatatgat                                                    >  2:12986/1‑141 (MQ=255)
                                                                                         tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTAGTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgatat                                                  <  1:356315/141‑1 (MQ=255)
                                                                                                ggCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAACACAATTAAGCAAAAGATATGATATACcgacg                                           >  1:293091/1‑141 (MQ=255)
                                                                                                  cGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATTATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAg                                         <  1:559839/141‑1 (MQ=255)
                                                                                                     tgttAAAAATTGAATGAAGTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATAAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                                      <  2:368513/139‑1 (MQ=11)
                                                                                                      tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaaa                                     <  1:223641/141‑1 (MQ=255)
                                                                                                           aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                                <  1:170376/141‑1 (MQ=255)
                                                                                                           aaaTTCGAAAAATTATAAAGACGAAGAGATTGCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                                <  1:330195/141‑1 (MQ=255)
                                                                                                           aaaTCTCAAAAATTATAAAGAGGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                                <  1:604285/141‑1 (MQ=255)
                                                                                                            aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGGCACACAATCACGGGAAAGATATGATATACCGACGTGGAACGCTCCa                               >  1:736477/1‑141 (MQ=17)
                                                                                                             aTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaa                              <  1:349076/141‑1 (MQ=255)
                                                                                                              tttaatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaaa                             <  2:740154/134‑1 (MQ=255)
                                                                                                               ttattaaaTTATAAAAATGTAGATATATCATATAATCAAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAg                            <  1:99987/136‑1 (MQ=17)
                                                                                                                 gAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGtt                          >  1:128523/1‑141 (MQ=255)
                                                                                                                     aaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATGTGATATACCGACGAGTAAAGCTCCAAAGTTATTg                      >  1:712139/1‑141 (MQ=255)
                                                                                                                      aTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGt                     <  1:79553/141‑1 (MQ=255)
                                                                                                                       ttATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGtt                    <  1:432030/141‑1 (MQ=255)
                                                                                                                          taAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaa                 >  2:706791/1‑141 (MQ=255)
                                                                                                                           aaaGACGGAGATATTTCATATAATCCAGAGGTGCCAAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaaa                <  1:71414/141‑1 (MQ=255)
                                                                                                                            aaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAg               >  1:503554/1‑141 (MQ=255)
                                                                                                                             aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTTCCAAGTTATTGTTAAAAgg              >  2:194656/1‑141 (MQ=255)
                                                                                                                             aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTAGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTACAAAGCTATTGTTAAAAgg              >  2:748847/1‑141 (MQ=255)
                                                                                                                                 ggAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCACAGTTATTGTTACAAGCTTCa          >  2:30099/1‑141 (MQ=255)
                                                                                                                                  gagaTATTTCATATAATACAGAGGTGACGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAg         <  2:571249/141‑1 (MQ=255)
                                                                                                                                   agaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAgg        >  2:443669/1‑141 (MQ=255)
                                                                                                                                   agaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAgg        >  2:262874/1‑141 (MQ=255)
                                                                                                                                    gaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAggg       <  1:192973/141‑1 (MQ=255)
                                                                                                                                     atatTTCATATAATCCAGAGGTGCCGAGATATTCGTCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACAGACGAGTAAAGCTCCAAAGTAATTGTTAAAAAGTTCAGGGa      >  2:407136/1‑141 (MQ=18)
                                                                                                                                      tataTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGaa     <  2:583302/141‑1 (MQ=255)
                                                                                                                                       atTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAAt    >  1:598046/1‑141 (MQ=255)
                                                                                                                                         ttCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAAttt  <  1:441788/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                            
GTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAATTT  >  USA300TCH1516_ALE/110523‑110800

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