BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,476,406	A→G	D124G (GAC→GGC)	yjgD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,476,406	0	A	G	80.8%	521.2 / 131.6	260	D124G (GAC→GGC)	yjgD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (29/21); new base G (98/112); total (127/133)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 270 total aligned reads displayed.

GTTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGA  >  NC_000913/4476315‑4476497

   gAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGAt                                                                                         >  1:1267777/1‑93 (MQ=255)

   gAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGAt                                                                                         >  2:6375850/1‑93 (MQ=255)

   gAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGAt                                                                                         >  2:4697397/1‑93 (MQ=255)

    aacaacTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATg                                                                                        <  1:6973917/93‑1 (MQ=255)

     acaacTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGa                                                                                       >  2:1033050/1‑93 (MQ=255)

     acaacTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGa                                                                                       >  2:4188849/1‑93 (MQ=255)

     acaacTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGa                                                                                       <  1:2207806/93‑1 (MQ=255)

     acaacTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGa                                                                                       <  1:2051304/93‑1 (MQ=255)

      caacTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGaa                                                                                      >  1:4018536/1‑93 (MQ=255)

        acTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGa                                                                                    >  2:2143856/1‑93 (MQ=255)

        acTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGa                                                                                    <  1:5428958/93‑1 (MQ=255)

         cTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGAt                                                                                   >  2:3222999/1‑93 (MQ=255)

          tgatgaCGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGAtt                                                                                  <  1:6064188/93‑1 (MQ=255)

            atgaCGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGAtttt                                                                                <  1:4620992/93‑1 (MQ=255)

             tgaCGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTg                                                                               <  2:6107081/93‑1 (MQ=255)

              gaCGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGt                                                                              >  1:2639121/1‑93 (MQ=255)

                  cTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGAt                                                                          <  2:1426206/93‑1 (MQ=255)

                   tGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATg                                                                         >  1:7104752/1‑93 (MQ=255)

                   tGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATg                                                                         <  1:3612287/93‑1 (MQ=255)

                   tGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATg                                                                         >  1:5497501/1‑93 (MQ=255)

                    ggCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATTGCGGCGATGAAGATTTTGTCGATGa                                                                        <  1:7443468/93‑1 (MQ=255)

                    ggCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGa                                                                        <  1:2845418/93‑1 (MQ=255)

                     gCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGaa                                                                       <  2:2006146/93‑1 (MQ=255)

                     gCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGaa                                                                       <  1:1939091/93‑1 (MQ=255)

                       agagAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGa                                                                     <  1:5877035/93‑1 (MQ=255)

                        gagaAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGAc                                                                    <  1:2504232/93‑1 (MQ=255)

                        gagaAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGAc                                                                    <  1:6647791/93‑1 (MQ=255)

                        gagaAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGAc                                                                    <  2:589689/93‑1 (MQ=255)

                        gagaAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGAc                                                                    >  1:5298268/1‑93 (MQ=255)

                        gagaAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGAc                                                                    >  2:4543523/1‑93 (MQ=255)

                          gaAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGa                                                                  <  1:115961/93‑1 (MQ=255)

                          gaAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGa                                                                  <  2:4828097/93‑1 (MQ=255)

                           aaaTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGAt                                                                 <  1:7229457/93‑1 (MQ=255)

                           aaaTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGAt                                                                 <  1:3760384/93‑1 (MQ=255)

                           aaaTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGAt                                                                 <  1:2001330/93‑1 (MQ=255)

                            aaTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATg                                                                >  1:2390944/1‑93 (MQ=255)

                             aTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGa                                                               <  2:5251632/93‑1 (MQ=255)

                             aTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGa                                                               <  2:3532948/93‑1 (MQ=255)

                              tttGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGAc                                                              >  2:5177158/1‑93 (MQ=255)

                                tGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACgg                                                            >  2:6218415/1‑93 (MQ=255)

                                 gACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGa                                                           <  1:6848165/93‑1 (MQ=255)

                                  aCGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAg                                                          >  1:5809497/1‑93 (MQ=255)

                                  aCGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAg                                                          <  2:7107727/93‑1 (MQ=255)

                                   cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGt                                                         >  2:1105453/1‑93 (MQ=255)

                                   cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGt                                                         >  1:1473635/1‑93 (MQ=255)

                                    gTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGtt                                                        >  1:5564536/1‑93 (MQ=255)

                                        aGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGcc                                                    <  2:768597/93‑1 (MQ=255)

                                         gTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCa                                                   <  2:2849974/93‑1 (MQ=255)

                                          tACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCAc                                                  <  2:915805/93‑1 (MQ=255)

                                          tACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCAc                                                  <  1:1895080/93‑1 (MQ=255)

                                            cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTa                                                >  2:1550120/1‑93 (MQ=255)

                                               cgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAAtt                                             >  1:5668211/1‑93 (MQ=255)

                                                ggATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTa                                            <  2:1867503/93‑1 (MQ=255)

                                                ggATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTa                                            >  1:6823167/1‑93 (MQ=255)

                                                 gATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTaa                                           <  2:7252177/93‑1 (MQ=255)

                                                 gATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTaa                                           <  2:1633922/93‑1 (MQ=255)

                                                 gATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTaa                                           <  2:6562369/93‑1 (MQ=255)

                                                  aTGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAg                                          <  1:4502701/93‑1 (MQ=255)

                                                   tGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGc                                         <  2:6019663/93‑1 (MQ=255)

                                                   tGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGc                                         >  2:5394872/1‑93 (MQ=255)

                                                      ggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTa                                      >  2:3960585/1‑93 (MQ=255)

                                                      ggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTa                                      >  2:3715420/1‑93 (MQ=255)

                                                       gCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAg                                     >  2:5730260/1‑93 (MQ=255)

                                                         acttacttTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAgcg                                   >  2:1508971/1‑93 (MQ=255)

                                                         acttacttTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAgcg                                   >  2:6805661/1‑93 (MQ=255)

                                                          cttacttTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAgcgc                                  >  1:4776083/1‑93 (MQ=255)

                                                          cttacttTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAgcgc                                  >  1:3051244/1‑93 (MQ=255)

                                                             acttTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATc                               <  1:3105937/93‑1 (MQ=255)

                                                              cttTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCa                              <  2:6036320/93‑1 (MQ=255)

                                                              cttTGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCa                              <  2:2675525/93‑1 (MQ=255)

                                                               tttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAg                             <  2:4483842/93‑1 (MQ=255)

                                                                 tGAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGc                           <  2:4971309/93‑1 (MQ=255)

                                                                  gAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGcc                          >  1:2900542/1‑93 (MQ=255)

                                                                  gAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGcc                          >  1:3498592/1‑93 (MQ=255)

                                                                  gAAGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGcc                          <  2:6299250/93‑1 (MQ=255)

                                                                    aGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCa                        <  2:2482155/93‑1 (MQ=255)

                                                                    aGATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCa                        <  2:2228587/93‑1 (MQ=255)

                                                                     gATCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAg                       <  2:526886/93‑1 (MQ=255)

                                                                      aTCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAgg                      <  2:135392/93‑1 (MQ=255)

                                                                       tCCTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACTATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                     >  1:3218076/1‑93 (MQ=255)

                                                                       tCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                     >  2:6930998/1‑93 (MQ=255)

                                                                        ccTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTa                    >  2:2407129/1‑93 (MQ=255)

                                                                         cTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAg                   >  1:7172085/1‑93 (MQ=255)

                                                                         cTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAg                   <  2:3415082/93‑1 (MQ=255)

                                                                         cTAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAg                   >  1:2946563/1‑93 (MQ=255)

                                                                          tAACGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGt                  >  1:389165/1‑93 (MQ=255)

                                                                          tAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGt                  <  1:3415118/93‑1 (MQ=255)

                                                                           aaCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGtt                 >  1:2365243/1‑93 (MQ=255)

                                                                            aCGGCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTc                <  2:5650609/93‑1 (MQ=255)

                                                                               gCGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTg             >  2:3151228/1‑93 (MQ=255)

                                                                               gCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTg             >  1:2112387/1‑93 (MQ=255)

                                                                                cGAAGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGGTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGc            >  1:4997752/1‑93 (MQ=255)

                                                                                   aGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAAc         <  1:6250716/93‑1 (MQ=255)

                                                                                   aGATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAAc         <  1:6409690/93‑1 (MQ=255)

                                                                                    gATGGCGGCGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAAcc        <  2:6822896/93‑1 (MQ=255)

                                                                                       ggcggcGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg     <  1:5222178/93‑1 (MQ=255)

                                                                                       ggcggcGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg     >  1:3338974/1‑93 (MQ=255)

                                                                                       ggCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg     >  2:1661533/1‑93 (MQ=255)

                                                                                        gcggcgATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGt    <  2:5147418/93‑1 (MQ=255)

                                                                                          gacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGa  <  1:6717636/93‑1 (MQ=255)

GTTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGA  >  NC_000913/4476315‑4476497

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