BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,538,398	T→G	5.9%	G168G (GGA→GGC)	aroC ←	Chorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,538,398	0	T	G	5.9%	405.9 / 14.3	136	G168G (GGA→GGC)	aroC	Chorismate synthase
Reads supporting (aligned to +/- strand): ref base T (57/71); new base G (4/4); total (61/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 141 total aligned reads displayed.

TTGAATCTCCTTCATTTTTAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGCACATTTCGTAAATCACGATGATT  >  USA300TCH1516_ALE/1538262‑1538538

ttGAATCTCCTTCATTTTTAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCAc                                                                                                                                          <  1:631281/141‑1 (MQ=255)

 tGAATCTCCTTCATTTTTAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGCTTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCAcc                                                                                                                                         >  2:249711/1‑141 (MQ=255)

    aTCTCCTTCATTTTTAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTACTTCAACCTAt                                                                                                                                      >  2:505510/1‑141 (MQ=255)

    aTCTCCTTCATTTTTAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTAt                                                                                                                                      <  2:884042/141‑1 (MQ=255)

    aTCTCCTTCATTTTTAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTAt                                                                                                                                      >  4:246245/1‑141 (MQ=255)

     tctcCTTCATTTTTAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTAtt                                                                                                                                     <  1:706336/141‑1 (MQ=255)

      ctcCTTCATTTTTAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTAttt                                                                                                                                    <  1:38019/141‑1 (MQ=255)

       tcCTTCATTTTTAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTc                                                                                                                                   <  3:345062/141‑1 (MQ=255)

          ttCATTTTTAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTcaac                                                                                                                                <  3:415794/141‑1 (MQ=255)

             aTTTTTAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCaacaac                                                                                                                             <  3:9814/141‑1 (MQ=255)

              tttttAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAacaaca                                                                                                                            >  1:754740/1‑141 (MQ=255)

              tttttAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATGCCACCTATTTCAacaaca                                                                                                                            <  2:957552/141‑1 (MQ=255)

               ttttAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAacac                                                                                                                           <  1:676263/141‑1 (MQ=255)

                  tAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGAc                                                                                                                        >  3:364310/1‑141 (MQ=255)

                  tAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGAc                                                                                                                        <  4:271609/141‑1 (MQ=255)

                   agtttCGTCAATTTTATTTAGCATTGCTTGTGAGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTTTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACt                                                                                                                       >  2:860198/4‑141 (MQ=255)

                   aGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACt                                                                                                                       >  1:290590/1‑141 (MQ=255)

                   aGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACt                                                                                                                       >  1:290630/1‑141 (MQ=255)

                   aGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCACAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACt                                                                                                                       >  4:424163/1‑141 (MQ=255)

                    gCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTg                                                                                                                      >  1:216751/1‑141 (MQ=255)

                    gCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTg                                                                                                                      <  1:874134/141‑1 (MQ=255)

                       tCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGtat                                                                                                                   >  1:74569/1‑141 (MQ=255)

                       tCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATGCCACCTATTTCAACAACACGACTGtat                                                                                                                   >  4:323108/1‑141 (MQ=255)

                            aaTTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTGCACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATGCCACCTATTTCAACAACACGACTGTATATATc                                                                                                              <  4:433459/141‑1 (MQ=255)

                            aaTTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATc                                                                                                              <  3:375006/141‑1 (MQ=255)

                               tttATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGat                                                                                                           <  3:15432/141‑1 (MQ=255)

                                 tATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGatat                                                                                                         <  1:597867/141‑1 (MQ=255)

                                  aTCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATc                                                                                                        >  2:1000720/1‑141 (MQ=255)

                                  aTCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATc                                                                                                        <  3:303591/141‑1 (MQ=255)

                                  aTCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTAAAAAATAATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATc                                                                                                        <  2:1071221/141‑1 (MQ=255)

                                   tctcGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCt                                                                                                       >  3:70375/1‑141 (MQ=255)

                                    ctcGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTa                                                                                                      <  2:627113/141‑1 (MQ=255)

                                    ctcGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTa                                                                                                      <  1:159138/141‑1 (MQ=255)

                                     tcGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGCATCATAAAAATATTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTaa                                                                                                     >  1:810916/1‑141 (MQ=255)

                                       gCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACt                                                                                                   >  1:455389/1‑141 (MQ=255)

                                             cTTGGGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGta                                                                                             <  3:47496/141‑1 (MQ=255)

                                              ttGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGtaa                                                                                            >  4:287113/1‑141 (MQ=255)

                                              ttGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGtaa                                                                                            >  3:99963/1‑141 (MQ=255)

                                              ttGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGtaa                                                                                            <  4:350173/141‑1 (MQ=255)

                                               tgtgCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGtaat                                                                                           >  2:65812/1‑141 (MQ=255)

                                                 tgCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTCTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGaaataa                                                                                         >  2:475483/1‑141 (MQ=255)

                                                   cGATACTTGCATCAATTACTCGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAaca                                                                                       <  3:63276/141‑1 (MQ=255)

                                                   cGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCCTAAAAATCTTTATCTTTAATTCCACCTTTTTCAACAACACGACTGTATATATCGAGATCTAACTGTTGTAATAaca                                                                                       >  4:61346/1‑141 (MQ=255)

                                                    gATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAacac                                                                                      <  3:136132/141‑1 (MQ=255)

                                                      tACTGTCATCAATTACACGAACGTCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTTATTCCGCCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACtt                                                                                    <  2:99118/141‑1 (MQ=255)

                                                            catcaATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATa                                                                              <  2:400235/141‑1 (MQ=255)

                                                            catcaATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATa                                                                              >  3:300761/1‑141 (MQ=255)

                                                              tcaATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAg                                                                            <  4:399676/141‑1 (MQ=255)

                                                               caATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAgg                                                                           <  1:998446/141‑1 (MQ=255)

                                                                  ttACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCAc                                                                        >  1:937356/1‑141 (MQ=255)

                                                                   tACACGAACATCAGTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCAcc                                                                       <  4:189378/141‑1 (MQ=255)

                                                                    acacGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCg                                                                      >  2:543687/1‑141 (MQ=255)

                                                                    acacGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATGCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCg                                                                      <  3:323108/141‑1 (MQ=255)

                                                                      acGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAGCAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGAc                                                                    <  2:66710/141‑1 (MQ=255)

                                                                      acGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGAc                                                                    >  1:411249/1‑141 (MQ=255)

                                                                         aaCATNATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATCTAGATAGCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGc                                                                 >  1:826748/1‑141 (MQ=255)

                                                                         aaCATCACTACGACCAAGATCTGCTTTACATGTTTCTGAATCATAAAAATCTTTATCTTTCATTCCACCTATTTCAACCCCACGCCTGTATATATCGATATCTAACTGTTGTAATAACACTTTACCTAAGGCACCGACTGc                                                                 <  4:154732/141‑1 (MQ=255)

                                                                          aCATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCa                                                                >  4:171056/1‑141 (MQ=255)

                                                                            atcatTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAAc                                                              <  1:1041376/141‑1 (MQ=255)

                                                                            atcatTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATGCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAAc                                                              <  2:850652/141‑1 (MQ=255)

                                                                              catTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTc                                                            <  2:327747/141‑1 (MQ=255)

                                                                               atTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCg                                                           <  2:411249/141‑1 (MQ=255)

                                                                               atTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCg                                                           <  1:471278/141‑1 (MQ=255)

                                                                               atTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGAAACTCg                                                           >  2:796120/1‑141 (MQ=255)

                                                                                 tACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAg                                                         >  2:707373/1‑141 (MQ=255)

                                                                                    gATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAgctg                                                      <  1:816461/141‑1 (MQ=255)

                                                                                      tCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCACCAACACGCCTGCATATCTCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAgctgct                                                    <  3:456449/141‑1 (MQ=255)

                                                                                      tCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACCCGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAgctgct                                                    <  3:456434/141‑1 (MQ=255)

                                                                                        aaGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGt                                                  <  4:85892/141‑1 (MQ=255)

                                                                                        aaGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATGCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGt                                                  >  2:103131/1‑141 (MQ=255)

                                                                                          gATTTGCTTTAAATGTTTCGGAATCATAAAAATCTTTATCTTTAATTCCACCTATCTCAACAACACGACTGTATATCTCGATCTCTAACTGCTGTAACAACACTTTACACACCGCACCGACTGCAACCCGAGctgctgccc                                                >  2:626439/1‑138 (MQ=255)

                                                                                           aTTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTtc                                               <  2:155904/141‑1 (MQ=255)

                                                                                               gCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTa                                           >  4:210666/1‑141 (MQ=255)

                                                                                               gCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTa                                           <  3:444187/141‑1 (MQ=255)

                                                                                                 tttAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGc                                         >  2:968519/1‑141 (MQ=255)

                                                                                                  ttAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCa                                        >  1:759561/1‑141 (MQ=255)

                                                                                                   tAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATACATCCATATCTAACTGTTGTAATAACACTTTACATAAGGCACCCACTGCAACCCGAGCTGCTGTTTCTCTCGCAg                                       >  2:515196/1‑141 (MQ=255)

                                                                                                      aTGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCTCCTATTCCAACAACACTACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTTGACCCGTTGTTTCTCTAGCAgatg                                    <  4:457556/141‑1 (MQ=255)

                                                                                                      aTGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAgatg                                    >  1:605785/1‑141 (MQ=255)

                                                                                                          ttCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCg                                <  4:187362/141‑1 (MQ=255)

                                                                                                          ttCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCg                                <  2:361920/141‑1 (MQ=255)

                                                                                                           tCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGc                               >  1:378480/1‑141 (MQ=255)

                                                                                                           tCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGc                               <  3:243566/141‑1 (MQ=255)

                                                                                                             tGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGctc                             >  3:115386/1‑141 (MQ=255)

                                                                                                                 tCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGc                         <  2:533848/141‑1 (MQ=255)

                                                                                                                  cATAAAAATCTTTATCTTTAATTACACCTATTTCAACAACAAGACTGTCTATATCAATATCTAAATGTTGTACTAACACTTACCAAAAGGGAAAGAGTGCAACACGCGCTGCGGTTTCTCTCGCAGATGCTCGCTCTAGcc                        >  1:158090/1‑140 (MQ=255)

                                                                                                                    tAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGcaca                      >  1:258901/1‑141 (MQ=255)

                                                                                                                      aaaaTCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGTACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGCACAtt                    >  1:68605/1‑141 (MQ=255)

                                                                                                                       aaaTCTTTATCTTTAATGCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGCACAttt                   >  2:892678/1‑141 (MQ=255)

                                                                                                                        aaTCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTAAATAGGGTACTGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGCACATTTc                  >  3:327158/1‑141 (MQ=255)

                                                                                                                           ctttatctttaATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGCACATTTCGTa               <  1:197786/141‑1 (MQ=255)

                                                                                                                           ctttatctttaATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGCACATTTCGTa               <  3:171754/141‑1 (MQ=255)

                                                                                                                             ttatctttaATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTATCTAGCAGATGATCGCTCTAGCACATTTCGTaaa             >  1:1006520/1‑141 (MQ=255)

                                                                                                                             ttatctttaATGCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGCACATTTCGTaaa             >  1:518546/1‑141 (MQ=255)

                                                                                                                                  tttaATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGCACATTTCGTAAATCACg        <  1:879342/141‑1 (MQ=255)

                                                                                                                                    taATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGCACATTTCGTAAATCACgat      <  4:8152/141‑1 (MQ=255)

                                                                                                                                       ttCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGCACATTTCGTAAATCACgatgat   <  1:788103/141‑1 (MQ=255)

                                                                                                                                       ttCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGCACATTTCGTAAATCACgatgat   <  2:290590/141‑1 (MQ=255)

                                                                                                                                       ttCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGCACATTTCGTAAATCACgatgat   <  2:290630/141‑1 (MQ=255)

                                                                                                                                        tCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTATCACATTTCGTAAATCACGATGAtt  >  1:1024808/1‑141 (MQ=255)

TTGAATCTCCTTCATTTTTAGCTTCGTCAATTTTATCTCGCATTGCTTGTGCGATACTGTCATCAATTACACGAACATCATTACGATCAAGATTTGCTTTAAATGTTTCTGAATCATAAAAATCTTTATCTTTAATTCCACCTATTTCAACAACACGACTGTATATATCGATATCTAACTGTTGTAATAACACTTTACATAAGGCACCGACTGCAACTCGAGCTGCTGTTTCTCTAGCAGATGATCGCTCTAGCACATTTCGTAAATCACGATGATT  >  USA300TCH1516_ALE/1538262‑1538538

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