BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,484,849	G→A	R355H (CGC→CAC)	yjgP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,484,849	0	G	A	96.9%	531.1 / ‑4.2	162	R355H (CGC→CAC)	yjgP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/4); new base A (64/92); total (66/96)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

GACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGTCCGCCGCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTCACCA  >  NC_000913/4484761‑4484937

gACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTg                                                                                      >  1:2715491/1‑93 (MQ=255)

gACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTg                                                                                      <  1:1770973/93‑1 (MQ=255)

   ccGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTgcgc                                                                                   <  2:3505417/93‑1 (MQ=255)

      aCGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGcc                                                                                <  1:1444022/93‑1 (MQ=255)

         cTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGt                                                                             >  2:2864021/1‑93 (MQ=255)

          tgtggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGtt                                                                            >  2:3565806/1‑93 (MQ=255)

          tgtggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGtt                                                                            >  1:4283089/1‑93 (MQ=255)

          tgtggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGtt                                                                            >  1:1730056/1‑93 (MQ=255)

                tgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTcgc                                                                      >  2:2525541/1‑93 (MQ=255)

                 gtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTcgcg                                                                     <  1:4467049/93‑1 (MQ=255)

                  tggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGt                                                                    >  1:1257686/1‑93 (MQ=255)

                   ggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTa                                                                   <  1:1949410/93‑1 (MQ=255)

                   ggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTa                                                                   >  1:3161689/1‑93 (MQ=255)

                   ggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTa                                                                   <  2:812721/93‑1 (MQ=255)

                     aCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTaaa                                                                 <  1:3585498/93‑1 (MQ=255)

                     aCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTaaa                                                                 >  1:4471609/1‑93 (MQ=255)

                      ccGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAg                                                                <  1:3750757/93‑1 (MQ=255)

                      ccGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAg                                                                <  1:297649/93‑1 (MQ=255)

                      ccGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAg                                                                >  1:709082/1‑93 (MQ=255)

                      ccGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCACGCCACCTGCGCGCCAGTTTTTCGCGTAAAg                                                                >  1:7157/1‑93 (MQ=255)

                       cGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCAGTCCGCCGCCTGCGCGCCAGTTTTTCGCGTAAAgg                                                               <  1:4280475/93‑1 (MQ=255)

                        gTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGTCCGCCGCCTGCGCGCCAGTTTTTCGCGTAAAGGa                                                              >  2:3824999/1‑93 (MQ=255)

                        gTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGa                                                              >  2:2538973/1‑93 (MQ=255)

                        gTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGa                                                              >  2:436589/1‑93 (MQ=255)

                         ttAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCTGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAg                                                             <  1:41060/93‑1 (MQ=255)

                         ttAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAg                                                             >  1:1796820/1‑93 (MQ=255)

                          tAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGc                                                            >  2:81863/1‑93 (MQ=255)

                           aaCCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCg                                                           >  1:3903662/1‑93 (MQ=255)

                             ccTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGgt                                                         <  2:2764715/93‑1 (MQ=255)

                               tGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGgtgt                                                       >  2:4231657/1‑93 (MQ=255)

                               tGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGgtgg                                                       >  2:4182574/1‑92 (MQ=255)

                                  tttATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGAt                                                    <  1:4352573/93‑1 (MQ=255)

                                   ttATCTGGCTTTAGCGATTGTTCTCATCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATg                                                   <  1:2145580/93‑1 (MQ=255)

                                     aTCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCa                                                 >  1:2057319/1‑93 (MQ=255)

                                      tCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCaa                                                <  1:3279459/93‑1 (MQ=255)

                                      tCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCaa                                                >  1:1004960/1‑93 (MQ=255)

                                         ggCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCt                                             <  2:1113313/93‑1 (MQ=255)

                                          gCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtt                                            <  1:3741948/93‑1 (MQ=255)

                                          gCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtt                                            <  2:1522953/93‑1 (MQ=255)

                                           cTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCttt                                           >  1:829325/1‑93 (MQ=255)

                                           cTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCttt                                           <  1:3776639/93‑1 (MQ=255)

                                            tttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtttt                                          <  1:861984/93‑1 (MQ=255)

                                            tttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtttt                                          <  1:731312/93‑1 (MQ=255)

                                                  gATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTa                                    <  1:535503/93‑1 (MQ=255)

                                                   aTTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTAc                                   <  1:4000963/93‑1 (MQ=255)

                                                   aTTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTAc                                   <  1:2360364/93‑1 (MQ=255)

                                                    ttGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACt                                  >  1:2222091/1‑93 (MQ=255)

                                                    ttGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACt                                  >  1:2608996/1‑93 (MQ=255)

                                                      gTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTg                                <  2:1524690/93‑1 (MQ=255)

                                                       ttCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGa                               <  2:727681/93‑1 (MQ=255)

                                                       ttCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGa                               <  1:3344298/93‑1 (MQ=255)

                                                        tctcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGAc                              <  2:3041289/93‑1 (MQ=255)

                                                        tctcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGAc                              >  1:1031800/1‑93 (MQ=255)

                                                         ctcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGAcc                             <  1:4661925/93‑1 (MQ=255)

                                                         ctcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGAcc                             >  1:3949671/1‑93 (MQ=255)

                                                          tcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCg                            >  1:2700222/1‑93 (MQ=255)

                                                           cAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGc                           <  2:2039560/93‑1 (MQ=255)

                                                            aaCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCt                          >  2:2635828/1‑93 (MQ=255)

                                                            aaCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCt                          <  1:3086658/93‑1 (MQ=255)

                                                            aaCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCt                          >  2:2179123/1‑93 (MQ=255)

                                                            aaCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCt                          <  2:3912871/93‑1 (MQ=255)

                                                             aCCTTTGGGACAGCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGTGTACTTGACCGCta                         <  1:401301/93‑1 (MQ=255)

                                                             aCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCta                         <  2:1291850/93‑1 (MQ=255)

                                                             aCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCta                         >  1:4232158/1‑93 (MQ=255)

                                                             aCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCta                         >  2:1884917/1‑93 (MQ=255)

                                                             aCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCta                         <  1:2519297/93‑1 (MQ=255)

                                                               cgttGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCtata                       <  2:2449671/91‑1 (MQ=255)

                                                               cTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCtata                       <  2:1643123/93‑1 (MQ=255)

                                                               cTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCtata                       <  2:3247236/93‑1 (MQ=255)

                                                                tttGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTatat                      <  2:4158253/93‑1 (MQ=255)

                                                                tttGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTatat                      <  1:1960524/93‑1 (MQ=255)

                                                                  tGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCg                    <  2:2304136/93‑1 (MQ=255)

                                                                   gggACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCgg                   <  2:3364858/93‑1 (MQ=255)

                                                                   gggACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCgg                   >  2:883714/1‑93 (MQ=255)

                                                                    ggACACCGTGCCGGCCCGCCGCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGt                  <  1:2793741/93‑1 (MQ=255)

                                                                    ggACACCGTGCCGGCCCGCCCCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGACCTTGCCCGCTCAATCgta                  >  1:1205773/1‑91 (MQ=255)

                                                                      acacCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaa                >  1:3007470/1‑93 (MQ=255)

                                                                       cacCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaaa               >  2:3615634/1‑93 (MQ=255)

                                                                        acCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaaaa              <  1:1378771/93‑1 (MQ=255)

                                                                        acCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaaaa              >  2:4730341/1‑93 (MQ=255)

                                                                        acCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaaaa              >  1:2827698/1‑93 (MQ=255)

                                                                        acCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaaaa              >  1:4387282/1‑93 (MQ=255)

                                                                        acCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTGCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTaaaa              <  1:4569115/93‑1 (MQ=255)

                                                                         ccGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAAc             <  2:2411205/93‑1 (MQ=255)

                                                                           gTGCCGGTCCGCCGCCGGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTa           <  1:1584802/93‑1 (MQ=255)

                                                                           gTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTa           >  2:1825550/1‑93 (MQ=255)

                                                                           gTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTa           <  1:299526/93‑1 (MQ=255)

                                                                           gTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTa           <  2:407875/93‑1 (MQ=255)

                                                                           gTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTa           <  2:4542236/93‑1 (MQ=255)

                                                                             gCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTAtt         <  2:3717415/93‑1 (MQ=255)

                                                                               cGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAACACTAtttc       >  1:4415566/1‑92 (MQ=255)

                                                                                ggcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTc      <  1:3396979/90‑1 (MQ=255)

                                                                                ggcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTc      <  2:596991/90‑1 (MQ=255)

                                                                                ggcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTc      <  2:1027677/90‑1 (MQ=255)

                                                                                ggcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTc      <  2:1596125/90‑1 (MQ=255)

                                                                                 gcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTca     >  1:3888814/3‑93 (MQ=255)

                                                                                 gcccgccaCCTGCGCGCCAGTTTGTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTca     <  2:661288/91‑1 (MQ=255)

                                                                                  cccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTcac    <  1:1366239/92‑1 (MQ=255)

                                                                                   ccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTcacc   >  2:4010841/1‑93 (MQ=255)

                                                                                    cgccacCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTcacca  <  2:697863/93‑1 (MQ=255)

GACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGTCCGCCGCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTCACCA  >  NC_000913/4484761‑4484937

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