Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,372,559 C→T 5.5% C406Y (TGT→TAT)  pbuG ← Guanine/hypoxanthine permease PbuG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,372,5590CT 5.5% 210.8 / 2.7 73C406Y (TGT→TAT) pbuGGuanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base T (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

AGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGG  >  USA300TCH1516_ALE/2372420‑2372698
                                                                                                                                           |                                                                                                                                           
aGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAGTATGAAGTACAGAAGGAATAGTACCATTAATACATACATTATCGGTTGTACTGCTTTATGTTGTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCcaca                                                                                                                                            >  2:574375/1‑141 (MQ=255)
 gTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAACCATTGTAATTGGATAGAAAATAAATCCccac                                                                                                                                           >  2:42178/1‑141 (MQ=255)
     ttttttGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCaa                                                                                                                                       >  2:398149/1‑141 (MQ=255)
     ttttttGTGTCTGGATTTAACCAAGAACAAATATGAAGTACAAAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCaa                                                                                                                                       <  1:738116/141‑1 (MQ=255)
      tttttGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAAt                                                                                                                                      >  1:499640/1‑141 (MQ=255)
      tttttGTGTCTGGATTTAACCATGAACAAATACGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTGTGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAAt                                                                                                                                      <  1:527089/141‑1 (MQ=255)
             gtCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGtt                                                                                                                               <  2:676857/141‑1 (MQ=255)
             gtCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGGACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGtt                                                                                                                               <  1:2293/141‑1 (MQ=255)
             gtCTGGATTTAACCACGAACAATTATTAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATCTTTTTCAGAGATACGCATTGTAATTGGATAGAAAATAGATCCACACGCAATACCAGtt                                                                                                                               <  1:479138/141‑1 (MQ=255)
              tCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTg                                                                                                                              <  1:175068/141‑1 (MQ=255)
              tCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTg                                                                                                                              >  1:159507/1‑141 (MQ=255)
               cTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGc                                                                                                                             <  1:268712/141‑1 (MQ=255)
                  caatcaaCCATGAACCACTATGACGGAAACAATGAGAATCACCATGAACACACACAGTATCGGATGTACTTTTTTATGTTTTTGAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAAt                                                                                                                          <  2:743661/136‑1 (MQ=255)
                   aTTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAAtt                                                                                                                         <  2:564663/141‑1 (MQ=255)
                           aTGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGa                                                                                                                 <  2:773212/141‑1 (MQ=255)
                             gAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATa                                                                                                               >  1:261282/1‑141 (MQ=255)
                              aaCAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATaa                                                                                                              <  2:203159/141‑1 (MQ=255)
                                cAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAg                                                                                                            <  2:108135/141‑1 (MQ=255)
                                     atGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGcatc                                                                                                       <  1:464629/141‑1 (MQ=255)
                                     atGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGcatc                                                                                                       <  1:448463/141‑1 (MQ=255)
                                      tGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCATACGCAATACCAGTTGCAATTGAATATGATAAAGGcatca                                                                                                      <  2:462800/141‑1 (MQ=255)
                                       gAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGcatcat                                                                                                     <  2:105817/141‑1 (MQ=255)
                                       gAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGAAATACCAGTTGCAATTGAATATGATAAATGcagcat                                                                                                     >  2:662945/1‑141 (MQ=255)
                                         aGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATaa                                                                                                   <  1:706090/141‑1 (MQ=255)
                                               gaatgaatAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCATACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTg                                                                                             <  1:281411/141‑1 (MQ=255)
                                                        gTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAgg                                                                                    <  2:680967/141‑1 (MQ=255)
                                                         tACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGc                                                                                   >  1:675828/1‑141 (MQ=255)
                                                           ccATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTAGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTg                                                                                 >  1:698186/1‑141 (MQ=255)
                                                                aaTACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGAAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGTTAAAGGCATCATAATGATTGTAATGAAAGCTGGCACt                                                                            >  1:700066/1‑141 (MQ=255)
                                                                    catacatTATCGGATGTACTTCTTTATGTTTTTTAGAGCTTATCATGGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCaa                                                                        >  2:720095/1‑141 (MQ=255)
                                                                          ttATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCaa                                                                  >  1:72041/1‑141 (MQ=255)
                                                                             tCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCAACACGCAATACCAGTTGCAATTGAATATGATATTGGCGTCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAAc                                                               >  1:298967/1‑141 (MQ=255)
                                                                               ggATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGTCACTGCAACTTCACAACtt                                                             >  1:17763/1‑141 (MQ=255)
                                                                               ggATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATGGTAATGAAGGCTGGCACTGCAACTTCAAAACtt                                                             >  1:17577/1‑141 (MQ=255)
                                                                                   gTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTcc                                                         >  2:69966/1‑141 (MQ=255)
                                                                                   gTACTACTTAATGTTTTTTTGATATAGGCATTGCAAGAGGGAAGAGACTAAATTCACAGGCAATACCAGTTGCAATTGAATATGATCAAGGCATAATAATAATTGTAAAGAAGGCTAGAATAGCAACTTCAAAACTCTTcc                                                         >  2:385952/1‑141 (MQ=255)
                                                                                    tACTTCTTTATGTTTTTTAGATATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCTGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGAAACTTCAAAACTTTTCCa                                                        <  1:315672/141‑1 (MQ=255)
                                                                                    tACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGAAACTTCAAAACTTTTCCa                                                        <  1:734356/141‑1 (MQ=255)
                                                                                    tACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAAAAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGCAGGCTGGCACTACGACTTCCAAAATTTTCCa                                                        >  2:512086/1‑141 (MQ=255)
                                                                                       ttcttTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGtt                                                     <  1:398149/141‑1 (MQ=255)
                                                                                                ttttttAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGc                                            >  2:335533/1‑141 (MQ=255)
                                                                                                    ttAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAg                                        >  2:30560/1‑141 (MQ=255)
                                                                                                    ttAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAg                                        <  2:216516/141‑1 (MQ=255)
                                                                                                      agagATAAGCATTGTAATTGGATAGAAAATAAATCCATACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGtt                                      >  1:444872/1‑141 (MQ=255)
                                                                                                        agaTAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTg                                    >  1:430764/1‑141 (MQ=255)
                                                                                                         gaTAAGCATTGTAATTTGATAGAAAATAAATCCACACGCAATACCAGCTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAGGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGc                                   <  2:325636/141‑1 (MQ=255)
                                                                                                              gCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCa                              <  2:157109/141‑1 (MQ=255)
                                                                                                               cATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAACATTTTCCGGTTTATTCCAGCACAGTTTGCAGCAAt                             >  1:322367/1‑141 (MQ=255)
                                                                                                               cATTATAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCAt                             >  2:546401/1‑141 (MQ=255)
                                                                                                                aTTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCAtt                            >  2:748342/1‑141 (MQ=255)
                                                                                                                 ttGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTa                           >  1:321212/1‑141 (MQ=255)
                                                                                                                 ttGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTa                           >  1:719782/1‑141 (MQ=255)
                                                                                                                 ttGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTAGCTCAGCAAAGTTTGCAGCCATTa                           >  2:321686/1‑141 (MQ=255)
                                                                                                                   gTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCGGCCATTAAt                         >  1:313136/1‑141 (MQ=255)
                                                                                                                   gTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAAt                         <  2:575583/141‑1 (MQ=255)
                                                                                                                     aaTTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATac                       <  2:646947/141‑1 (MQ=255)
                                                                                                                      aTTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATaca                      <  1:651278/141‑1 (MQ=255)
                                                                                                                       tttgATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATacac                     <  2:232650/138‑1 (MQ=255)
                                                                                                                       ttGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATacac                     >  1:707379/1‑141 (MQ=255)
                                                                                                                       ttGGATAGAAAATAAAACCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATacac                     <  2:643009/141‑1 (MQ=255)
                                                                                                                       ttGGATACAAAATAAATCCACACGCAATACCAGTTGTAATTCAATATGCTAAAGGCATCATAATGATTGTAACGAAGGCTGGCCCTGCAACTTCAAAACTTCTCCAATTTATTTCTGCAAAGTTTGCAGCCATTAATacaa                     <  2:5563/141‑2 (MQ=255)
                                                                                                                       ttGGAGAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAATGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATacac                     <  2:545710/141‑1 (MQ=255)
                                                                                                                        tGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACAcc                    >  1:560927/1‑141 (MQ=255)
                                                                                                                        tGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAAAACAcc                    >  1:245730/1‑141 (MQ=255)
                                                                                                                        tGGATAGAAAACAAATCCACACGGAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACAcc                    <  1:230256/141‑1 (MQ=255)
                                                                                                                         ggATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAGCGGCATCAGAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACGTTTCCATTTTGTCTCAGCAAAGTTTGCAGCCATTCATACACCt                   >  1:263931/1‑141 (MQ=255)
                                                                                                                         ggATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCt                   <  1:693566/141‑1 (MQ=255)
                                                                                                                             gtaaaaTAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTacca               <  2:740191/139‑1 (MQ=255)
                                                                                                                             aGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTacca               <  1:42178/141‑1 (MQ=255)
                                                                                                                               aaaaTAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTaccacc             <  2:624652/141‑1 (MQ=255)
                                                                                                                                aaataaatCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTAccacca            >  1:201193/1‑141 (MQ=255)
                                                                                                                                      atCCATACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTCTTCCAGTTTATCTCAGCAACGTTTGCAGCCATTAATACAACTACCACCACTAAgg      >  2:713386/1‑141 (MQ=255)
                                                                                                                                      atCCACACGCAAGTCCGGTTGCAATCGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAgg      <  2:196523/141‑1 (MQ=255)
                                                                                                                                       tCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGc     <  2:735243/141‑1 (MQ=255)
                                                                                                                                        ccACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCt    >  2:767657/1‑141 (MQ=255)
                                                                                                                                          acacGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGAATGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTGCCACCACTACGGCTgg  >  2:510110/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                           
AGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGG  >  USA300TCH1516_ALE/2372420‑2372698

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: