BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,854,119	C→T	R587C (CGC→TGC)	fhlA →	FhlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,854,119	0	C	T	84.0%	767.9 / 98.3	293	R587C (CGC→TGC)	fhlA	FhlA
Reads supporting (aligned to +/- strand): ref base C (24/23); new base T (141/105); total (165/128)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 296 total aligned reads displayed.

CGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGC  >  NC_000913/2854027‑2854209

   aaGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCa                                                                                         >  2:5326019/1‑93 (MQ=255)

     gATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCAcc                                                                                       >  1:4994883/1‑93 (MQ=255)

       tatTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCCACAGCATTCCTGCCGAGACGCTGTGCACCtt                                                                                     >  1:5363968/1‑93 (MQ=255)

        atTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTg                                                                                    >  1:3808840/1‑93 (MQ=255)

         ttCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGa                                                                                   <  2:7065047/93‑1 (MQ=255)

           ccGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGc                                                                                 >  2:3875834/1‑93 (MQ=255)

            cGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCa                                                                                >  2:5432094/1‑93 (MQ=255)

             gctgctGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCaa                                                                               >  1:7214113/1‑93 (MQ=255)

              ctgctgGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAAc                                                                              >  2:4057471/1‑93 (MQ=255)

              ctgctgGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAAc                                                                              <  1:6372719/93‑1 (MQ=255)

               tgctgGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACa                                                                             >  2:7307780/1‑93 (MQ=255)

               tgctgGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACa                                                                             >  1:6929912/1‑93 (MQ=255)

               tgctgGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACa                                                                             >  2:3075770/1‑93 (MQ=255)

                gctgGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACAt                                                                            <  2:112682/93‑1 (MQ=255)

                gctgGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACAt                                                                            <  2:6707262/93‑1 (MQ=255)

                gctgGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACAt                                                                            >  2:1361895/1‑93 (MQ=255)

                 ctgGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATg                                                                           <  1:4685426/93‑1 (MQ=255)

                  tgGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACAtnn                                                                          >  1:2371989/1‑91 (MQ=255)

                  tgGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATgg                                                                          >  1:3481523/1‑93 (MQ=255)

                   ggCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGa                                                                         <  1:3207394/93‑1 (MQ=255)

                    gCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAg                                                                        >  1:6028341/1‑93 (MQ=255)

                     cGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGt                                                                       <  2:452150/93‑1 (MQ=255)

                      gAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTg                                                                      >  1:1034294/1‑93 (MQ=255)

                      gAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTg                                                                      >  1:2558144/1‑93 (MQ=255)

                       aaaGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTgg                                                                     >  1:5831124/1‑93 (MQ=255)

                        aaGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACGTGGAGTGGc                                                                    <  1:6751671/93‑1 (MQ=255)

                         aGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGcc                                                                   <  2:6305425/93‑1 (MQ=255)

                         aGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGcc                                                                   >  1:5542829/1‑93 (MQ=255)

                          gCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCg                                                                  >  2:724091/1‑93 (MQ=255)

                            cTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCggg                                                                >  1:6804253/1‑93 (MQ=255)

                              ttACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTa                                                              <  2:1861648/93‑1 (MQ=255)

                                aCCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAAc                                                            >  1:7090724/1‑93 (MQ=255)

                                aCCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGCGGCCGGGGAAc                                                            >  2:760671/1‑93 (MQ=255)

                                 ccTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACg                                                           <  1:5703776/93‑1 (MQ=255)

                                 ccTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACg                                                           >  1:4514601/1‑93 (MQ=255)

                                  cTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGt                                                          <  2:871200/93‑1 (MQ=255)

                                  cTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGt                                                          >  2:2160775/1‑93 (MQ=255)

                                  cTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGt                                                          <  2:4270928/93‑1 (MQ=255)

                                   ttCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTa                                                         <  2:6912191/93‑1 (MQ=255)

                                   ttCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGGAACGTa                                                         >  1:4938437/1‑93 (MQ=255)

                                    tCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTAc                                                        <  1:1842362/93‑1 (MQ=255)

                                     cAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTAcg                                                       >  2:198/1‑93 (MQ=255)

                                      aaaaTTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTAcgc                                                      >  1:2070907/1‑93 (MQ=255)

                                      aaaaTTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTAcgc                                                      >  1:7305912/1‑93 (MQ=255)

                                      aaaaTTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTAcgc                                                      <  2:1619139/93‑1 (MQ=255)

                                         aTTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGaa                                                   >  2:5604485/1‑93 (MQ=255)

                                         aTTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGaa                                                   >  2:1916945/1‑93 (MQ=255)

                                         aTTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGaa                                                   <  2:1224013/93‑1 (MQ=255)

                                          ttGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAAc                                                  >  2:3303960/1‑93 (MQ=255)

                                           tGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACt                                                 <  2:1693571/93‑1 (MQ=255)

                                           tGCCCGCCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCAGGGTAACGTACGCGAACt                                                 <  1:1321127/93‑1 (MQ=255)

                                            gCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTg                                                <  2:2587435/93‑1 (MQ=255)

                                            gCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTg                                                <  2:338689/93‑1 (MQ=255)

                                              ccGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGa                                              >  1:6180925/1‑93 (MQ=255)

                                              ccGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGa                                              >  1:3589419/1‑93 (MQ=255)

                                               cgtcgtCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGaa                                             >  2:4422044/1‑93 (MQ=255)

                                                gtcgtcTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGaaa                                            >  1:940218/1‑93 (MQ=255)

                                                 tcgtcTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGaaaa                                           <  1:1721171/93‑1 (MQ=255)

                                                   gtcTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACg                                         >  1:6545623/1‑93 (MQ=255)

                                                    tcTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGt                                        >  2:3643187/1‑93 (MQ=255)

                                                       ggggCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCAt                                     >  1:2956152/1‑93 (MQ=255)

                                                       ggggCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCAt                                     >  2:2972194/1‑93 (MQ=255)

                                                         ggCGCAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTg                                   >  2:4851730/1‑93 (MQ=255)

                                                         ggCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTg                                   <  1:3857829/93‑1 (MQ=255)

                                                          gggcAATATCGACAGCATGCCTGCCGAGACGCTGTGCACCTTGAGCAACATGAAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGa                                  <  1:6578978/91‑1 (MQ=255)

                                                           cgcAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAg                                 <  2:4218921/93‑1 (MQ=255)

                                                            gcAATATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAgc                                >  1:5906079/1‑93 (MQ=255)

                                                              aaTATCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAgcgc                              >  1:1557852/1‑93 (MQ=255)

                                                              aaTATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAgcgc                              >  2:5741127/1‑93 (MQ=255)

                                                               atatCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAgcgcg                             <  1:3979497/93‑1 (MQ=255)

                                                               atatCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAgcgcg                             >  1:6335240/1‑93 (MQ=255)

                                                                tatCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAgcgcgc                            >  2:3859035/1‑93 (MQ=255)

                                                                 atCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGcgcgcg                           >  1:5751551/1‑93 (MQ=255)

                                                                  tCGACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCgg                          <  2:4134348/93‑1 (MQ=255)

                                                                    gACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTa                        >  1:5231045/1‑93 (MQ=255)

                                                                    gACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTa                        >  1:3507724/1‑93 (MQ=255)

                                                                    gACAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTa                        <  2:1076678/93‑1 (MQ=255)

                                                                     aCAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTAt                       <  2:3503019/93‑1 (MQ=255)

                                                                     aCAGCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTAt                       >  1:5190976/1‑93 (MQ=255)

                                                                      cAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTAtt                      <  2:2352134/93‑1 (MQ=255)

                                                                      cAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTAtt                      >  2:3750190/1‑93 (MQ=255)

                                                                        gCATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGc                    >  2:5532276/1‑93 (MQ=255)

                                                                         cATTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCt                   <  2:5237676/93‑1 (MQ=255)

                                                                          aTTCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTa                  >  2:2673083/1‑93 (MQ=255)

                                                                           ttCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTaa                 >  1:1441409/1‑93 (MQ=255)

                                                                           ttCCTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATGGCTaa                 <  2:4592419/93‑1 (MQ=255)

                                                                              cTGCCGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAacac              >  2:3012710/1‑93 (MQ=255)

                                                                                  cGAGACGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCgg          >  2:6631694/1‑93 (MQ=255)

                                                                                    agaCGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTa        <  1:7359484/93‑1 (MQ=255)

                                                                                     gaCGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTaa       <  2:4745081/93‑1 (MQ=255)

                                                                                     gaCGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTaa       <  2:4469871/93‑1 (MQ=255)

                                                                                     gaCGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTaa       >  1:2843212/1‑93 (MQ=255)

                                                                                     gaCGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTaa       <  1:4227159/93‑1 (MQ=255)

                                                                                     gaCGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTaa       <  2:6066438/93‑1 (MQ=255)

                                                                                      aCGCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAAc      >  2:95233/1‑93 (MQ=255)

                                                                                        gCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGt    >  2:2817298/1‑93 (MQ=255)

                                                                                        gCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGt    >  2:2095463/1‑93 (MQ=255)

                                                                                        gCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGt    <  1:1062130/93‑1 (MQ=255)

                                                                                        gCTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGAGGTAACGt    <  1:5120961/93‑1 (MQ=255)

                                                                                         cTGTGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGtg   <  1:788184/93‑1 (MQ=255)

CGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGC  >  NC_000913/2854027‑2854209

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