Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,372,586 G→T 19.5% S397* (TCA→TAA)  pbuG ← Guanine/hypoxanthine permease PbuG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,372,5860GT19.5% 147.1 / 19.4 69S397* (TCA→TAA) pbuGGuanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base T (6/8);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACA  >  USA300TCH1516_ALE/2372449‑2372723
                                                                                                                                         |                                                                                                                                         
gAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGTGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGAt                                                                                                                                         >  2:213518/1‑140 (MQ=255)
gAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATTCCATTTGCCTTTTAATATgtt                                                                                                                                         >  2:378623/1‑138 (MQ=255)
     aaTATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAgg                                                                                                                                    <  1:610607/140‑1 (MQ=255)
     aaTATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAgg                                                                                                                                    <  1:419730/140‑1 (MQ=255)
     aaTATGAAGGACAGAATGAATAGTACCATTAATAAATATATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAgg                                                                                                                                    <  1:295184/140‑1 (MQ=255)
           aaGTACATAATTACTAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGcatcat                                                                                                                              <  2:678877/140‑1 (MQ=255)
           aaGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGCATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGcatcat                                                                                                                              >  2:207588/1‑140 (MQ=255)
              tACAGAATGAATAGTACCATTAATACATTCATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATGTGATAAAGGCATCATAAt                                                                                                                           <  1:311833/140‑1 (MQ=255)
               aCAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGCCATCAAAATg                                                                                                                          >  2:683087/1‑140 (MQ=255)
                cAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGa                                                                                                                         <  2:450442/140‑1 (MQ=255)
                  gaatgaatAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGAtt                                                                                                                       <  2:27799/140‑1 (MQ=255)
                  gaatgaatAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGAtt                                                                                                                       <  2:366969/140‑1 (MQ=255)
                   gatgaatAGTAACATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATTATAAAGGCATCATAATGATTg                                                                                                                      <  2:660368/139‑1 (MQ=255)
                    atgaatAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGt                                                                                                                     >  1:692643/1‑140 (MQ=255)
                          aGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGAAACATAATGATTGTACTgca                                                                                                               >  2:419397/1‑138 (MQ=255)
                           gTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAg                                                                                                              >  1:5789/1‑140 (MQ=255)
                           caaCCACTAATACATACATTATCAGATGTACTTCTTTATGATTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATTATAAAGGCATCATAATGATTGTAATGAAg                                                                                                              <  1:236531/138‑1 (MQ=255)
                                 ttAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGc                                                                                                        >  1:500427/1‑140 (MQ=255)
                                     tacatacaTTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATTATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTg                                                                                                    >  2:208119/1‑140 (MQ=255)
                                          ccatTATCGGATGTACGTCATTATGTTATTTAGAGAAAAGCATTGTAATTGGCTAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATTAAGGCTGGCACTGCAACt                                                                                               <  2:461090/139‑1 (MQ=255)
                                          acatTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTTTATATAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATTAAGGCTGGCACTGCAACt                                                                                               <  2:497462/140‑1 (MQ=255)
                                          acatTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATTAAGGCTGGCACTGCAACt                                                                                               <  2:483743/140‑1 (MQ=255)
                                           catTATCGGATGTGCTTTTTTATGTTGTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACtt                                                                                              <  2:544495/140‑1 (MQ=255)
                                             ttATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCa                                                                                            >  2:625195/1‑140 (MQ=255)
                                                 cGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAAc                                                                                        <  2:687517/140‑1 (MQ=255)
                                                  ggATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACt                                                                                       >  1:657861/1‑140 (MQ=255)
                                                     tGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACtttt                                                                                    <  1:176440/140‑1 (MQ=255)
                                                        aCTTATTTATGTTTTTTATAGATTTGTATTGTAATTGGATATAAATTATTTCCACACGCAATACCATTTGCAATTGAATATTATCAAGGCCTCATAATGATTGTAATCAAGGCTGGCACTGCAACTTCAAACATTTTcta                                                                                 >  2:255716/1‑138 (MQ=255)
                                                         cttcttTATGTTTTTTAGAGATAAGTATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAg                                                                                <  1:146519/140‑1 (MQ=255)
                                                         cttcttTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAg                                                                                <  1:374263/140‑1 (MQ=255)
                                                         cttcttGAGGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAg                                                                                <  2:558911/140‑1 (MQ=255)
                                                            cttTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATTATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGttt                                                                             >  1:537622/1‑140 (MQ=255)
                                                            cttGATGTTTTTTAGAGAGAAGCATTGTAATTGGATAAAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTATTGTAGGCTGGCACTGCAACTTCAAAACTTTTCCAGttt                                                                             <  2:126694/140‑1 (MQ=255)
                                                                   ttttttAGAGATAAGCATTGTAATTGGATTGAAAATAAATCTACCCGCGATACCATTTGCAATTGCATATGATAAAGACATCATAACGGCTGTCTTGAGGGCAAAAGCTTCAAATTCTAAACTTTGCCAGTTTAactccc                                                                      >  2:646328/1‑138 (MQ=255)
                                                                      tttCGACATAAGCGATGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATTATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCaa                                                                   <  2:415993/140‑1 (MQ=255)
                                                                      tttAGATATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACAAGTTGCAATTGAATATTATAAAGGCATCATAATGATTGTAATTAAGGCTGGCACTGCAACTTAAAAAAATTTCCAGTTTAAATAAGCaa                                                                   >  2:40271/1‑140 (MQ=255)
                                                                        tAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAg                                                                 >  1:310898/1‑140 (MQ=255)
                                                                                aGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATTATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGc                                                         <  1:400391/140‑1 (MQ=255)
                                                                                 gCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGcc                                                        <  1:482232/140‑1 (MQ=255)
                                                                                  cATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCCCTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAACGTTCGCAGCCa                                                       >  2:214358/1‑140 (MQ=255)
                                                                                   aTTGTAATTGGATACAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATTATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCAt                                                      <  2:695126/140‑1 (MQ=255)
                                                                                            ggATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACAcc                                             <  2:328918/140‑1 (MQ=255)
                                                                                             gATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATTATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCt                                            >  2:89273/1‑140 (MQ=255)
                                                                                              aTAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATTATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACAACTa                                           >  2:345069/1‑140 (MQ=255)
                                                                                               tAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTac                                          >  1:235153/1‑140 (MQ=255)
                                                                                                   aaataaatCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTaccacc                                      >  1:177157/1‑140 (MQ=255)
                                                                                                    aataaacCCACAAGCAATACCAGTTGCAATTGAATATTATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTAccacca                                     <  1:47623/140‑1 (MQ=255)
                                                                                                       aaatCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTACTACACCTACCCCCACTa                                  >  2:658294/1‑140 (MQ=255)
                                                                                                       aaatCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGAAGCCATTAATACACCTACCACCACTa                                  >  2:16723/1‑140 (MQ=255)
                                                                                                         atCCACACGCAATACAAGTTGCAATTGAATATGAGAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAg                                >  1:485276/1‑140 (MQ=255)
                                                                                                           ccACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTCATACACCTACCACCACTAAGGc                              >  1:89371/1‑140 (MQ=255)
                                                                                                            cacaCGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTATCCAGGTTATCTCAGCAAAGTTTGCAGCATCGGATACACCTACCACCTCTATGGCt                             >  1:292592/1‑140 (MQ=255)
                                                                                                             acacGCAATACCAGTTGCAATTGAATATGAGAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAAGTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTCGCAGCCATTAATACACCTACCACCACTAAGGCTg                            <  1:467653/140‑1 (MQ=255)
                                                                                                                cGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAGAGTTTGCAGCCATTAATCAACCTACCACCCCTAAGGctgctg                         >  2:6047/1‑140 (MQ=255)
                                                                                                                cGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGgtg                         >  1:654488/1‑140 (MQ=255)
                                                                                                                cGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGgtg                         >  1:179270/1‑140 (MQ=255)
                                                                                                                       gcAGTTGCAATTGAATATGATAAAGGCATCAGAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAAc                  <  2:542318/139‑1 (MQ=255)
                                                                                                                       ccggtcgCAATTGAATATGATAAAGGCACCATGATGATTGTAATGAAGGCTGGCACTGCAACCACAAAACTTTTCCAGTGTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAAc                  <  1:442723/134‑1 (MQ=255)
                                                                                                                        cAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTaagt                 >  2:438700/1‑138 (MQ=255)
                                                                                                                         aGTTGCAATTGAATATCATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAATGCTGGTGTCGTAACTg                <  1:270312/140‑1 (MQ=255)
                                                                                                                          gTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTATCCAGTTTATCTCAGCAAAGTTTACAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGc               >  2:567573/1‑140 (MQ=255)
                                                                                                                            tGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTGTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCAc             <  2:44457/140‑1 (MQ=255)
                                                                                                                            tGCAATTCGACATGAGAAAGGCATCATAATGATTGTAATGAAGGCTGGCACCGCAACTTCAAAACTTTTCCAGTCTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCAc             <  1:416547/140‑1 (MQ=255)
                                                                                                                             gCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACt            >  1:163690/1‑140 (MQ=255)
                                                                                                                                tatGAATATTATATACGCATCATAATGATTATAATGAAGGCTGTCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCATCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGt         <  2:664429/138‑1 (MQ=255)
                                                                                                                                aTTGAATATGATAAAGGCATCATAATGATTGTAATGCAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCTGCAAATTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGt         <  1:65458/140‑1 (MQ=255)
                                                                                                                                aTTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCAATGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACATACCACCACTAAGGCTGGTGTCGTAACTGCACTTGt         <  2:246920/140‑1 (MQ=255)
                                                                                                                                  tggatatTATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTCCACAAATACAACCTTTTTTCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTa       <  1:514593/137‑1 (MQ=255)
                                                                                                                                     atatGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTacaa    <  1:683087/140‑1 (MQ=255)
                                                                                                                                       atGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCCCTTGTTccaaca  >  2:151436/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                         
GAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTCTTTATGTTTTTTAGAGATAAGCATTGTAATTGGATAGAAAATAAATCCACACGCAATACCAGTTGCAATTGAATATGATAAAGGCATCATAATGATTGTAATGAAGGCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACA  >  USA300TCH1516_ALE/2372449‑2372723

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