BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,507,420	T→C	T18A (ACG→GCG)	b1437 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,507,420	0	T	C	99.7%	1038.0 / ‑5.1	326	T18A (ACG→GCG)	b1437	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (135/190); total (136/190)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 343 total aligned reads displayed.

ACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGTCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAATGA  >  NC_000913/1507330‑1507509

  gTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCAcc                                                                                       >  2:648635/1‑93 (MQ=255)

    tGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCAcccc                                                                                     <  2:6110383/93‑1 (MQ=255)

       cTCGAATCTCGGGGCGTCGCTGTAGCGAATGACAGCAACCATATGAAACTCAGATTTCATGGGAGGCGCAGTGTCATGCAGCGCCACCCCTGc                                                                                  <  2:8955730/93‑1 (MQ=255)

         cGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGa                                                                                <  2:3383944/93‑1 (MQ=255)

         cGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGa                                                                                <  1:7286283/93‑1 (MQ=255)

         cGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGa                                                                                >  1:1677462/1‑93 (MQ=255)

         cGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGa                                                                                >  1:6223424/1‑93 (MQ=255)

         cGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGa                                                                                >  1:7193586/1‑93 (MQ=255)

         cGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGa                                                                                <  1:982163/93‑1 (MQ=255)

           aaTCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATg                                                                              <  1:7844343/93‑1 (MQ=255)

           aaTCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATg                                                                              >  2:457930/1‑93 (MQ=255)

             tctcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATgag                                                                            <  1:1895846/93‑1 (MQ=255)

             tctcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATgag                                                                            <  2:4482123/93‑1 (MQ=255)

             tctcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATgag                                                                            >  1:5383913/1‑93 (MQ=255)

              ctcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATgaga                                                                           >  2:1863237/1‑93 (MQ=255)

              ctcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATgaga                                                                           <  2:1927340/93‑1 (MQ=255)

                cAGGGCGTCGATGTAGAGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTAATGCCGCGCCACCCCTGCGATGAGAtt                                                                         <  2:8821508/93‑1 (MQ=255)

                  gggCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTaa                                                                       >  2:1370953/1‑93 (MQ=255)

                         gATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCa                                                                <  2:1509337/93‑1 (MQ=255)

                         gATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCa                                                                >  2:9569784/1‑93 (MQ=255)

                         gATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCa                                                                <  1:5672759/93‑1 (MQ=255)

                            gTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTg                                                             <  1:1080008/93‑1 (MQ=255)

                            gTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTg                                                             <  2:4867275/93‑1 (MQ=255)

                             tAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGc                                                            >  1:5055905/1‑93 (MQ=255)

                             tAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGc                                                            >  2:5922906/1‑93 (MQ=255)

                             tAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGc                                                            >  1:6518796/1‑93 (MQ=255)

                                cGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTa                                                         <  1:9547933/93‑1 (MQ=255)

                                cGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTa                                                         >  2:4216698/1‑93 (MQ=255)

                                  aaTGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTaaa                                                       >  2:5725220/1‑93 (MQ=255)

                                  aaTGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTaaa                                                       >  2:8685354/1‑93 (MQ=255)

                                  aaTGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCAACCCCTGCGATGAGATTAAAGAACCATTGCGTaaa                                                       <  2:4065550/93‑1 (MQ=255)

                                   aTGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAg                                                      <  1:3372282/93‑1 (MQ=255)

                                   aTGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAg                                                      >  2:9137524/1‑93 (MQ=255)

                                    tGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGc                                                     <  1:5936661/93‑1 (MQ=255)

                                    tGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGc                                                     <  2:2862622/93‑1 (MQ=255)

                                      gcagcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCaa                                                   <  2:4886604/93‑1 (MQ=255)

                                        agcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATc                                                 <  1:4917597/93‑1 (MQ=255)

                                        agcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATc                                                 <  2:8012198/93‑1 (MQ=255)

                                        agcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATc                                                 >  1:46147/1‑93 (MQ=255)

                                         gcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGTCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATcc                                                >  1:9197478/1‑93 (MQ=255)

                                         gcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATcc                                                <  1:3851956/93‑1 (MQ=255)

                                            aCCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGa                                             >  1:7343069/1‑93 (MQ=255)

                                            aCCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAACAACCATTGCGTAAAGCAATCCTGa                                             <  2:7701720/93‑1 (MQ=255)

                                             ccATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGaa                                            >  1:4868349/1‑93 (MQ=255)

                                             ccATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGaa                                            <  1:4533107/93‑1 (MQ=255)

                                               aTTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAaac                                          <  2:8063827/93‑1 (MQ=255)

                                                  tGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAaacaac                                       >  1:8032902/1‑93 (MQ=255)

                                                  tGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAaacaac                                       <  1:9042062/93‑1 (MQ=255)

                                                   gAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACt                                      <  2:7034470/93‑1 (MQ=255)

                                                   gAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACt                                      <  2:6134729/93‑1 (MQ=255)

                                                   gAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACt                                      <  2:552550/93‑1 (MQ=255)

                                                    aaaCTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTc                                     <  2:8511806/93‑1 (MQ=255)

                                                      aCTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCgg                                   <  2:3791924/93‑1 (MQ=255)

                                                       cTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGt                                  <  1:2569803/93‑1 (MQ=255)

                                                       cTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGt                                  >  2:10009307/1‑93 (MQ=255)

                                                            gTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAg                             <  1:4123421/93‑1 (MQ=255)

                                                            gTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAg                             <  1:1498163/93‑1 (MQ=255)

                                                             tttCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGt                            <  1:7453830/93‑1 (MQ=255)

                                                             tttCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGt                            >  1:7892132/1‑93 (MQ=255)

                                                              ttCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGtt                           <  1:5850897/93‑1 (MQ=255)

                                                               tCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTa                          >  1:2237385/1‑93 (MQ=255)

                                                               tCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTa                          >  1:7728123/1‑93 (MQ=255)

                                                               tCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTa                          <  1:1397177/93‑1 (MQ=255)

                                                               tCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTa                          >  1:6314952/1‑93 (MQ=255)

                                                                taTGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTaa                         <  1:2033361/92‑1 (MQ=255)

                                                                cATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTaa                         >  1:298544/1‑93 (MQ=255)

                                                                 aTGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAAt                        >  2:5027752/1‑93 (MQ=255)

                                                                  tGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATc                       >  1:8522215/1‑93 (MQ=255)

                                                                     gAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGcc                    <  2:5612250/93‑1 (MQ=255)

                                                                     gAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGcc                    <  2:7213571/93‑1 (MQ=255)

                                                                       ggCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCaa                  >  2:705703/1‑93 (MQ=255)

                                                                       ggCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCaa                  <  2:5448557/93‑1 (MQ=255)

                                                                       ggCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCaa                  >  1:6097435/1‑93 (MQ=255)

                                                                       ggCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCaa                  <  2:6598543/93‑1 (MQ=255)

                                                                       ggCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCaa                  <  1:8222773/93‑1 (MQ=255)

                                                                        gcgcAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAAt                 <  2:217165/93‑1 (MQ=255)

                                                                         tgcAGTGTCATGCGGCGCCTCCCCTGCGATGAGATTAAAGAGCCATTGCTTAAAACAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAAtt                <  2:7573936/92‑1 (MQ=255)

                                                                         ggcAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAAtt                <  1:7149436/92‑1 (MQ=255)

                                                                         cgcAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAAtt                <  1:5566498/93‑1 (MQ=255)

                                                                         cgcAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAAtt                <  1:400575/93‑1 (MQ=255)

                                                                          gcAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTa               >  1:2741969/1‑93 (MQ=255)

                                                                           cAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaa              >  1:1038714/1‑93 (MQ=255)

                                                                             gtgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTcaaa            >  1:3321958/1‑93 (MQ=255)

                                                                             gtgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaaa            <  1:6586542/93‑1 (MQ=255)

                                                                             gtgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaaa            <  1:6945337/93‑1 (MQ=255)

                                                                              tgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGACACAACTCGGTTTGAGTTAATCGCCAATTaaaaa           >  1:9684819/1‑93 (MQ=255)

                                                                              tgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaaaa           >  2:2609278/1‑93 (MQ=255)

                                                                              tgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaaaa           <  1:3667603/93‑1 (MQ=255)

                                                                              tgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaaaa           >  2:3546869/1‑93 (MQ=255)

                                                                                tCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAgg         >  1:3082460/1‑93 (MQ=255)

                                                                                 cATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGt        <  2:9470297/93‑1 (MQ=255)

                                                                                  aTGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGtt       <  1:8377017/93‑1 (MQ=255)

                                                                                   tGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTa      >  1:8881729/1‑93 (MQ=255)

                                                                                   tGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTa      <  1:5794218/93‑1 (MQ=255)

                                                                                    gCCGCGCCACCCCTGCGATGAGATTAAAGACCCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTaa     <  1:4341429/93‑1 (MQ=255)

                                                                                    gCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTaa     <  2:2427525/93‑1 (MQ=255)

                                                                                     ccGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAAt    <  1:9920240/93‑1 (MQ=255)

                                                                                      cgcgCCACCCCTGCGATGGGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTTATCGCCAATTAAAAAGGTTAATg   <  2:9692496/93‑1 (MQ=255)

                                                                                       gcgcCACCCCTGCGATGAGATTAAATAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAATGa  <  2:8225080/93‑1 (MQ=255)

                                                                                       gcgcCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAATGa  <  2:3044783/93‑1 (MQ=255)

ACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGTCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAATGA  >  NC_000913/1507330‑1507509

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