BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,408,560	G→A	A144T (GCC→ACC)	yjfJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,408,560	0	G	A	99.7%	1075.8 / NA	301	A144T (GCC→ACC)	yjfJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (170/129); total (170/131)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 310 total aligned reads displayed.

CTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTA  >  NC_000913/4408470‑4408649

 tGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCa                                                                                        >  2:6190948/1‑93 (MQ=255)

 tGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCa                                                                                        >  1:397943/1‑93 (MQ=255)

  ggAAGTCTCCCGTGATGGCGTGGAAAAGGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCAc                                                                                       <  2:5657896/93‑1 (MQ=255)

  ggAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCAc                                                                                       >  1:4757466/1‑93 (MQ=255)

  ggAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCAc                                                                                       <  2:2475392/93‑1 (MQ=255)

  ggAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCAc                                                                                       <  2:1842615/93‑1 (MQ=255)

   gAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACt                                                                                      >  2:1046107/1‑93 (MQ=255)

   gAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACt                                                                                      >  2:3428839/1‑93 (MQ=255)

    aaGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTg                                                                                     >  2:4108868/1‑93 (MQ=255)

      gTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGTGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGaa                                                                                   >  1:3317403/1‑93 (MQ=255)

      gTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGaa                                                                                   >  2:3446012/1‑93 (MQ=255)

      gTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGaa                                                                                   >  2:4970404/1‑93 (MQ=255)

      gTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGaa                                                                                   >  2:2216902/1‑93 (MQ=255)

        ctcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGc                                                                                 >  1:5290358/1‑93 (MQ=255)

         tcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGcc                                                                                <  1:6018972/93‑1 (MQ=255)

         tcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGcc                                                                                >  2:2783817/1‑93 (MQ=255)

         tcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGcc                                                                                >  2:5840443/1‑93 (MQ=255)

           ccGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCAt                                                                              <  2:3733472/93‑1 (MQ=255)

            cGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATg                                                                             >  2:6030660/1‑93 (MQ=255)

                aTGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGc                                                                         <  2:4680579/93‑1 (MQ=255)

                  ggCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGt                                                                       >  2:3347035/1‑93 (MQ=255)

                   gCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTg                                                                      <  2:5267797/93‑1 (MQ=255)

                   gCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTg                                                                      <  1:3606413/93‑1 (MQ=255)

                   gCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTg                                                                      <  1:3138294/93‑1 (MQ=255)

                     gTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGca                                                                    >  1:4512159/1‑93 (MQ=255)

                      tGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGGTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGGGcac                                                                   >  1:5914513/1‑93 (MQ=255)

                        gAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCacaa                                                                 >  1:3677879/1‑93 (MQ=255)

                           aaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAg                                                              <  2:3268723/93‑1 (MQ=255)

                           aaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAg                                                              >  1:580716/1‑93 (MQ=255)

                            aaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAgg                                                             >  1:2958329/1‑93 (MQ=255)

                             aGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGc                                                            <  1:2570985/93‑1 (MQ=255)

                             aGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGc                                                            <  1:6666731/93‑1 (MQ=255)

                             aGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGc                                                            <  1:5061823/93‑1 (MQ=255)

                              gCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCt                                                           <  2:1615784/93‑1 (MQ=255)

                               cAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGCAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTg                                                          >  2:2061143/1‑93 (MQ=255)

                                   tACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaac                                                      <  2:3328133/93‑1 (MQ=255)

                                    acagcgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaaca                                                     >  2:682704/1‑93 (MQ=255)

                                     cagcgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaacaa                                                    <  2:4028273/93‑1 (MQ=255)

                                       gcgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAAcc                                                  <  1:5430691/93‑1 (MQ=255)

                                       gcgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAAcc                                                  <  2:3868031/93‑1 (MQ=255)

                                        cgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACct                                                 >  1:2462640/1‑93 (MQ=255)

                                         gacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctc                                                <  2:2913389/93‑1 (MQ=255)

                                         gacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctc                                                <  1:4134525/93‑1 (MQ=255)

                                         gacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctc                                                <  2:1373488/93‑1 (MQ=255)

                                          acagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctct                                               >  2:6153162/1‑93 (MQ=255)

                                           cagcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTa                                              <  2:2686892/93‑1 (MQ=255)

                                            agcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTAc                                             <  1:3247646/93‑1 (MQ=255)

                                            agcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTAc                                             <  2:5529793/93‑1 (MQ=255)

                                             gcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTAcc                                            >  1:4649894/1‑93 (MQ=255)

                                             gcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTAcc                                            >  1:99061/1‑93 (MQ=255)

                                               gcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGt                                          >  1:5611708/1‑93 (MQ=255)

                                                cagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGtt                                         >  1:766812/1‑93 (MQ=255)

                                                cagcgTATCGCTCAGGTTGGGCAGCAAATGGAAGTCCTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGTCTGTAGCAACCTCTACCGtt                                         >  1:1219253/1‑93 (MQ=255)

                                                 agcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTg                                        >  1:2363842/1‑93 (MQ=255)

                                                 agcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTg                                        >  1:2647411/1‑93 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       >  1:2811616/1‑93 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       <  1:2876210/93‑1 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       >  1:5843829/1‑93 (MQ=255)

                                                    gTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGgcg                                     >  2:308225/1‑93 (MQ=255)

                                                     tATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGgcgc                                    >  1:3192876/1‑93 (MQ=255)

                                                      aTCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCa                                   >  2:4799457/1‑93 (MQ=255)

                                                      aTCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCa                                   >  1:3397923/1‑93 (MQ=255)

                                                        cGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATc                                 >  1:1148957/1‑93 (MQ=255)

                                                        cGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATc                                 <  2:3947433/93‑1 (MQ=255)

                                                         gCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCt                                >  1:4381211/1‑93 (MQ=255)

                                                          cTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCtt                               <  2:4357468/93‑1 (MQ=255)

                                                           tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTc                              >  1:489116/1‑93 (MQ=255)

                                                           tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTc                              >  1:2193951/1‑93 (MQ=255)

                                                           tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTc                              <  1:6442429/93‑1 (MQ=255)

                                                           tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTc                              <  2:4287390/93‑1 (MQ=255)

                                                           tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTc                              >  2:1196149/1‑93 (MQ=255)

                                                           tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTc                              >  1:4272360/1‑93 (MQ=255)

                                                            cAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTcc                             <  1:1564544/93‑1 (MQ=255)

                                                            cAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTcc                             <  1:2024007/93‑1 (MQ=255)

                                                            cAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTcc                             <  2:5226805/93‑1 (MQ=255)

                                                            cAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTcc                             >  2:4715105/1‑93 (MQ=255)

                                                            cAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTcc                             >  1:1415070/1‑93 (MQ=255)

                                                             aGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCa                            >  2:5130118/1‑93 (MQ=255)

                                                             aGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCa                            <  1:1253209/93‑1 (MQ=255)

                                                              gTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAg                           <  2:2611508/93‑1 (MQ=255)

                                                              gTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAg                           >  1:3849524/1‑93 (MQ=255)

                                                                ttGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCg                         >  2:628340/1‑93 (MQ=255)

                                                                   agcagcAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTTCAGCGttt                      <  1:6571038/93‑1 (MQ=255)

                                                                   agcagcAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGATGGCGCATCTTCCAGCGttt                      >  1:3950568/1‑93 (MQ=255)

                                                                       gcaAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGAc                  >  2:5476485/1‑93 (MQ=255)

                                                                          aaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGc               >  2:4451512/1‑93 (MQ=255)

                                                                          aaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGc               >  1:2744284/1‑93 (MQ=255)

                                                                           aTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCg              <  1:6238252/93‑1 (MQ=255)

                                                                             ggAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGc            >  2:3541080/1‑93 (MQ=255)

                                                                                 gTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGaa        >  1:4037481/1‑93 (MQ=255)

                                                                                 gTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGaa        <  2:6399287/93‑1 (MQ=255)

                                                                                 gTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGaa        <  1:1851450/93‑1 (MQ=255)

                                                                                 gTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGaa        <  2:5997485/93‑1 (MQ=255)

                                                                                  tCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAAt       >  1:1036129/1‑93 (MQ=255)

                                                                                  tCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAAt       <  1:6207737/93‑1 (MQ=255)

                                                                                  tCGTTAAAACCACTGAAGCCATGAAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGCTTCGACAGCGGCAGAAt       >  1:512166/1‑93 (MQ=255)

                                                                                    gTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATcc     <  2:4527334/93‑1 (MQ=255)

                                                                                    gTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATcc     >  2:211748/1‑93 (MQ=255)

                                                                                     ttAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCt    <  1:6371428/93‑1 (MQ=255)

                                                                                     ttAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCt    >  2:379989/1‑93 (MQ=255)

CTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTA  >  NC_000913/4408470‑4408649

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: `‑`