BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	638,614	G→T	T14K (ACA→AAA)	pdxK ←	Pyridoxine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	638,614	0	G	T	99.2%	397.4 / NA	120	T14K (ACA→AAA)	pdxK	Pyridoxine kinase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base T (53/66); total (53/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

TCAAATACATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATGAATCATT  >  USA300TCH1516_ALE/638470‑638758

tCAAATACATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTttg‑t                                                                                                                                               <  1:178228/148‑1 (MQ=255)

   aaTACATCCACTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAg                                                                                                                                            <  1:408538/148‑1 (MQ=255)

      aCATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCacgac                                                                                                                                         >  2:298686/1‑144 (MQ=255)

       cATCCATTGTTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAAcc                                                                                                                                        <  2:753069/148‑1 (MQ=255)

        ctCCTTTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCg                                                                                                                                       <  1:963858/147‑1 (MQ=255)

        aTCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCg                                                                                                                                       <  1:762469/148‑1 (MQ=255)

         tCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGGACCgg                                                                                                                                      >  1:397849/1‑148 (MQ=255)

          ccATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCgcc                                                                                                                                     >  2:204322/1‑146 (MQ=255)

          ccATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGc                                                                                                                                     >  1:148596/1‑148 (MQ=255)

           ccttGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCa                                                                                                                                    <  2:966698/146‑1 (MQ=255)

           cATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCa                                                                                                                                    >  1:609629/1‑148 (MQ=255)

             ttGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAAt                                                                                                                                  <  2:181440/148‑1 (MQ=255)

              tGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTCTTG‑TCAGAACCGGCAAtt                                                                                                                                 >  1:833541/1‑148 (MQ=255)

                gTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTGGGATCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGt                                                                                                                               <  1:659655/148‑1 (MQ=255)

                 tAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGtt                                                                                                                              >  2:325147/1‑148 (MQ=255)

                  aaTGTTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTa                                                                                                                             <  2:131320/148‑1 (MQ=255)

                  aaTGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTACTTCTGGGGACGTTTTCAAATCTGCTTGCATACCTGCGCCCGCACATTTA‑TCAGAACCGACAATTGTTa                                                                                                                             >  1:129748/1‑148 (MQ=255)

                     ggTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTaaaa                                                                                                                          >  1:186053/1‑148 (MQ=255)

                     ggTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTaaaa                                                                                                                          >  1:910710/1‑148 (MQ=255)

                        gtAACCTCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATCCCTGCGCCAGCACTTTTG‑TCAGCACCGGCAATTTTTAAAActt                                                                                                                       <  2:856960/148‑1 (MQ=255)

                        gtAACATCGTGTGACCATGTATCTTTATCCCTAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCATTTCACTGG‑TCAGAACCGGCAATTGTTAAAActt                                                                                                                       <  2:955770/148‑1 (MQ=255)

                        gtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAActt                                                                                                                       <  1:67126/148‑1 (MQ=255)

                        gtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAActt                                                                                                                       <  2:134886/148‑1 (MQ=255)

                         tAACCTCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAActtt                                                                                                                      <  2:272583/148‑1 (MQ=255)

                          aaCATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAActttc                                                                                                                     <  1:407283/148‑1 (MQ=255)

                          aaCATCGTCTGACCATGTATGTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAActttc                                                                                                                     <  2:973455/148‑1 (MQ=255)

                           acacCGCGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAActttct                                                                                                                    <  2:202844/148‑1 (MQ=255)

                           aCATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCTGGA‑GTGTGTTCAGAACCGGCAATTGTTAAAActttct                                                                                                                    <  1:729457/148‑1 (MQ=255)

                           aCATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAActttct                                                                                                                    <  2:189140/148‑1 (MQ=255)

                           aCATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAActttct                                                                                                                    <  2:345074/148‑1 (MQ=255)

                           aCATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTACGGTCATTCCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTGAAActttct                                                                                                                    >  1:276567/1‑148 (MQ=255)

                              tCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTa                                                                                                                 <  1:171152/148‑1 (MQ=255)

                               cGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTaa                                                                                                                >  1:957828/1‑148 (MQ=255)

                                  ttgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGc                                                                                                             <  2:471498/147‑1 (MQ=255)

                                      ccATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGGATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCAtt                                                                                                         >  2:844358/1‑148 (MQ=255)

                                        aTGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGa                                                                                                       <  2:897350/148‑1 (MQ=255)

                                            aTCTTTATCCATACTCACGCTGGCAGTTAAAGCGACCATGCCATACGTATCTAACTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGttt                                                                                                   <  1:74081/148‑1 (MQ=255)

                                                ttATCCACAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCCTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACt                                                                                               <  1:267481/148‑1 (MQ=255)

                                                 tATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTc                                                                                              <  1:686540/148‑1 (MQ=255)

                                                   tcaaTAGTACCGATGGCAGTTAATGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCt                                                                                            <  2:537654/145‑1 (MQ=255)

                                                    ccATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTa                                                                                           <  2:344808/148‑1 (MQ=255)

                                                    ccATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTa                                                                                           >  2:390551/1‑148 (MQ=255)

                                                    ccATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTa                                                                                           >  1:210178/1‑148 (MQ=255)

                                                     cATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAACTGTTACTAATGTCATTAAAGCCCTTGCGCCTTACTCCTAc                                                                                          >  1:469539/1‑148 (MQ=255)

                                                       caGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACat                                                                                        <  2:737744/147‑1 (MQ=255)

                                                        aGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCAAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATGGTTAAAACGTTCTTTAAAGCCATTAGGCTGCACTCCTACata                                                                                       >  1:248995/1‑148 (MQ=255)

                                                        aGTAACGAGGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGTGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACata                                                                                       <  2:386100/148‑1 (MQ=255)

                                                         gTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataa                                                                                      >  2:110014/1‑148 (MQ=255)

                                                         gTAACGATGGCAGTGAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataa                                                                                      >  2:482890/1‑148 (MQ=255)

                                                            aCGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                                   >  2:382556/1‑148 (MQ=255)

                                                             cGATGTCAGTTAAAGCGGCCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAACCCATTGAGCTTCACTCCTACataataat                                                                                  <  1:196939/148‑1 (MQ=255)

                                                             cGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTAaataataat                                                                                  >  2:576903/1‑148 (MQ=255)

                                                             cGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataat                                                                                  >  2:100007/1‑148 (MQ=255)

                                                               aTGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAatat                                                                                >  1:914926/1‑148 (MQ=255)

                                                               aTGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAatat                                                                                >  1:524233/1‑148 (MQ=255)

                                                                tGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATAtt                                                                               >  1:548110/1‑148 (MQ=255)

                                                                       tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtca                                                                        <  1:466315/148‑1 (MQ=255)

                                                                        aaaGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcat                                                                       <  1:99097/148‑1 (MQ=255)

                                                                          aGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcatca                                                                     >  2:586084/1‑148 (MQ=255)

                                                                           gCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTcatcat                                                                    >  2:537268/1‑148 (MQ=255)

                                                                            cGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑GCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcata                                                                   <  2:400099/148‑1 (MQ=255)

                                                                             gACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTGAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatat                                                                  >  1:338223/1‑148 (MQ=255)

                                                                                           tATCTAATTCCTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTGCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTCAcatca                                                    >  2:738030/1‑144 (MQ=255)

                                                                                              cTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTg                                                 <  2:96512/148‑1 (MQ=255)

                                                                                                   tCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaaaaa                                            <  1:265297/148‑1 (MQ=255)

                                                                                                     ttGGAACGTTTTCAAAACTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAata                                          >  2:393232/1‑148 (MQ=255)

                                                                                                       tgAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGACCCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATAtt                                        <  1:146956/147‑1 (MQ=255)

                                                                                                       ggAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATAtt                                        <  2:914926/148‑1 (MQ=255)

                                                                                                       ggAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATAtt                                        >  1:352998/1‑148 (MQ=255)

                                                                                                       ggAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATAtt                                        <  2:524233/148‑1 (MQ=255)

                                                                                                       ggAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATAtt                                        <  1:298686/148‑1 (MQ=255)

                                                                                                        gAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTa                                       >  1:692465/1‑148 (MQ=255)

                                                                                                         aaCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTaa                                      >  1:830204/1‑148 (MQ=255)

                                                                                                         aaCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTaa                                      <  2:960471/148‑1 (MQ=255)

                                                                                                          aCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGCGCATCGAACTAACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAg                                     >  1:229710/1‑148 (MQ=255)

                                                                                                          aCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATCTCATTTTTAACCTAATTGAAAAATATTAAg                                     >  1:376944/1‑148 (MQ=255)

                                                                                                              tttCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCAtt                                 <  2:129748/148‑1 (MQ=255)

                                                                                                                tCAAATCTGCTTGCATACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCa                               >  1:897958/1‑148 (MQ=255)

                                                                                                                              aTACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttg                 <  1:541462/148‑1 (MQ=255)

                                                                                                                              aTACCTGCGCCAGCACCTTTG‑CCAGAACCGGAAATTGTTAAACCTTTCTTTAAAGCCATTGAGCTTCACTTCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttg                 <  1:325147/148‑1 (MQ=255)

                                                                                                                               tACCTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt                <  2:148596/148‑1 (MQ=255)

                                                                                                                               tACCTGCGCAAGATCTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATACTATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt                <  1:204322/148‑1 (MQ=255)

                                                                                                                               tACCCGCGCCACCACTTTTT‑CCATATCCGGCAATTGTTAAAACTTTTTTTAAAGCAACTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt                <  2:46989/148‑1 (MQ=255)

                                                                                                                                tccTGCGCCAGCACCTTTG‑TCACAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgtt               <  2:62922/147‑1 (MQ=255)

                                                                                                                                aCCTGCGCCAGCGCGTTTG‑TCAGAGCCGGCCATGGCTAACACTTCCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgtt               <  1:110014/148‑1 (MQ=255)

                                                                                                                                  cTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGa             >  1:22866/1‑148 (MQ=255)

                                                                                                                                  cTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGa             >  1:159234/1‑148 (MQ=255)

                                                                                                                                  cTGCGCCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGa             >  1:921044/1‑148 (MQ=255)

                                                                                                                                    gcgcCAGCACTTTTG‑TCAGAACCGGCGATTGTTAAAACTTTCTTTAAAGCGATTGAGCTTCACGCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAAGTGAAAAATATTCAGCATTGAATTTGTGATGAGGGTGGaaa           >  1:85273/1‑148 (MQ=255)

                                                                                                                                    gcgcCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa           >  1:817918/1‑148 (MQ=255)

                                                                                                                                    gcgcCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCAGGCGGCCCACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa           <  1:721992/148‑1 (MQ=255)

                                                                                                                                    gcgcCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCAATTGAGCATCACTCCTACATAATAATATTGAATACCTCCTATAATTTTTAAACTAATTGAAAAATATCAAGCATTCAATATTTAATGATAGTTGaaa           >  1:203533/1‑148 (MQ=255)

                                                                                                                                     cgcCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGTTTGTTGaatt          >  1:10869/1‑146 (MQ=255)

                                                                                                                                      gcCAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATg         >  1:19980/1‑148 (MQ=255)

                                                                                                                                        cAGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCTTTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATGaa       <  2:830204/148‑1 (MQ=255)

                                                                                                                                         aGCACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATGAAt      >  1:755180/1‑148 (MQ=255)

                                                                                                                                          gCACTTTCG‑CCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATGAAtc     <  2:809096/148‑1 (MQ=255)

                                                                                                                                           cACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATGAAtca    <  1:181570/148‑1 (MQ=255)

                                                                                                                                           cACTTTTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATGAAtca    <  1:992026/148‑1 (MQ=255)

                                                                                                                                             cTTTTG‑GCAGAACCGGCAATTAGTAAAACGTTCTTTGAAGTCATTGAGCTGCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATTATTGTTGAAATGAAtcatt  <  2:176587/148‑1 (MQ=255)

TCAAATACATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTG‑TCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATGAATCATT  >  USA300TCH1516_ALE/638470‑638758

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