Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,231,446 C→T 6.6% Q31* (CAA→TAA)  stp → Serine/threonine phosphatase stp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,231,4460CT 6.6% 202.3 / 2.2 69Q31* (CAA→TAA) stpSerine/threonine phosphatase stp
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base T (2/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGT  >  USA300TCH1516_ALE/1231308‑1231583
                                                                                                                                          |                                                                                                                                         
ggTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTGTTTATAATCAAACTAATcaacca                                                                                                                                         >  1:118181/1‑139 (MQ=255)
   cAATTAAGAGCAAAGGAACGACAAGTAGGAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACtt                                                                                                                                      <  2:658791/141‑1 (MQ=255)
                aGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCtgtgtg                                                                                                                         <  1:863496/141‑1 (MQ=255)
                    aCGACAAGTAGATACGAGGTAAAGACAAATGCTAGAGGCACAAGTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATACTCAAACTAATCAACAGCCTTCAGTTCTGCGTGAtag                                                                                                                     <  2:790344/141‑3 (MQ=255)
                    aCGACAAGTAGAAACGAGGTAAAGACAAATACTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAACCACCCTAATCATCAACCTTTAGTCCTGTGTGATgg                                                                                                                     <  2:85848/141‑1 (MQ=255)
                      aaCAAGAAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTa                                                                                                                   <  2:232004/140‑1 (MQ=255)
                       aggaGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAATAACTTTTAGTTCTGTGTGATGGTAt                                                                                                                  <  2:231647/138‑1 (MQ=255)
                       aaaaGTAGAAACCGGGTAAAGACAAATGCTAGAGGCACAATTTTTTACGGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTTGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTAt                                                                                                                  <  2:31981/139‑1 (MQ=255)
                       aCAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTAt                                                                                                                  >  1:208187/1‑141 (MQ=255)
                       aCAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATAAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTAt                                                                                                                  >  2:116114/1‑141 (MQ=255)
                       aCAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGAAACTGTACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTAt                                                                                                                  <  1:26727/141‑1 (MQ=255)
                        cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAATAACTTTTAGTTCTGTGTGATGGTATg                                                                                                                 <  2:870968/141‑1 (MQ=255)
                        cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATg                                                                                                                 >  2:282117/1‑141 (MQ=255)
                        cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATg                                                                                                                 >  1:392580/1‑141 (MQ=255)
                        cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCCGTGTGATGGTATg                                                                                                                 >  1:766846/1‑141 (MQ=255)
                        cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAAATTTTAGTTCTGTGTGATGGTATg                                                                                                                 >  2:425447/1‑141 (MQ=255)
                        cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAAAAGAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATg                                                                                                                 >  2:287129/1‑141 (MQ=255)
                        cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATGTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGAAGCGGGTGGTATTTTTTATAATCAAACTAATAANCAANTNTTAGTCTTGTGTGATGGTATg                                                                                                                 >  1:486228/1‑141 (MQ=255)
                        cAAGTAGAAACGAGCTAAAGACAAATGCTAGCGGCCCAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAATAACTTTTAGTTCTGTGTGATGGTATg                                                                                                                 <  2:867881/141‑1 (MQ=255)
                           gTAGAAACGGGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGt                                                                                                              <  2:143985/141‑1 (MQ=255)
                           gTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGt                                                                                                              <  2:567734/141‑1 (MQ=255)
                             aGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTgg                                                                                                            <  1:781462/141‑1 (MQ=255)
                              gAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGc                                                                                                           >  1:552862/1‑141 (MQ=255)
                              gAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGc                                                                                                           >  1:372698/1‑141 (MQ=255)
                              gAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTGTTAGTTCTGTGTGATGGTATGGGTGGc                                                                                                           >  1:403299/1‑141 (MQ=255)
                                     ggTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGATATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAg                                                                                                    >  1:202900/1‑141 (MQ=255)
                                     ggTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGAATGGGTGGCCATaaac                                                                                                    >  2:284948/1‑140 (MQ=255)
                                      gTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGc                                                                                                   <  1:322265/141‑1 (MQ=255)
                                        aaaGACAAATGCTATAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTAGCCATAAAGCAg                                                                                                 <  1:606125/141‑1 (MQ=255)
                                        aaaGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAg                                                                                                 >  1:43822/1‑141 (MQ=255)
                                           gACAAATGCTAGCGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGgag                                                                                              <  2:479314/141‑1 (MQ=255)
                                           gACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGgag                                                                                              <  2:641979/141‑1 (MQ=255)
                                                   cTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCa                                                                                      >  2:470547/1‑141 (MQ=255)
                                                   cTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCa                                                                                      <  2:678433/141‑1 (MQ=255)
                                                   cTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTAGGGGTGGCCATAAATCAGGAGAAGTTGcg                                                                                      >  1:195954/1‑140 (MQ=255)
                                                     agagGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAg                                                                                    <  2:43822/141‑1 (MQ=255)
                                                          cacaATTTTTTACTGATACTGGACAACATAGAGATCAGACTGACGAAAAGGGCGGTATTCTTTCTACTCAAAATAATAGAAGATAATCACGTTTGCGCAGTGGCACGGGTGCCCATAAAGCATGATAAGGTGCAAGTAAAt                                                                               >  2:315396/1‑141 (MQ=255)
                                                          cacaATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAAt                                                                               >  1:79735/1‑141 (MQ=255)
                                                            caATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAAttt                                                                             >  1:177733/1‑141 (MQ=255)
                                                                 ttttACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCGTAAAGCAGGAGAAGTTGCAAGTAAATTTTTTAc                                                                        >  1:436399/1‑141 (MQ=255)
                                                                   ttACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAAGCAAACTAATCAACAACTTGTAGTTCTGTGTGATGGTAGGGGTGGCCATAAAGCGGGAGAAGTTGCCAGTAAATTTGTTACAg                                                                      >  1:902276/1‑141 (MQ=255)
                                                                     aCTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGATATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGAt                                                                    <  2:675182/141‑1 (MQ=255)
                                                                      cTGATACTGGACAACATAGAGATAAGCATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTACATAAACGTAAAGGTATGGGAGGCCAGAAAGCAGGAGAAGTTACAAGTAAATTTGTTACAGATg                                                                   >  1:926276/1‑141 (MQ=255)
                                                                         aTACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAAATAATCAACAACTGTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGt                                                                <  2:335669/141‑1 (MQ=255)
                                                                          tgcTGTACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTATAATCAAACTAATGAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTg                                                              <  2:231848/139‑1 (MQ=255)
                                                                                acaacaTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCCGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGACATc                                                         >  1:206823/1‑141 (MQ=255)
                                                                                acaacaTAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATc                                                         >  1:41336/1‑141 (MQ=255)
                                                                                  aacaTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGTGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATccc                                                       >  1:264043/1‑141 (MQ=255)
                                                                                     aTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTTATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGtt                                                    >  1:151463/1‑141 (MQ=255)
                                                                                     aTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTATTAGTTCTGTGTTATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTCAATTTGTTACAGATGAGTTGAAATCCCgct                                                    >  1:724845/1‑139 (MQ=255)
                                                                                        gagaTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAATAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTg                                                 >  1:300826/1‑141 (MQ=255)
                                                                                        gagaTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAATAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTg                                                 >  1:74564/1‑141 (MQ=255)
                                                                                        gagaTAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGCGAGGTATTAAAGCAGGATAAGTTGCAAGAAAACTTGTTACAGATGCGTTGAAATCCCGTTTTg                                                 <  2:429967/141‑1 (MQ=255)
                                                                                           agaaGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAACCACCAACTTTTAGTTGTGTGTGATGGTATGGGTGGCCAGAAAGCAGGAGAAGTTGCAAGTAAATTTGGTACAGAGGAGTTGAGATCCCTTGTTTAAg                                              <  2:309154/139‑1 (MQ=255)
                                                                                           aTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAg                                              <  2:448847/141‑1 (MQ=255)
                                                                                           aTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAg                                              <  1:882744/141‑1 (MQ=255)
                                                                                            tAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGc                                             <  1:331907/141‑1 (MQ=255)
                                                                                             aaGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCg                                            <  1:776831/141‑1 (MQ=255)
                                                                                             aaGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCg                                            <  2:177733/141‑1 (MQ=255)
                                                                                                aaTAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGACGCGGaa                                         >  2:124701/1‑141 (MQ=255)
                                                                                                 aTGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACGTTTAGTTACGTGTGATGGTATGGGTGACCATACAGCAGGAGAAGTTGCAAGTAAAGTTATTACAGATGAGTGGAAATCGCGTATTGAGGCGGaaa                                        >  2:411099/1‑141 (MQ=255)
                                                                                                     aGATGCGGGTGGTATTTTTTATAATCAAACTAGTCAACAACTTTTATGTTTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTATAAATCCCGTTTTGAAGCGGAAAATCt                                    <  1:225441/141‑1 (MQ=255)
                                                                                                          cGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAg                               <  2:401254/141‑1 (MQ=255)
                                                                                                          cGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAg                               <  2:2955/141‑1 (MQ=255)
                                                                                                               ggtATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGaacaa                          >  1:671937/1‑141 (MQ=255)
                                                                                                                  aTTTTTTATAATCAAACTAATCAACAACTTTTAGTTCCGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAACAACTTATACAACAACAt                       >  2:649489/1‑141 (MQ=255)
                                                                                                                       ttATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGGTTTGGGGGGGAAAATCTTATGGAACAAAATcgcga                  >  1:18913/1‑137 (MQ=255)
                                                                                                                                     aTCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAGGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGc    >  2:593227/1‑141 (MQ=255)
                                                                                                                                     aTCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGc    >  1:717126/1‑141 (MQ=255)
                                                                                                                                     aTCAACAACTTTTAGTTCCGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTTGGTATTGGTTGc    <  2:206823/141‑1 (MQ=255)
                                                                                                                                       caacaaCTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGt  <  2:552862/141‑1 (MQ=255)
                                                                                                                                       caacaaCTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGt  <  2:372698/141‑1 (MQ=255)
                                                                                                                                       caacaaCTTTTAGTTCGGTGTGATTGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGt  <  2:403299/141‑1 (MQ=255)
                                                                                                                                       caacaaCTTTTAGTTCCGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGGAAATTGGTTGCGt  >  2:776320/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                         
GGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGT  >  USA300TCH1516_ALE/1231308‑1231583

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: