Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 2,856,087 T→A P80P (CCT→CCA rpoA_2 → DNA‑directed RNA polymerase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,856,0870TA98.0% 346.9 / ‑1.8 98P80P (CCT→CCArpoA_2DNA‑directed RNA polymerase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (54/42);  total (56/42)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GGAGATGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAAT  >  USA300TCH1516_ALE/2855949‑2856221
                                                                                                                                          |                                                                                                                                      
ggAGATGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCaaa                                                                                                                                      >  4:119252/1‑141 (MQ=255)
  agaTGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAAc                                                                                                                                    >  2:81912/1‑141 (MQ=255)
   gaTGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCACCTACATTTTAAAACGGAAGTGCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACg                                                                                                                                   <  1:261161/141‑1 (MQ=255)
     tGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTc                                                                                                                                 >  3:99662/1‑141 (MQ=255)
     tGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTc                                                                                                                                 >  3:99656/1‑141 (MQ=255)
           gCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACaaa                                                                                                                           <  2:508557/141‑1 (MQ=255)
            cGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAAt                                                                                                                          <  1:121509/141‑1 (MQ=255)
               tGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGAt                                                                                                                       >  1:486393/1‑141 (MQ=255)
                  tCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGa                                                                                                                    <  2:444529/141‑1 (MQ=255)
                   ccTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAt                                                                                                                   <  1:113445/141‑1 (MQ=255)
                   ccTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAt                                                                                                                   >  1:306240/1‑141 (MQ=255)
                     tAAAGCTATATCAATAGTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAttt                                                                                                                 <  1:236404/141‑1 (MQ=255)
                             tatCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAgttg                                                                                                         >  2:551529/1‑141 (MQ=255)
                              atCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGTAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAgttgt                                                                                                        >  3:47127/1‑141 (MQ=255)
                              atCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAgttgt                                                                                                        >  3:189191/1‑141 (MQ=255)
                              atCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATAATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTAAACAAATGATAGATTTTATAgttgt                                                                                                        >  2:348514/1‑141 (MQ=255)
                               tCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAgttgtt                                                                                                       >  1:310154/1‑141 (MQ=255)
                                    atTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTAcagc                                                                                                  <  3:5678/141‑1 (MQ=255)
                                     ttGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATTCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTCTTAcagca                                                                                                 >  4:244169/1‑141 (MQ=255)
                                     ttGGAACAAGCTTTATTTCAGCACCAACTACAGTTTAAAACGGAAGTTCAGTCATCATGACCGTGTAAGATAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGACAACAAAAGATAGATTTTATAGTTGTTAcagca                                                                                                 >  1:630440/1‑141 (MQ=255)
                                         aaCAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCAt                                                                                             >  2:683147/1‑141 (MQ=255)
                                         aaCAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCAt                                                                                             >  1:565614/1‑141 (MQ=255)
                                            aaGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAg                                                                                          >  4:297307/1‑141 (MQ=255)
                                            aaGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAg                                                                                          >  2:43141/1‑141 (MQ=255)
                                               cTTTATTTCAGCACTAACTACATTTTAAAACGGAAGTTCAGTCATCATTACTTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGata                                                                                       <  1:259318/141‑1 (MQ=255)
                                               cTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGata                                                                                       <  1:167919/141‑1 (MQ=255)
                                                 ttATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATAtt                                                                                     >  2:539970/1‑141 (MQ=255)
                                                  tATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACCGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTg                                                                                    <  2:404089/141‑1 (MQ=255)
                                                     ttCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAat                                                                                 >  4:157319/1‑141 (MQ=255)
                                                       cAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAatta                                                                               >  2:407922/1‑141 (MQ=255)
                                                       cAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTa                                                                                                      >  2:150418/1‑118 (MQ=255)
                                                       cAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTa                                                                                                      <  1:150418/118‑1 (MQ=255)
                                                       cAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAatta                                                                               >  1:657670/1‑141 (MQ=255)
                                                        aGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCAtt                                                                                            >  3:121442/1‑127 (MQ=255)
                                                        aGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCAtt                                                                                            <  4:121442/127‑1 (MQ=255)
                                                        aGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAattat                                                                              >  2:653696/1‑141 (MQ=255)
                                                         gCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAattatt                                                                             >  3:129401/1‑141 (MQ=255)
                                                         gCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAattatt                                                                             <  2:539381/141‑1 (MQ=255)
                                                          cACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACg                                                                                                                                   >  2:130828/1‑86 (MQ=255)
                                                          cACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACg                                                                                                                                   <  1:130828/86‑1 (MQ=255)
                                                              aaCTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGtt                                                                        <  1:349572/141‑1 (MQ=255)
                                                               aCTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTc                                                                       >  3:100401/1‑141 (MQ=255)
                                                                 tACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCAc                                                                     <  4:248030/141‑1 (MQ=255)
                                                                 tACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCAc                                                                     <  3:89275/141‑1 (MQ=255)
                                                                 tACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCAc                                                                     >  4:209592/1‑141 (MQ=255)
                                                                 aaagggttAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCAc                                                                     <  1:136226/135‑1 (MQ=255)
                                                                      tttAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGtt                                                                >  1:139546/1‑141 (MQ=255)
                                                                         aaaaCGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTaca                                                             <  4:207826/141‑1 (MQ=255)
                                                                          aaaCGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTacaa                                                            >  2:265642/1‑141 (MQ=255)
                                                                            aCGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTacaaca                                                          <  2:139546/141‑1 (MQ=255)
                                                                            aCGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTacaaca                                                          >  3:80077/1‑141 (MQ=255)
                                                                            aCGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTacaaca                                                          <  2:649056/141‑1 (MQ=255)
                                                                             cGGAAGTTCAGTCATCNTTACCTTTTAAGNTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAg                                                         >  1:502600/1‑141 (MQ=255)
                                                                             cGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAg                                                         >  1:288380/1‑141 (MQ=255)
                                                                             cGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAg                                                         <  1:655275/141‑1 (MQ=255)
                                                                              ggAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGa                                                        >  1:275153/1‑141 (MQ=255)
                                                                              ggAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGa                                                        >  2:661903/1‑141 (MQ=255)
                                                                              ggAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGa                                                        >  1:257476/1‑141 (MQ=255)
                                                                               gAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGaa                                                       >  4:217337/1‑141 (MQ=255)
                                                                                aaGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAat                                                      <  4:165177/141‑1 (MQ=255)
                                                                                aaGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAat                                                      <  3:9682/141‑1 (MQ=255)
                                                                                aaGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAat                                                      <  2:219636/141‑1 (MQ=255)
                                                                                  gTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAattt                                                    <  1:159154/141‑1 (MQ=255)
                                                                                     cAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAatttatt                                                 <  2:446585/141‑1 (MQ=255)
                                                                                       gTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTg                                               >  2:622056/1‑141 (MQ=255)
                                                                                            catTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTg                                          >  4:253883/1‑141 (MQ=255)
                                                                                               tACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTcc                                       >  3:231138/1‑141 (MQ=255)
                                                                                                aCCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTccc                                      <  1:683147/141‑1 (MQ=255)
                                                                                                 ccTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTacaac                                                           >  4:19072/1‑119 (MQ=255)
                                                                                                 ccTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTacaac                                                           <  3:19072/119‑1 (MQ=255)
                                                                                                  cTTTTAAGTTAACAGGAGATGTCTCATGCAATCATGTGCCAAATCGTCAACAAATTATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAGTTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTg                                    <  3:252923/141‑1 (MQ=255)
                                                                                                   ttttAGGTTAACAGGAGATGTGTCATGCAATTATGCGCCAAAACGTCAACAAATTATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTgg                                   <  3:244169/141‑1 (MQ=255)
                                                                                                    tttAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGc                                  >  2:194700/1‑141 (MQ=255)
                                                                                                       aaGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGc                               >  4:282207/1‑141 (MQ=255)
                                                                                                            aaCAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGa                          >  1:212893/1‑141 (MQ=255)
                                                                                                             aCAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAatttat                                                  <  2:621793/116‑1 (MQ=255)
                                                                                                             aCAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAatttat                                                  >  1:621793/1‑116 (MQ=255)
                                                                                                             aCAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGAc                         <  2:698170/141‑1 (MQ=255)
                                                                                                             aCAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGAc                         >  4:239884/1‑141 (MQ=255)
                                                                                                               aGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACat                       >  1:256798/1‑141 (MQ=255)
                                                                                                                ggAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACata                      <  2:275080/141‑1 (MQ=255)
                                                                                                                ggAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACata                      >  1:569694/1‑141 (MQ=255)
                                                                                                                  agaTGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACAtata                    >  1:300994/1‑141 (MQ=255)
                                                                                                                   gaTGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATAc                   >  2:105780/1‑141 (MQ=255)
                                                                                                                    aTGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACg                  <  1:626111/141‑1 (MQ=255)
                                                                                                                     tgtgTCATGCAATCATGCGCCAAAACGTCAACAAATGAt                                                                                                                       <  1:296615/39‑1 (MQ=255)
                                                                                                                     tgtgTCATGCAATCATGCGCCAAAACGTCAACAAATGAt                                                                                                                       >  2:296615/1‑39 (MQ=255)
                                                                                                                        gtCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACTGAcc              >  2:261467/1‑141 (MQ=255)
                                                                                                                         tCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCt             >  2:266949/1‑141 (MQ=255)
                                                                                                                          cATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTc            <  4:129401/141‑1 (MQ=255)
                                                                                                                          cATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTc            <  2:143209/141‑1 (MQ=255)
                                                                                                                          cATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTc            <  4:78686/141‑1 (MQ=255)
                                                                                                                           aTGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATCGTTGTTATAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCa           <  1:245560/141‑1 (MQ=255)
                                                                                                                           aTGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTATTACAGCAGCATTAGATATTGAATTATTACGTTAACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATTAGGACCTCa           >  4:97294/1‑141 (MQ=255)
                                                                                                                                aTCATTCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACGACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAatat      >  3:299597/1‑141 (MQ=255)
                                                                                                                                aTCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAatat      <  2:466229/141‑1 (MQ=255)
                                                                                                                                 tCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATAtt     <  3:297307/141‑1 (MQ=255)
                                                                                                                                  cATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTa    >  2:478258/1‑141 (MQ=255)
                                                                                                                                    tGCGCCAAAACGACAACAACAGATAGATTTTACAGTTGTTACAGCAGCATTAGATAATGCATTATTACGTTCACTGGTTACAACAGAATTTATATGGGCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAAt  >  2:121908/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                      
GGAGATGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAAT  >  USA300TCH1516_ALE/2855949‑2856221

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