Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 396,540 G→A 25.9% A208T (GCA→ACA)  efeM → putative iron uptake system component EfeM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE396,5400GA25.9% 140.4 / 44.8 88A208T (GCA→ACA) efeMputative iron uptake system component EfeM
Reads supporting (aligned to +/- strand):  ref base G (36/28);  new base A (15/8);  total (51/37)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATACATTAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATGAGTCA  >  USA300TCH1516_ALE/396400‑396677
                                                                                                                                            |                                                                                                                                          
aTACATTAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGgcg                                                                                                                                            >  1:423695/1‑141 (MQ=255)
 tACATTAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGgcgc                                                                                                                                           >  1:603482/1‑141 (MQ=255)
 tACATTAGATATCATACAAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGgcgc                                                                                                                                           >  1:72616/1‑141 (MQ=255)
         atatCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATCCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACaaaaaa                                                                                                                                   <  1:926515/141‑1 (MQ=255)
          tatCATACCAAAAGTAATGCTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAAt                                                                                                                                  <  2:465371/141‑1 (MQ=255)
           ataaTACCAAAAATAAGGTTACAAGGTTCTGTTGACCTATTATATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAAtt                                                                                                                                 <  1:1145986/138‑1 (MQ=255)
           atCATACCAAAATTAATGTTACAAGGGTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAAtt                                                                                                                                 >  2:718263/1‑141 (MQ=255)
              aTACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTAt                                                                                                                              >  1:988921/1‑141 (MQ=255)
              aTACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTACAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTAt                                                                                                                              >  1:424007/1‑141 (MQ=255)
                aCCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGa                                                                                                                            >  1:96312/1‑141 (MQ=255)
                aCCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGa                                                                                                                            <  1:820555/141‑1 (MQ=255)
                aCCAAAATTAATGTTACAAGGTTCTGTGGACATATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAGGCAGAAATTTATTCACGTACAGATTTATATGATTTTAAAGCGAAGGTTGAAGGCGGGGAAAAAATTTATGa                                                                                                                            >  1:424531/1‑141 (MQ=255)
                  cAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACt                                                                                                                          <  1:998183/141‑1 (MQ=255)
                         aaTGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTaa                                                                                                                   >  2:1164081/1‑141 (MQ=255)
                           ttttACAAGGGTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGCTTTCATTGGATAATGAGAAGGCGCACAAAAAATTTATGACTTATTTAAAc                                                                                                                 >  2:539980/3‑141 (MQ=255)
                           tGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAAc                                                                                                                 >  1:230818/1‑141 (MQ=255)
                           tGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAAc                                                                                                                 >  1:682896/1‑141 (MQ=255)
                            gTTACGAGGTGCTGCTGACATATTAAATTAAGTTGCAACTTCTAAAATCCTAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAAcc                                                                                                                <  2:681977/141‑1 (MQ=255)
                            gTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGGGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAAcc                                                                                                                >  1:796042/1‑141 (MQ=255)
                                   ggTTCTGTTGACCTATGAACTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTa                                                                                                         <  2:977671/141‑1 (MQ=255)
                                      tCTTTTGACCTATTAACTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGACATTTATTCACATACAGATTAATATGCTTTTAAAGCGCACGTTCAAGGCACACAAAAAATTTCTGACTTATTTAAACCTATTTTagag                                                                                                      >  2:204624/1‑141 (MQ=255)
                                       cTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAgaga                                                                                                     >  1:310465/1‑141 (MQ=255)
                                          ttGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaa                                                                                                  >  2:388278/1‑141 (MQ=255)
                                          ttGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaa                                                                                                  >  1:702514/1‑141 (MQ=255)
                                              ccTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGa                                                                                             >  2:616188/1‑141 (MQ=255)
                                              ccTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGa                                                                                             <  2:72616/141‑1 (MQ=255)
                                              ccTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATGTATATGCTTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGa                                                                                             >  1:1178775/1‑141 (MQ=255)
                                                  ttAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATaaa                                                                                         >  2:1091644/1‑141 (MQ=255)
                                                      aTGAAGTAGCAACTACTAAAAAAACAAGAGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAAAATAAAAGAAAAACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAAt                                                                                     <  2:962357/141‑1 (MQ=255)
                                                        gAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTa                                                                                   >  2:765122/1‑141 (MQ=255)
                                                         aaGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTaa                                                                                  >  1:712091/1‑141 (MQ=255)
                                                          aGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAg                                                                                 >  1:454187/1‑141 (MQ=255)
                                                          aGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAg                                                                                 >  1:605414/1‑141 (MQ=255)
                                                           ggtGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGTCGGTTCAGAATTAATTTAGTACTTATTTCAGGGGTTTTTAGAGAAGAAAGTTAAAAAATTAAGt                                                                                <  2:543882/139‑1 (MQ=255)
                                                           gTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGt                                                                                <  1:1004287/141‑1 (MQ=255)
                                                              gCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgat                                                                             <  1:188791/141‑1 (MQ=255)
                                                               cAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgatg                                                                            <  1:1081188/141‑1 (MQ=255)
                                                                aaCTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGAACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgatga                                                                           <  1:333082/141‑1 (MQ=255)
                                                                 aCTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTgatgat                                                                          >  1:960774/1‑141 (MQ=255)
                                                                   ttCTAAACACACAGGAGTAGAAGAAATTTATTCACATACAGATTTATAGCCTTTGAAAGTGAACGTTGAAGGGGGCCAAAACATGTAGGACTTATTGAAGAATATTGTAGAGAAAAAAGATAAAAAATTAAGGACTGatat                                                                        >  2:314237/1‑141 (MQ=255)
                                                                   ttCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGatat                                                                        >  1:134857/1‑141 (MQ=255)
                                                                    tCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGCCTTATTTAAAACTATTTTAGAGAAAAAAGATAAAAAAATAAGTGATGatata                                                                       >  2:415265/1‑140 (MQ=255)
                                                                    tCTAAAATCACAGGCGAAGAAGAAATTTATTCACATACAGATTCATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATc                                                                       <  2:616718/141‑1 (MQ=255)
                                                                     cTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATcc                                                                      >  2:154066/1‑141 (MQ=255)
                                                                      tAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCa                                                                     <  2:230818/141‑1 (MQ=255)
                                                                       aaaaTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACACAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGGTAAAACATTAATTGATGATATCCaa                                                                    >  1:324750/1‑141 (MQ=255)
                                                                         aaTCACAGGTGAAGAAGAAATTTATACACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATg                                                                 >  1:655792/1‑141 (MQ=255)
                                                                          aTCAAAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACAAAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATg                                                                 <  2:600765/141‑1 (MQ=255)
                                                                           tCACAGGTGAAGAAGAAATTTATTCACAAACAGATTCATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGACGATATCCAAATGa                                                                >  2:589413/1‑141 (MQ=255)
                                                                             acaGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAAc                                                              >  2:440920/1‑141 (MQ=255)
                                                                             acaGGTCAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCAATAATAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAAc                                                              <  2:674285/141‑1 (MQ=255)
                                                                             acaGGCGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAAc                                                              <  1:248453/141‑1 (MQ=255)
                                                                              caGGTGAAGAAGAAATTTATTCACATACAGATTTATGTGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACt                                                             <  1:663215/141‑1 (MQ=255)
                                                                               aGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTCCATTAATCCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACtt                                                            >  1:274528/1‑141 (MQ=255)
                                                                               aGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACtt                                                            >  1:1136221/1‑141 (MQ=255)
                                                                               aGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTACTTGAGAGAAAACAGATAAAAAATTAAGTGATGATATCCAAATGAACtt                                                            >  1:865765/1‑141 (MQ=255)
                                                                               aGGTGAAGAAGAAATTTATTCACATACAGATTTATATGACTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTGAAACCTATTTTAGAGAAAAAAGATAAAAAATCAAATGATGATATCCAAATGAACtc                                                            >  2:409839/1‑140 (MQ=255)
                                                                                    aagaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATGCGACTCGTAAAAGATACAAAATTCAGTGATGATATCCAAATGAACTTCGATa                                                       >  1:1171850/1‑141 (MQ=255)
                                                                                    aagaagaaATTTATTCACACACAGATTTATATGATTTTAAAGCGAACGTTGCAGGCACACAAAAAATTTATAACTTATTCAAACCTATTTCAGAAAAAAAAGACAAAAAATTAAGTGATGATATCCAAAAGAACGTCGATa                                                       >  2:1008283/1‑141 (MQ=255)
                                                                                    aagaagaaAGTTATTCACATACAGATTTATATGATTTTAAAGCTAACGTTGAAGGCGCACAAAAAAATTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGAGATCCAAATGAACTTCGATa                                                       <  2:44659/141‑1 (MQ=255)
                                                                                     agaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATaa                                                      >  2:841407/1‑141 (MQ=255)
                                                                                      aaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCAACCAAAATCTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATaaa                                                     <  1:568446/140‑1 (MQ=255)
                                                                                         aaaaTTTATTCACATACAGATTTATATGATTTTAAAGAGAACGTTGAAGTTGCACAAAAAATTTATTAATTATTTAAACCTATTTTATAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTg                                                  <  2:486148/140‑1 (MQ=255)
                                                                                          aaaTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATATTAGAAAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGa                                                 >  2:1023319/1‑141 (MQ=255)
                                                                                             tttATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATc                                              >  1:378380/1‑141 (MQ=255)
                                                                                                  tCACATACAGATCTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAATCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAatta                                         <  2:1090177/141‑1 (MQ=255)
                                                                                                   cacaTACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAattat                                        <  1:575511/141‑1 (MQ=255)
                                                                                                   cacaTACAGATTTATATGATTTTAAAGCGAACATTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAattat                                        <  2:733452/141‑1 (MQ=255)
                                                                                                       tACAGATTTATATGATTTTAAAGCGAACGTTTAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGa                                    <  1:150267/141‑1 (MQ=255)
                                                                                                           gATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATaaa                                >  2:671059/1‑141 (MQ=255)
                                                                                                            aTTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAat                               <  2:423695/141‑1 (MQ=255)
                                                                                                               tataTGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAAtata                            >  2:369984/1‑141 (MQ=255)
                                                                                                                  atGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAg                         >  2:745092/1‑141 (MQ=255)
                                                                                                                    gATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGAt                       <  2:274528/141‑1 (MQ=255)
                                                                                                                     aTTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATCTCCAAAGGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATa                      >  1:270840/1‑141 (MQ=255)
                                                                                                                     aTTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATa                      >  2:487504/1‑141 (MQ=255)
                                                                                                                     aTTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATa                      >  1:916571/1‑141 (MQ=255)
                                                                                                                     aTTGTAAAGCGAACGTGGAAGGCGCACAAAAAATTTATGACTTTTTTAAACCTACTTTAGAGAAAAAAGATAAAAAATTAATTGATGATATCCAAATGAACTTTGATAAAGTGAATCAATTATTGGATAAATATAAAGATa                      <  2:1183968/141‑1 (MQ=255)
                                                                                                                       tttAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATaac                    <  1:371368/141‑1 (MQ=255)
                                                                                                                       tttAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGCGAATCAATTATTGGATAAATATAAAGATaac                    >  2:653846/1‑141 (MQ=255)
                                                                                                                            aGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAagg               >  2:692617/1‑141 (MQ=255)
                                                                                                                            aGCGAACGTAAAAGGAGCAAAAAAAAAACATCATTTATACAAACAAATTTCACACAAAAAAGATAAAAAATCAAGTGATGATCTCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcgg               <  1:539980/141‑1 (MQ=255)
                                                                                                                              cGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcggcg             >  2:914492/1‑141 (MQ=255)
                                                                                                                                aaagTTGAAAGCAAAAAAAAAATTTATGACCTATTTAAACCTATTTTAAAGAAAAAAGAAAAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcggcgg            <  1:1164081/138‑1 (MQ=255)
                                                                                                                                  cGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTa         <  2:157781/141‑1 (MQ=255)
                                                                                                                                  cGTTGAAGGCAAACAAAAAATTTATAACTTATTTAAACCTATTTGAGAGAAAAAAAAAAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTa         <  2:1062327/141‑1 (MQ=255)
                                                                                                                                    ttGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATACAAAATTAAGTGATGATATCCAAATGAACTTCGAAAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATg       >  2:393842/1‑141 (MQ=255)
                                                                                                                                      gAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAGATATAAAGATAACAACGGCGGttgttag     >  1:76725/1‑136 (MQ=255)
                                                                                                                                      gAAGGCGCACAAAAAATTTATGACTTATTTAAACCGATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATGAg     >  1:1123627/1‑141 (MQ=255)
                                                                                                                                         ggCGCACAAAAACTTTATGACTTATTTAACCCCATTTAAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACCTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATGAGTCa  <  1:314237/141‑1 (MQ=255)
                                                                                                                                         ggCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATGATTGGATAAATATAAAGATAACAACGGCGGTTATGAGTCa  <  2:134857/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
ATACATTAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATGAGTCA  >  USA300TCH1516_ALE/396400‑396677

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