BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	740,328:1	+G	5.6%	intergenic (+555/+33)	pitA_1 → / ← sle1_2	Low‑affinity inorganic phosphate transporter 1/N‑acetylmuramoyl‑L‑alanine amidase sle1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	740,328	1	.	G	5.6%	2008.8 / 117.6	429	intergenic (+555/+33)	pitA_1/sle1_2	Low‑affinity inorganic phosphate transporter 1/N‑acetylmuramoyl‑L‑alanine amidase sle1
Reads supporting (aligned to +/- strand): ref base . (154/251); new base G (15/9); total (169/260)
Fisher's exact test for biased strand distribution p-value = 2.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 450 total aligned reads displayed.

AATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCATA  >  USA300TCH1516_ALE/740237‑740418

aatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAAT‑TCACATCATGTGa                                                                                             <  1:3988323/91‑1 (MQ=255)

 atatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGa                                                                                             <  1:512750/91‑1 (MQ=255)

    tgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑Att                                                                                         >  1:4770653/1‑91 (MQ=255)

    tgtatgATAGAGAAGGGTGAGATAAATTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑Att                                                                                         <  2:4662699/91‑1 (MQ=255)

     gtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTg                                                                                        <  1:2603597/91‑1 (MQ=255)

       atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCt                                                                                      <  2:2551039/91‑1 (MQ=255)

        tggtAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTa                                                                                     >  1:549934/4‑91 (MQ=255)

        tgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAA‑TTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAg                                                                                    <  2:3462029/91‑1 (MQ=255)

         gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAg                                                                                    <  1:6589288/91‑1 (MQ=255)

         gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTAAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAg                                                                                    <  2:918382/91‑1 (MQ=255)

          aTAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGt                                                                                   <  1:3736117/91‑1 (MQ=255)

           tAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGta                                                                                  <  2:5983837/91‑1 (MQ=255)

              ggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAAT‑TCACATCATGTGAGATTGCTAGttttt                                                                               <  2:3695957/91‑5 (MQ=25)

               gAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGtatata                                                                              <  1:4906615/91‑1 (MQ=255)

                aaGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTatatat                                                                             <  1:6285926/91‑1 (MQ=255)

                 aGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATc                                                                            >  1:6815879/1‑91 (MQ=255)

                 aGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATc                                                                            >  1:2002628/1‑91 (MQ=255)

                 aGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATc                                                                            <  2:3767409/91‑1 (MQ=255)

                 aGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATc                                                                            >  2:315854/1‑91 (MQ=255)

                  gggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCa                                                                           >  1:4263354/1‑91 (MQ=255)

                  gggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTAAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCa                                                                           <  2:1818149/91‑1 (MQ=255)

                   ggTGAGATAATTTTTTGTACTAAAATGAAAATTTTTTTCCCAAATAAAAAACTCGCAATCTCACATGATGTGC‑CTTGCTCGCACATATCAg                                                                          >  2:7006059/1‑91 (MQ=255)

                   ggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAg                                                                          >  2:1575793/1‑91 (MQ=255)

                   ggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAg                                                                          <  1:4961774/91‑1 (MQ=255)

                    gTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGt                                                                         <  2:2353153/91‑1 (MQ=255)

                     tGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTa                                                                        >  2:5740314/1‑91 (MQ=255)

                      gagaTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTAc                                                                       <  1:583742/91‑1 (MQ=255)

                      gagaTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTAc                                                                       >  2:6790642/1‑91 (MQ=255)

                      gagaTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTAc                                                                       <  1:123460/91‑1 (MQ=255)

                       agaTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACa                                                                      >  2:2154294/1‑91 (MQ=255)

                        gaTAATTTTTTGTACTTAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACaa                                                                     >  2:5116271/1‑91 (MQ=255)

                           aaTTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAattt                                                                  <  2:622402/91‑1 (MQ=255)

                             ttttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAatttat                                                                <  1:5211521/91‑1 (MQ=255)

                               ttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAatttattt                                                              <  2:5214638/91‑1 (MQ=255)

                               ttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAatttattt                                                              <  1:2940391/91‑1 (MQ=255)

                                 ttGTACCAATATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaa                                                            <  1:2556181/91‑1 (MQ=255)

                                 ttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaa                                                            <  2:4426687/91‑1 (MQ=255)

                                 ttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaa                                                            <  2:6130249/91‑1 (MQ=255)

                                 ttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaa                                                            <  2:5893798/91‑1 (MQ=255)

                                 ttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaa                                                            >  2:12501/1‑91 (MQ=255)

                                  tGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaat                                                           >  1:5045018/1‑91 (MQ=255)

                                  tGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaat                                                           >  1:5283927/1‑91 (MQ=255)

                                   gTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaatt                                                          >  1:5763818/1‑91 (MQ=255)

                                    tACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaatta                                                         >  1:2294842/1‑91 (MQ=255)

                                    tACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaatta                                                         >  1:3170780/1‑91 (MQ=255)

                                     aCCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaattaa                                                        <  1:5580379/91‑1 (MQ=255)

                                     aCCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaattaa                                                        >  1:2373613/1‑91 (MQ=255)

                                     aCCAAAATGAAAATATTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaattaa                                                        <  2:3570024/91‑1 (MQ=255)

                                     aCCAAAATGAAAATAATATTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATttaattaa                                                        <  2:1325773/91‑1 (MQ=255)

                                      ccAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTtaattaat                                                       >  1:2757754/1‑91 (MQ=255)

                                        aaaaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATgg                                                     <  1:321647/91‑1 (MQ=255)

                                         aaaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGa                                                    <  1:3542843/91‑1 (MQ=255)

                                           aTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATg                                                  <  1:3328520/91‑1 (MQ=255)

                                             gAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGaa                                                <  1:918612/91‑1 (MQ=255)

                                             gAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTACATATCAGTACAATTTATTTAATTAATGGATGaa                                                <  1:6946961/91‑1 (MQ=255)

                                              aaaaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAAt                                               >  2:6433201/1‑91 (MQ=255)

                                               aaaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATg                                              >  1:6920199/1‑91 (MQ=255)

                                               aaaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATg                                              <  2:388025/91‑1 (MQ=255)

                                                aaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAattt                                                                  >  1:5886777/1‑70 (MQ=255)

                                                aaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGc                                             >  2:1028604/1‑91 (MQ=255)

                                                aaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGc                                             >  1:2171807/1‑91 (MQ=255)

                                                 aTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCa                                            >  2:6518802/1‑91 (MQ=255)

                                                 aTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCa                                            <  1:462425/91‑1 (MQ=255)

                                                  tttttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCAt                                           <  1:345970/91‑1 (MQ=255)

                                                     ttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGc                                        <  2:4741761/91‑1 (MQ=255)

                                                      tttCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCt                                       >  1:6013621/1‑91 (MQ=255)

                                                       ttCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTa                                      >  2:5011854/1‑91 (MQ=255)

                                                        tCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAg                                     <  1:4121176/91‑1 (MQ=255)

                                                          aaaaaTAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaa                                   >  2:2589056/1‑91 (MQ=255)

                                                          aaaaaTAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaa                                   <  2:2757611/91‑1 (MQ=255)

                                                          aaaaaTAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaa                                   <  2:2203585/91‑1 (MQ=255)

                                                           aaaaTAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaaa                                  >  2:2863926/1‑91 (MQ=255)

                                                           aaaaTAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaaa                                  <  1:6479743/91‑1 (MQ=255)

                                                           aaaaTAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaaa                                  <  1:444939/91‑1 (MQ=255)

                                                              aTAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACtt                               <  1:6117228/91‑1 (MQ=255)

                                                              aTAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACtt                               <  1:2221774/91‑1 (MQ=255)

                                                               tAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTc                              <  2:2662476/91‑1 (MQ=255)

                                                                aaaaaaCTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCt                             <  2:3301982/91‑1 (MQ=255)

                                                                   aaaCTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGaa                          <  1:5375953/91‑1 (MQ=255)

                                                                   aaaCTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGaa                          >  2:4246786/1‑91 (MQ=255)

                                                                    aaCTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg                         >  1:1350338/1‑91 (MQ=255)

                                                                     aCTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGc                        <  1:2427678/91‑1 (MQ=255)

                                                                      cTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCt                       <  1:4317354/91‑1 (MQ=255)

                                                                         gCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCt                                       >  2:6653717/1‑72 (MQ=255)

                                                                         gCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGa                    <  2:5052281/91‑1 (MQ=255)

                                                                         gCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGa                    <  1:3644144/91‑1 (MQ=255)

                                                                         gCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGa                    >  1:1626243/1‑91 (MQ=255)

                                                                           aaTCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAAt                  >  1:3847019/1‑91 (MQ=255)

                                                                             attcACATCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAAtt                 >  1:1017329/3‑91 (MQ=39)

                                                                              ctcACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGt               <  1:1108580/91‑1 (MQ=255)

                                                                               tcACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa              <  2:4894349/91‑1 (MQ=255)

                                                                               tcACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa              <  1:1925273/91‑1 (MQ=255)

                                                                                cacaTGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTAc             >  1:3356669/1‑91 (MQ=255)

                                                                                cacaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa              >  1:6224931/1‑91 (MQ=255)

                                                                                 acaTGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACg            <  1:3673119/91‑1 (MQ=255)

                                                                                   atgatgTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGt          <  2:224610/91‑1 (MQ=255)

                                                                                     gatgTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACgggtg        >  2:5798820/1‑88 (MQ=255)

                                                                                     gatgTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAg        <  2:4781041/91‑1 (MQ=255)

                                                                                        gtgA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTc     <  1:724665/91‑1 (MQ=255)

                                                                                           a‑aTTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCata  >  2:4031220/1‑91 (MQ=255)

AATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGA‑ATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCATA  >  USA300TCH1516_ALE/740237‑740418

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