BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,035,428	A→G	noncoding (330/2905 nt)	rrlA →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,035,428	0	A	G	96.4%	243.3 / 0.3	111	noncoding (330/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (1/3); new base G (52/55); total (53/58)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 115 total aligned reads displayed.

CCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG  >  NC_000913/4035337‑4035506

ccAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGAc                                                                               <  1:3730789/93‑1 (MQ=33)

 cAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCa                                                                              <  2:3917974/93‑1 (MQ=21)

 cAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCa                                                                              >  1:612764/1‑93 (MQ=18)

 cAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCa                                                                              >  2:1287201/1‑93 (MQ=18)

  aGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGGGGAAGCGCCTGGAAAGGCGTGCGATACAGGGTGGCAg                                                                             >  2:3968991/1‑93 (MQ=255)

   gTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGc                                                                            >  2:1844579/1‑93 (MQ=18)

   gTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGc                                                                            >  1:3368796/1‑93 (MQ=18)

     aGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGccc                                                                          <  1:2611644/93‑1 (MQ=21)

     aGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGccc                                                                          <  1:4311816/93‑1 (MQ=21)

         gcgagcgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTa                                                                      <  2:4626614/93‑1 (MQ=21)

            agcgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTacac                                                                   >  1:3263316/1‑93 (MQ=18)

              cgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACaa                                                                 >  2:3380257/1‑93 (MQ=18)

                aaCGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACaaaa                                                               >  1:3860310/1‑93 (MQ=18)

                aaCGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACaaaa                                                               >  1:4516133/1‑93 (MQ=33)

                 aCGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACaaaaa                                                              <  1:1805216/93‑1 (MQ=21)

                  cGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAAt                                                             <  2:2943170/93‑1 (MQ=18)

                  cGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAAt                                                             <  1:2620325/93‑1 (MQ=18)

                     ggAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGca                                                          >  2:3026769/1‑93 (MQ=18)

                      gAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGcac                                                         >  2:595372/1‑93 (MQ=18)

                        gcagcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACAt                                                       <  1:2249429/93‑1 (MQ=18)

                        gcagcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACAt                                                       >  2:3746725/1‑93 (MQ=18)

                        gcagcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACAt                                                       <  1:2188627/93‑1 (MQ=21)

                        gcagcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACAt                                                       >  1:1904352/1‑93 (MQ=18)

                        gcagcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACAt                                                       <  2:691657/93‑1 (MQ=21)

                          agcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGc                                                     >  1:301748/1‑93 (MQ=18)

                          agcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGc                                                     <  1:4252424/93‑1 (MQ=21)

                          agcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGc                                                     <  1:414710/93‑1 (MQ=18)

                          agcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGc                                                     >  2:3804266/1‑93 (MQ=18)

                          agcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGc                                                     >  2:812998/1‑93 (MQ=18)

                          agcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGc                                                     <  2:1507318/93‑1 (MQ=21)

                          agcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGc                                                     <  2:4040838/93‑1 (MQ=33)

                           gcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCt                                                    >  1:1506608/1‑93 (MQ=18)

                            cccAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCtg                                                   >  2:4112840/1‑93 (MQ=18)

                            cccAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCtg                                                   >  1:3298661/1‑93 (MQ=18)

                            cccAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCtg                                                   <  1:2905996/93‑1 (MQ=21)

                             ccAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCtgt                                                  <  1:2067062/93‑1 (MQ=18)

                              cAGAGTCTGAGCCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTCACAGCCCCGTACACAAAAATGCACATGCtgtg                                                 <  1:2128055/93‑1 (MQ=11)

                               agagCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGa                                                <  1:606060/93‑1 (MQ=21)

                               agagCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGa                                                >  2:4175408/1‑93 (MQ=18)

                               agagCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGa                                                <  2:4413985/93‑1 (MQ=21)

                                 agCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGc                                              >  2:1892888/1‑93 (MQ=18)

                                 agCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGc                                              <  2:2772259/93‑1 (MQ=18)

                                  gCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCt                                             <  1:2050925/93‑1 (MQ=18)

                                  gCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCt                                             <  1:3871755/93‑1 (MQ=21)

                                   ccTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTc                                            <  2:2013779/93‑1 (MQ=18)

                                    cTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCCCg                                           >  2:2673808/1‑93 (MQ=16)

                                     tGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGa                                          <  1:40729/93‑1 (MQ=18)

                                     tGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGa                                          >  1:1358865/1‑93 (MQ=21)

                                     tGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGa                                          >  2:3407761/1‑93 (MQ=18)

                                       aaTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATg                                        >  2:1906164/1‑93 (MQ=21)

                                       aaTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATg                                        <  1:76147/93‑1 (MQ=18)

                                       aaTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATg                                        <  1:57293/93‑1 (MQ=18)

                                       aaTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATg                                        <  1:3816762/93‑1 (MQ=18)

                                       aaTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATg                                        <  2:855254/93‑1 (MQ=18)

                                        aTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGa                                       >  1:102037/1‑93 (MQ=18)

                                        aTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGa                                       <  2:2199086/93‑1 (MQ=18)

                                        aTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGa                                       >  1:2431468/1‑93 (MQ=18)

                                           aGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTa                                    <  2:2961325/93‑1 (MQ=18)

                                           aGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTa                                    <  2:3790112/93‑1 (MQ=18)

                                           aGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTa                                    >  2:1696249/1‑93 (MQ=18)

                                            gtgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAg                                   >  2:1281201/1‑93 (MQ=18)

                                            gtgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAg                                   <  2:342776/93‑1 (MQ=18)

                                            gtgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAg                                   >  2:3828302/1‑93 (MQ=18)

                                             tgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAgg                                  <  2:3417554/93‑1 (MQ=18)

                                             tgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAgg                                  <  2:111556/93‑1 (MQ=18)

                                             tgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAgg                                  <  2:1411776/93‑1 (MQ=18)

                                             tgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAgg                                  <  1:786971/93‑1 (MQ=18)

                                             tgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAgg                                  >  1:2219785/1‑93 (MQ=18)

                                              gtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAggg                                 <  2:1388412/93‑1 (MQ=18)

                                               tgtgtgtTAGTGGAAGCGTCTTGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGc                                <  1:4596769/93‑1 (MQ=16)

                                               tgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGc                                <  1:4596772/93‑1 (MQ=18)

                                                  gtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCggg                             >  2:1231866/1‑93 (MQ=21)

                                                     ttAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGaca                          >  2:1183453/1‑93 (MQ=18)

                                                      tAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGacac                         >  2:2556741/1‑93 (MQ=18)

                                                       aGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACg                        <  1:1910752/93‑1 (MQ=18)

                                                        gTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGt                       <  2:3853072/93‑1 (MQ=18)

                                                          ggAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTgg                     >  1:1816236/1‑93 (MQ=18)

                                                           gAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGt                    >  2:4360913/1‑93 (MQ=18)

                                                           gAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGt                    >  2:1788644/1‑93 (MQ=21)

                                                            aaGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTa                   <  2:2900550/93‑1 (MQ=18)

                                                             aGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTac                  >  2:1270964/1‑92 (MQ=18)

                                                             aGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTAt                  >  2:1507347/1‑93 (MQ=18)

                                                              gCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATc                 <  1:391318/93‑1 (MQ=18)

                                                               cGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATcc                <  1:2224247/93‑1 (MQ=18)

                                                                 tCTGGAAAGTCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTg              >  2:3143305/1‑93 (MQ=16)

                                                                 tCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTg              >  2:827961/1‑93 (MQ=18)

                                                                  cTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGt             <  1:1482340/93‑1 (MQ=18)

                                                                  cTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGt             <  1:4055449/93‑1 (MQ=18)

                                                                   tGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTc            <  1:1716457/93‑1 (MQ=18)

                                                                    ggAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCt           >  1:3848439/1‑93 (MQ=18)

                                                                    ggAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCt           >  1:3848440/1‑93 (MQ=18)

                                                                    ggAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCt           >  2:3881830/1‑93 (MQ=18)

                                                                     gAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTg          >  1:1062644/1‑93 (MQ=18)

                                                                      aaaGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGa         >  2:4500448/1‑93 (MQ=21)

                                                                         ggCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAata      >  1:4375232/1‑93 (MQ=18)

                                                                         ggCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAata      <  2:1052096/93‑1 (MQ=18)

                                                                         ggCGTGCGATACAGGGTGGCACCCCCGTACAAAAAAATGCACATGCTGTGAGCTCGATGAGCAGGGCGGGACACGTGGTATCCTGTCTGAata      >  1:3797100/1‑93 (MQ=11)

                                                                            gTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgg   <  1:4468979/93‑1 (MQ=17)

                                                                            gTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgg   <  2:2898824/93‑1 (MQ=17)

                                                                             tGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATggg  <  1:3384141/93‑1 (MQ=14)

CCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG  >  NC_000913/4035337‑4035506

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