BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BX571857	2,489,372	G→C	6.9%	A172P (GCC→CCC)	SAS2328 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BX571857	2,489,372	0	G	C	6.9%	892.7 / 45.5	218	A172P (GCC→CCC)	SAS2328	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (91/112); new base C (7/8); total (98/120)
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.16e-01

Only 100 of 223 total aligned reads displayed.

ATATTCACGTAACGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGGAAGAAAGTGCATTTATGTGCAAC  >  BX571857/2489277‑2489468

atatTCACGTAACGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGc                                                                                                 >  2:612412/1‑97 (MQ=255)

atatTCACGTAACGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGc                                                                                                 <  2:2678562/97‑1 (MQ=255)

  atTCACGTAACGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCca                                                                                               >  1:2455512/1‑97 (MQ=255)

      aCGTAACGTTTCTCATTGGGGACGATGATTACTATTACTATTATATCGGTCTACGTTGGAGTCATATTCGTCGTGGCCTAGTTTTACGTGCcatcat                                                                                           <  2:2338549/97‑1 (MQ=255)

       cGTAACGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCAtt                                                                                          <  2:1339626/97‑1 (MQ=255)

          aaCGTTACTCGTTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCa                                                                                       <  2:96692/97‑1 (MQ=255)

          aaCGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCa                                                                                       <  2:1047708/97‑1 (MQ=255)

          aaCGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTCCCATCATTCCa                                                                                       <  2:2726759/97‑1 (MQ=255)

           aCGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCCCTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTCCCATCATTCCaa                                                                                      >  2:482484/1‑97 (MQ=255)

           aCGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCaa                                                                                      <  2:144260/97‑1 (MQ=255)

           aCGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCaa                                                                                      >  1:406309/1‑97 (MQ=255)

           aCGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCaa                                                                                      <  1:2688161/97‑1 (MQ=255)

            cGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCaaa                                                                                     >  2:2788830/1‑97 (MQ=255)

            cGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCaaa                                                                                     >  2:1904066/1‑97 (MQ=255)

              ttACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATc                                                                                   >  1:441595/1‑97 (MQ=255)

                aCTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGc                                                                                 >  2:2585983/1‑97 (MQ=255)

                aCTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGc                                                                                 >  1:956687/1‑97 (MQ=255)

                  tCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGcaa                                                                               <  1:1414973/97‑1 (MQ=255)

                  tCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTCCCATCATTCCAAATCGcaa                                                                               <  2:1611424/97‑1 (MQ=255)

                   cATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGcaaa                                                                              >  1:1370771/1‑97 (MQ=255)

                     ttGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGcaaaca                                                                            <  2:66525/97‑1 (MQ=255)

                     ttGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGcaaaca                                                                            <  1:391294/97‑1 (MQ=255)

                     ttGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGcaaaca                                                                            <  1:256613/97‑1 (MQ=255)

                      tGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGcaaacaa                                                                           <  1:608523/97‑1 (MQ=255)

                        ggTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAat                                                                         <  1:2770977/97‑1 (MQ=255)

                         gTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAatt                                                                        <  1:2635636/97‑1 (MQ=255)

                           aCGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAattat                                                                      <  1:209684/97‑1 (MQ=255)

                           aCGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAattat                                                                      >  1:2964029/1‑97 (MQ=255)

                            cGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAattatt                                                                     >  1:1083543/1‑97 (MQ=255)

                                gatTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTaaaa                                                                 <  1:2130754/97‑1 (MQ=255)

                                gatTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTaaaa                                                                 >  2:118642/1‑97 (MQ=255)

                                 attactattactATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAAt                                                                <  2:12482/97‑1 (MQ=255)

                                  ttactattactaTTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAAtt                                                               <  2:1114225/97‑1 (MQ=255)

                                   tactattactatTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTc                                                              >  1:285949/1‑97 (MQ=255)

                                   tactattactatTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTc                                                              >  2:2482688/1‑97 (MQ=255)

                                   tactattactatTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTCCCATCATTCCAAATCGCAAACAAATTATTAAAATTc                                                              >  1:2195776/1‑97 (MQ=255)

                                   tactattactatTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTCTACCTCCCATCATTCCAAATCGCAAACAAATTATTAAAAtcc                                                              >  2:1555780/1‑95 (MQ=255)

                                    actattactattATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCg                                                             >  1:738974/1‑97 (MQ=255)

                                     ctattactattaTATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGt                                                            <  1:1958650/97‑1 (MQ=255)

                                      tattactattaTATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTa                                                           <  2:1957287/97‑1 (MQ=255)

                                         tactattaTATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAag                                                        <  2:941049/97‑1 (MQ=255)

                                           ctattaTATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAagag                                                      >  1:1434109/1‑97 (MQ=255)

                                           ctattaTATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAagag                                                      >  2:2532508/1‑97 (MQ=255)

                                            tattatATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAgaga                                                     >  1:194680/1‑97 (MQ=255)

                                            tattatATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAgaga                                                     >  1:2546693/1‑97 (MQ=255)

                                            tattatATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTCCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAgaga                                                     >  1:2472378/1‑97 (MQ=255)

                                            tattatATCGATCCACTTTGGAGTCATATTCGTCGTGGACTAGTTTTACCAGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAgaga                                                     <  1:1585651/97‑1 (MQ=255)

                                             attatATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGaa                                                    >  2:1508712/1‑97 (MQ=255)

                                             attatATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGaa                                                    <  2:2610025/97‑1 (MQ=255)

                                              ttatATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGaaa                                                   <  2:1997385/97‑1 (MQ=255)

                                              ttatATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGaaa                                                   >  1:1868921/1‑97 (MQ=255)

                                               tataTCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGaaaa                                                  <  1:828052/97‑1 (MQ=255)

                                                  atCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATgg                                               <  2:942686/97‑1 (MQ=255)

                                                   tCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGa                                              <  1:1751138/97‑1 (MQ=255)

                                                   tCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGa                                              >  2:180131/1‑97 (MQ=255)

                                                   tCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGa                                              >  2:1426566/1‑97 (MQ=255)

                                                   tCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGa                                              >  2:1174582/1‑97 (MQ=255)

                                                   tCGATCCACTTTGGAGTCATATTCGTCGTGGACTAGGGTTGCCTGCCATCATTCCAAGTCGCAGACAAATTATTAAAATTCGGAAAGATAAAAGGGa                                              >  2:68166/1‑97 (MQ=255)

                                                    cGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGaa                                             <  1:1125814/97‑1 (MQ=255)

                                                      aTCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATa                                           <  2:2242918/97‑1 (MQ=255)

                                                       tCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATAc                                          >  1:148211/1‑97 (MQ=255)

                                                        ccACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTCCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACa                                         >  2:2880572/1‑97 (MQ=255)

                                                         cACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAg                                        <  1:1216600/97‑1 (MQ=255)

                                                         cACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAg                                        <  1:2368927/97‑1 (MQ=255)

                                                         cACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAg                                        >  1:912616/1‑97 (MQ=255)

                                                         cACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAg                                        >  1:2716900/1‑97 (MQ=255)

                                                          aCTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTCCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGc                                       <  1:178535/97‑1 (MQ=255)

                                                             ttGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCttt                                    <  2:138403/97‑1 (MQ=255)

                                                              tGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTa                                   >  2:2378600/1‑97 (MQ=255)

                                                              tGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTa                                   <  2:1575504/97‑1 (MQ=255)

                                                               ggAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTaa                                  <  1:41882/97‑1 (MQ=255)

                                                               ggAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTaa                                  <  2:450293/97‑1 (MQ=255)

                                                               ggAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGGAAAGAGAAAATGGAATACAGCTTTaa                                  <  2:79559/97‑1 (MQ=255)

                                                                gAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAAt                                 <  2:2746094/97‑1 (MQ=255)

                                                                 aGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAAtt                                <  1:2378458/97‑1 (MQ=255)

                                                                 aGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAAtt                                >  1:2409759/1‑97 (MQ=255)

                                                                   tCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCa                              >  1:1042262/1‑97 (MQ=255)

                                                                    cATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCAc                             <  2:723059/97‑1 (MQ=255)

                                                                    cATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCAc                             >  2:2895157/1‑97 (MQ=255)

                                                                     atatTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCAcc                            >  1:2020833/1‑97 (MQ=255)

                                                                     atatTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCAcc                            >  1:1642238/1‑97 (MQ=255)

                                                                     atatTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCAcc                            >  1:3001856/1‑97 (MQ=255)

                                                                        ttCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAgg                         <  2:2473247/97‑1 (MQ=255)

                                                                        ttCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAgg                         >  2:1074041/1‑97 (MQ=255)

                                                                           gtcgtGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGgaag                      <  2:2959902/97‑1 (MQ=255)

                                                                           gtcgtGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGgaag                      >  1:1613131/1‑97 (MQ=255)

                                                                             cgtGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGgaagaa                    <  1:1274614/97‑1 (MQ=255)

                                                                              gtGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGGAAGaaa                   <  1:544839/97‑1 (MQ=255)

                                                                              gtGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGGAAGaaa                   >  1:27679/1‑97 (MQ=255)

                                                                              gtGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGGAAGaaa                   <  1:2587640/97‑1 (MQ=255)

                                                                               tGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGGAAGAAAg                  <  2:2166961/97‑1 (MQ=255)

                                                                                     aGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGGAAGAAAGTGCAtt            <  1:589841/97‑1 (MQ=255)

                                                                                     aGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGGAAGAAAGTGCAtt            >  1:2919241/1‑97 (MQ=255)

                                                                                     aGTTTTACCTGCCATCATTCCAAAGCGCAAACAAATTATTAAAATGCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGGAAGAAAGTGCAtt            <  2:1199558/97‑1 (MQ=255)

                                                                                      gTTTTACCTGCCATCATTACAAATCGCAAACAGATTATTACATTTAGTAAAGAGAAAATGGACTACAGCTTTAATTCACCAGGACGAACGCGCAttt           >  2:892028/1‑97 (MQ=255)

                                                                                        tttACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGGAAGAAAGTGCATTTAt         <  1:2183134/97‑1 (MQ=255)

                                                                                         ttACCTGCCAGCATTCCAAATCGCAGGCAAATTATTGAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGGAAGAAAGTGCATTTAtg        <  2:394347/97‑1 (MQ=255)

                                                                                           aCCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGGAAGAAAGTGCATTTAtgtg      <  1:362606/97‑1 (MQ=255)

                                                                                              tGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATGCACCAGGAAGAAAGTGCATTTAAGTGCaa   >  2:575328/1‑97 (MQ=255)

                                                                                               gCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGGAAGAAAGTGCATTTATGTGCAAc  >  2:784277/1‑97 (MQ=255)

ATATTCACGTAACGTTACTCATTGGGTACGATGATTACTATTACTATTATATCGATCCACTTTGGAGTCATATTCGTCGTGGCCTAGTTTTACCTGCCATCATTCCAAATCGCAAACAAATTATTAAAATTCGTAAAGAGAAAATGGAATACAGCTTTAATTCACCAGGAAGAAAGTGCATTTATGTGCAAC  >  BX571857/2489277‑2489468

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