Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,979,300 G→T 6.3% A89D (GCT→GAT)  USA300HOU_RS09860 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,979,3000GT 6.3% 236.9 / 2.8 79A89D (GCT→GAT) USA300HOU_RS09860hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (44/29);  new base T (2/3);  total (46/33)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CTCGGACTCTGTCGGCAGTCTATTCCTCTATTTAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTGAGTTTCACGGAACTCGT  >  USA300TCH1516_ALE/1979161‑1979438
                                                                                                                                           |                                                                                                                                          
ctcGGACTCTGTCGGCAGTCTATTCCTCTATTTAAAGGATTGAAAAATATTATTTTAATATTAGGCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATAATTGACCCATTTTTTGTTCAGc                                                                                                                                           >  2:607716/1‑141 (MQ=255)
 tcGGACTCTGTCGGCAGTCTATTCCTCTATTTAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGca                                                                                                                                          >  1:455524/1‑141 (MQ=255)
       tctGTCGGCAGTCTATTCCTCTATTTAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGctttcatt                                                                                                                                    >  2:11804/1‑141 (MQ=255)
               cAGTCTATTCCTCTATTTAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCaga                                                                                                                            <  1:313213/141‑1 (MQ=255)
               cAGTCTATTCCTCTATTTAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCCTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCATCATTCATTAATGCaga                                                                                                                            <  1:388724/141‑1 (MQ=255)
                aggcTATTCCTCTATTTAAAGGATTGAAAAATATTATTTCACTATTAGGCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCagag                                                                                                                           <  1:128829/138‑1 (MQ=255)
                aGTCTATTCCTCTATTTAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCagag                                                                                                                           <  1:691915/141‑1 (MQ=255)
                       tCCTCTATTTAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTATTGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCATCATTCATTAATGCAGAGATGtttt                                                                                                                    <  2:498376/141‑1 (MQ=255)
                       tCCTCAATTAAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGtttt                                                                                                                    <  2:601108/141‑1 (MQ=255)
                          tctATTTAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTATTCTTCATTCATTACTGCAGAGATGTTTTCAt                                                                                                                 <  2:635506/141‑1 (MQ=255)
                              tttAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCttt                                                                                                             >  1:575022/1‑141 (MQ=255)
                               ttAACGGATTGAAAAATATTATTTTAATAACAGTCAGCGTAACGCTCGTCTAATGGATTACCGGTAATTTGGTTGATCTGTTGGAATACTTGACGCCTTGTGTGTTCAGCATTATTTAATGCAGTAGTCTTTTTATCtttt                                                                                                            >  2:294706/1‑141 (MQ=255)
                               ttAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTTTTTTGTATTCATTAATGCAGAGATGTTTTCATCtttt                                                                                                            >  2:442100/1‑141 (MQ=255)
                                tAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCttttt                                                                                                           >  1:627039/1‑141 (MQ=255)
                                tAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCttttt                                                                                                           >  2:472072/1‑141 (MQ=255)
                                    ggATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAAt                                                                                                       >  1:378311/1‑141 (MQ=255)
                                          aaaaaTATTATTTTATTATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAAttgcttg                                                                                                 <  1:11781/141‑1 (MQ=255)
                                          aaaaaTATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAAttgcttg                                                                                                 <  1:556964/141‑1 (MQ=255)
                                          aaaaaTATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAAttgcttg                                                                                                 <  2:116561/141‑1 (MQ=255)
                                           aaaaTATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTTTTCAGCATTCATTAATGAAAAGATCTACTCATCTTTTTCAAttgcttgc                                                                                                >  2:344941/1‑141 (MQ=255)
                                           aaaaTATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAAttgcttgc                                                                                                >  2:47330/1‑141 (MQ=255)
                                                 ttattTTAACATTAGTCAGTGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCCTCTTTTTCAATTGCTTGCGCTTGt                                                                                          <  1:443442/141‑1 (MQ=255)
                                                   attTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTc                                                                                        >  2:357447/1‑141 (MQ=255)
                                                   attTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTc                                                                                        >  1:192713/1‑141 (MQ=255)
                                                   attTTAATACTAGTCAGCGTAAATTTCGTCTAATGGTTGAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCATGCGCTTGtta                                                                                        <  2:251860/141‑2 (MQ=255)
                                                      ttAGTATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTg                                                                                     <  2:108616/141‑1 (MQ=255)
                                                      ttAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTg                                                                                     <  1:173364/141‑1 (MQ=255)
                                                        aaTATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTTTTTTTTTTTTTTTTATGTAGAGTTGTTTTGATTTTTTTCAATTGCTTTCGCTTGtggcggcg                                                                                   >  2:134687/1‑134 (MQ=255)
                                                          tatTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCCTGCGCt                                                                                 <  2:438474/141‑1 (MQ=255)
                                                                     gTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCGTTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTCTTCTTGCGGttttttgcactc                                                                      >  1:252904/1‑134 (MQ=255)
                                                                     gTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAAtc                                                                      >  2:184075/1‑141 (MQ=255)
                                                                     gTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATAATTGAATCATTTTTTGTTCAGCATTCATTAATGAAGAGATGTTTTCATCTTTTTCAATTGCTTGAGCTTGTTCTTGCGCTTTTTGTAAAtt                                                                      >  2:453355/1‑140 (MQ=255)
                                                                       aaaTTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAAtctt                                                                    >  1:912/1‑141 (MQ=255)
                                                                        aaTTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTCAAtcttc                                                                   >  2:263869/1‑141 (MQ=255)
                                                                        aaTTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAAtcttc                                                                   >  2:591814/1‑141 (MQ=255)
                                                                         aTTTCGTCTAATTGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAAtcttct                                                                  <  2:42702/141‑1 (MQ=255)
                                                                         aTTTCGTCTAATTGTTTAACGATAATGTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAAtcttct                                                                  <  2:125627/141‑1 (MQ=255)
                                                                           tttctCTAATGGTTGAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCttcttc                                                                <  1:632938/137‑1 (MQ=255)
                                                                                  aaTGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAAtt                                                         >  1:11278/1‑141 (MQ=255)
                                                                                   aTGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAAttt                                                        >  1:75288/1‑141 (MQ=255)
                                                                                   aTGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAAttt                                                        >  1:252305/1‑141 (MQ=255)
                                                                                     ggTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATttct                                                      <  1:111459/141‑1 (MQ=255)
                                                                                     ggTCTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATttct                                                      <  1:351142/141‑1 (MQ=255)
                                                                                        ttAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTTTTTTTTTTTCATTAATGCAGAGCTGTTTTCATTTTTTTCAATTGCTTGCGGTTGTTCTTGGGCTTTTTTTTCTTCTTCTTCAGCATTtttttct                                                   >  2:536469/1‑140 (MQ=255)
                                                                                        ttAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCACTAATGCAGAGATGTTTTCATCTTTTTCTATTTCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCATCAATttcttta                                                   >  2:190766/1‑139 (MQ=255)
                                                                                         tAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTTCGCTTGTTCTTGCGCTTTTTGTACATCTTCTTCAGCTATTTCTTCAc                                                  >  2:291932/1‑141 (MQ=255)
                                                                                           aCGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTTTTTTTTTTTTTCTTATTAGTAGGCGTTTTGTAAATTTTCTTCAGCCTTCTCTTCACCt                                                >  2:96140/1‑141 (MQ=255)
                                                                                            cGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTTAATCTTCTTCAGCAATTTTTTCTCCtt                                               >  2:332427/1‑141 (MQ=255)
                                                                                            cGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCtt                                               >  1:685192/1‑141 (MQ=255)
                                                                                            cGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCtt                                               >  2:230862/1‑141 (MQ=255)
                                                                                             gATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTg                                              >  1:117278/1‑141 (MQ=255)
                                                                                             gATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTg                                              <  2:351484/141‑1 (MQ=255)
                                                                                               tAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGct                                            >  1:39024/1‑140 (MQ=255)
                                                                                               tAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCa                                            >  1:692804/1‑141 (MQ=255)
                                                                                                    tGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTCCTTGCGGTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCAtttgt                                       >  1:160059/1‑141 (MQ=255)
                                                                                                      gTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCAtttgttt                                     <  2:34564/141‑1 (MQ=255)
                                                                                                         taTTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCAtttgtttttg                                  <  2:613652/140‑1 (MQ=255)
                                                                                                              cTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGa                             >  1:580063/1‑141 (MQ=255)
                                                                                                               ttGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGaa                            >  1:600447/1‑141 (MQ=255)
                                                                                                               ttGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGaa                            >  1:383650/1‑141 (MQ=255)
                                                                                                                 ggAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGCAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGATGt                          >  1:201783/1‑141 (MQ=255)
                                                                                                                     tACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAAt                      >  1:238122/1‑141 (MQ=255)
                                                                                                                      ttctGACTCATTTTTTGGTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAAtt                     <  1:181522/138‑1 (MQ=255)
                                                                                                                      aCTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTCTTCATTTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAAtt                     <  1:535967/141‑1 (MQ=255)
                                                                                                                      aCTTGACTCATTTTTTGTTCAACATTCAGTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAAtt                     <  1:525881/141‑1 (MQ=255)
                                                                                                                        ttGACTCATTTTTTTTTCAGCATTCATGAATGCAGAGATGTTTTCAGCTTTTTCAATTGCTTGCGCTTTTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTg                   <  1:142169/141‑1 (MQ=255)
                                                                                                                          gACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTgcg                 >  2:117881/1‑139 (MQ=255)
                                                                                                                          gACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTGAg                 >  1:374126/1‑141 (MQ=255)
                                                                                                                          gACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTGAg                 >  1:218816/1‑141 (MQ=255)
                                                                                                                                 tttttGTTCAGCATTCATTAATTCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTGAGTTTCACg          <  2:564084/141‑1 (MQ=255)
                                                                                                                                  ttttGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTTGTTTTTCTTATATGACTGCTTCACCATGCGCCTTTTTTTGTTGGAAGTTAATTTGAGTTTCACgg         <  2:146797/141‑1 (MQ=255)
                                                                                                                                  ttttGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTGAGTTTCACgg         <  2:624265/141‑1 (MQ=255)
                                                                                                                                  ttttGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTGAGTTTCACgg         >  1:323410/1‑141 (MQ=255)
                                                                                                                                    ttGTTCAGCATTCTTTAATTCAGAGTTGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTCTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTGAGTTTCACGGaa       <  2:252904/141‑1 (MQ=255)
                                                                                                                                    ttGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTCGCATTTGTTTTTGTTGGAAGTTAATTTGAGTTTCACGGaa       <  2:114969/141‑1 (MQ=255)
                                                                                                                                     tGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTGAGTTTCACGGAAc      <  2:192713/141‑1 (MQ=255)
                                                                                                                                      gTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTGAGTTTCACGGAACt     >  1:223854/1‑141 (MQ=255)
                                                                                                                                      gTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTGAGTTTCACGGAACt     >  1:297244/1‑141 (MQ=255)
                                                                                                                                         cAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTGAGTTTCACGGAACtcgt  >  1:46078/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                          
CTCGGACTCTGTCGGCAGTCTATTCCTCTATTTAAAGGATTGAAAAATATTATTTTAATATTAGTCAGCGTAAATTTCGTCTAATGGTTTAACGATAATTTGGTTGATTTCTTGGAATACTTGACTCATTTTTTGTTCAGCATTCATTAATGCAGAGATGTTTTCATCTTTTTCAATTGCTTGCGCTTGTTCTTGCGCTTTTTGTAAATCTTCTTCAGCAATTTCTTCACCTTGCATTTGTTTTTGTTGGAAGTTAATTTGAGTTTCACGGAACTCGT  >  USA300TCH1516_ALE/1979161‑1979438

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