Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,260,448 G→T 6.7% P260H (CCT→CAT)  gmuF ← putative mannose‑6‑phosphate isomerase GmuF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,260,4480GT 6.7% 207.2 / 2.0 72P260H (CCT→CAT) gmuFputative mannose‑6‑phosphate isomerase GmuF
Reads supporting (aligned to +/- strand):  ref base G (36/31);  new base T (2/3);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

ATGTGAAATCACCTTCAAAGACACTATCCAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTGCATTGTACTGA  >  USA300TCH1516_ALE/2260308‑2260585
                                                                                                                                            |                                                                                                                                         
aTGTGAAATCACCTTCAAAGACACTATCCAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAg                                                                                                                                           <  2:595124/141‑1 (MQ=255)
  ttgAAATCCCCTTCAAAGACACGATCCAAATCTTCAGAATTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGc                                                                                                                                         <  1:643203/140‑1 (MQ=255)
  gtgAAATCACCTTCAAAGACACTATCCAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTTCATCACTGTAACTAAACAGAACTCTCTAGGc                                                                                                                                         <  1:592401/141‑1 (MQ=255)
        tgtccggcAAAGAACCTATCCAAATCTTCAGAAGTCAAAATATAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCata                                                                                                                                   <  2:495590/134‑1 (MQ=255)
        tCACCTTCAAAGACACTATCCAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAAAAGAACTCTCTAGGCTTCata                                                                                                                                   >  1:287265/1‑141 (MQ=255)
                  aGACACTATCCAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACg                                                                                                                         <  2:567425/141‑1 (MQ=255)
                   gACACTATCCAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGt                                                                                                                        <  2:546397/141‑1 (MQ=255)
                     cacTATCCAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTATGCTTCATATAATTTAACGTGc                                                                                                                      >  2:89989/1‑141 (MQ=255)
                     cacTATCCAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGc                                                                                                                      >  2:177999/1‑141 (MQ=255)
                      acTATCCAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGcc                                                                                                                     <  1:509532/141‑1 (MQ=255)
                      acTATCCAAATCTTCAGAAGCCAAAATAAAGTTTGTACCCGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGcc                                                                                                                     <  1:16705/141‑1 (MQ=255)
                         aTCCAAATCTTCAGAAGTCAAAATAAAGTTTGTACCTGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCGTCATATAATTTAACGTGCCAGa                                                                                                                  >  1:461246/1‑141 (MQ=255)
                         aTCCAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGa                                                                                                                  >  1:162832/1‑141 (MQ=255)
                          tCCAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGTACTCTCTATGCTTCATATAATTTAACGTGCAAGaa                                                                                                                 >  2:586543/1‑141 (MQ=255)
                           ccAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGaaa                                                                                                                >  2:459541/1‑141 (MQ=255)
                            cAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAAt                                                                                                               <  1:578220/141‑1 (MQ=255)
                            cAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAAt                                                                                                               <  1:192013/141‑1 (MQ=255)
                             aaaTCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAAtt                                                                                                              >  1:113199/1‑141 (MQ=255)
                               aTCTTCAGAAGTCACAATAAAGTTTGTACCAGTAGTCAGTTTGAACATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACTCTGTAACTAAACAGAACTCTCTATGCTTCATATAATTTAACGTGCCAGAAATTTc                                                                                                            <  1:595682/141‑1 (MQ=255)
                                 cTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTccc                                                                                                          >  1:78761/1‑141 (MQ=255)
                                   tCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATAGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCAt                                                                                                        >  2:605216/1‑141 (MQ=255)
                                          tCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCa                                                                                                 >  1:215276/1‑141 (MQ=255)
                                           cAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTGCCAACACTGTAACTAAACAGAGCGCCCTAGGCTTCTTATAATTTAACCCGCCAGAAATTTCCCATTTAACCaa                                                                                                >  1:252142/1‑141 (MQ=255)
                                              aaTAAAGTTTGGACCAGTCGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACTCTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGt                                                                                             <  1:217227/141‑1 (MQ=255)
                                                  aaGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAg                                                                                         <  1:604789/141‑1 (MQ=255)
                                                   aGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGa                                                                                        >  2:528630/1‑141 (MQ=255)
                                                   aGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGa                                                                                        >  2:512449/1‑141 (MQ=255)
                                                     tttGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATATCCCATTTAACCAATGTAAAGaaa                                                                                      <  1:426092/141‑1 (MQ=255)
                                                              gtagtCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGAt                                                                             >  1:649056/1‑141 (MQ=255)
                                                              gtagtCAGTTTGAAAATTTCACCAACGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGAt                                                                             >  1:423255/1‑141 (MQ=255)
                                                              gtagtCAATTTGAAAATTTCCCCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTACATATAATTTAACGTGCCAGAAATTTCCCATTTAGCCAATGTAAAGAAATCATTCGAt                                                                             >  1:419369/1‑141 (MQ=255)
                                                                ggtCAGTTTGAAAATTTCACCATCGACAATCATTCGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACTTGCCAGAAATTTCCCATGTAACCAATGTAAAGAAATCATTCGATAc                                                                           <  2:10044/140‑1 (MQ=255)
                                                                     gTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAGCACAGTAACTAAACAGAACTCTCTAGGCTTCCTCTAATTTAACGTGCCAGAAATTTCCCATTTACCCAATGTAAAGAAATCATTCGATCCAAtgg                                                                      >  2:2851/1‑140 (MQ=255)
                                                                     gTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATGATTCGATACAAtgt                                                                      >  1:317246/1‑141 (MQ=255)
                                                                           aaaTTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACACTGTGTGTaa                                                                >  1:54714/1‑140 (MQ=255)
                                                                           aaaTTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTac                                                                >  1:233207/1‑141 (MQ=255)
                                                                              tttCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTATGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACt                                                             <  2:522680/141‑1 (MQ=255)
                                                                                tCACCATCGACTATCATTTGCCCTTCGCCTTCCAACACTGTACATAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATCTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTa                                                           <  1:193961/141‑1 (MQ=255)
                                                                                tCACCATCGACAATCATTTGCCCTTCGCCTTCCGACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCACTCGATACAATGTGTGTACACTTa                                                           <  2:370732/141‑1 (MQ=255)
                                                                                tCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTACCTAACCAGAACTCTCTAGGCTTCATATGATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTa                                                           <  1:354839/141‑1 (MQ=255)
                                                                                 ccccATCTGCAATCATTTGCCCTGCGCCTTCCAACACTGTAACTAATCAGAACTCTCTAGGCTTCATAGAATGTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTAt                                                          <  1:298995/139‑1 (MQ=255)
                                                                                  aCCATCGACAATCATTTGCCCTTCGCCTACCAACACGATAACTAAACGGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTCAACTTATg                                                         >  1:421831/1‑141 (MQ=255)
                                                                                  aCCATCGACAATCATCTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATg                                                         >  1:494592/1‑141 (MQ=255)
                                                                                    cATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGt                                                       >  1:430043/1‑141 (MQ=255)
                                                                                                tttGCCCTTCGCCTTCCGATACTGCAACTAACCAGAACTCTCTAGGCTTCATATAATCTAACGTGCCAGAAATTTCCCATTGAACCAATGTAAGGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAAttt                                           <  2:348976/141‑1 (MQ=255)
                                                                                                     ccTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCtt                                      >  1:601693/1‑141 (MQ=255)
                                                                                                         cGCCTTCCCACTCTGTAACTACACAGAACTCTCTCTGCTCCGTATAATTTAACGTGCCAGAAATTTCCCATTTAATTAATGTCAAGAAATCATTCGATACAATGTGTGCACACTTATGGTTTTCAATAATTTCGCTTTCAg                                  <  2:625187/141‑1 (MQ=255)
                                                                                                           ccTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGc                                <  1:528630/141‑1 (MQ=255)
                                                                                                            cTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCa                               <  2:494592/141‑1 (MQ=255)
                                                                                                              tgcAACACTGTAACGAAACAGAACTCTCTAGGCTTCATATAAATTAACGTGCCAGAAATTTCCCATTTAACCAATGGAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCaaa                             <  2:476718/139‑1 (MQ=255)
                                                                                                              tCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCaaa                             <  2:530277/141‑1 (MQ=255)
                                                                                                                    ccTGTAACTAAACAGAACCCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTa                       <  1:371125/140‑1 (MQ=255)
                                                                                                                     ctttaactaacCAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATATAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATTTTAg                      <  2:406960/138‑1 (MQ=255)
                                                                                                                     cTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGACATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAg                      <  2:403623/141‑1 (MQ=255)
                                                                                                                       gTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGt                    <  1:223318/141‑1 (MQ=255)
                                                                                                                              aaaagtaCTCTCTAGGCTTCATATAATTTAACGTTCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTg             <  2:22291/135‑1 (MQ=255)
                                                                                                                              aaCAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTg             <  1:73042/141‑1 (MQ=255)
                                                                                                                               aCAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTGc            >  1:622815/1‑141 (MQ=255)
                                                                                                                                cAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTGCa           >  1:49039/1‑141 (MQ=255)
                                                                                                                                 aGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGAGTGTACACTTCTGGTTTTCAATAATTTCACTGTCAGGCAAAATATTAGGTAATGGTGCAt          >  2:148304/1‑141 (MQ=255)
                                                                                                                                  gAACTCTCTAGGCTTCATATAGTTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTGAGGCAAAATATTAGGTAATGGTGCAtt         >  2:285623/1‑141 (MQ=255)
                                                                                                                                  gAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTGCAtt         >  1:597088/1‑141 (MQ=255)
                                                                                                                                     ctctctAGGCTTCATATAATTTAATGTGCCAGAAATTTCCCGTTTAACCGATGTAAATAAATCATTCGATACAGTGTGTGTACACTTTTGGGTTTAGATAATTTAGCTTTCAGGCAAAATATTAAGTAATGGTGCATAGtc      >  2:285071/1‑140 (MQ=255)
                                                                                                                                     ctctctAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTGCATTGTa      <  2:543764/141‑1 (MQ=255)
                                                                                                                                     ctctctAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTGCATTGTa      >  1:98491/1‑141 (MQ=255)
                                                                                                                                     ctctctAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTACGGAAATATAAGAGGTAATGGTGCATTGTa      >  2:136687/1‑141 (MQ=255)
                                                                                                                                     ctctctAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGGGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTGCATTGTa      <  2:84965/141‑1 (MQ=255)
                                                                                                                                       ctctAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTGCATTGTACt    >  2:73254/1‑141 (MQ=255)
                                                                                                                                       ctctAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTGCATTGTACt    >  2:116415/1‑141 (MQ=255)
                                                                                                                                       ctctAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGGGTACACTTATGGTATTCACTAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTGCATTGTACt    >  2:202650/1‑141 (MQ=255)
                                                                                                                                         ctAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAACAATTAGGTGATGATGCATTGTACtta  >  2:609775/1‑139 (MQ=255)
                                                                                                                                         ctAGGCTTCATATAATTTAACGTGCCAGAAATTTCACATTTAACAAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTAGCTTTCAGGCAAAATATTAGGTACTGGTGCATTGTACTGa  >  1:231447/1‑141 (MQ=255)
                                                                                                                                         ccaGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTGCATTGTACTGa  <  2:215276/139‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                         
ATGTGAAATCACCTTCAAAGACACTATCCAAATCTTCAGAAGTCAAAATAAAGTTTGTACCAGTAGTCAGTTTGAAAATTTCACCATCGACAATCATTTGCCCTTCGCCTTCCAACACTGTAACTAAACAGAACTCTCTAGGCTTCATATAATTTAACGTGCCAGAAATTTCCCATTTAACCAATGTAAAGAAATCATTCGATACAATGTGTGTACACTTATGGTTTTCAATAATTTCGCTTTCAGGCAAAATATTAGGTAATGGTGCATTGTACTGA  >  USA300TCH1516_ALE/2260308‑2260585

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: