BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,173,983	C→T	G136R (GGG→AGG)	ycfT ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,173,983	0	C	T	99.5%	791.9 / NA	211	G136R (GGG→AGG)	ycfT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (96/114); minor base A (1/0); total (97/114)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

GCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCAC  >  NC_000913/1173893‑1174072

gCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGt                                                                                         <  2:4007626/93‑1 (MQ=255)

gCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGt                                                                                         <  1:4623253/93‑1 (MQ=255)

gCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGt                                                                                         <  2:2102624/93‑1 (MQ=255)

gCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGt                                                                                         <  1:87814/93‑1 (MQ=255)

gCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGt                                                                                         <  2:5016129/93‑1 (MQ=255)

 cAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCtgtg                                                                                        <  1:5641204/93‑1 (MQ=255)

  aGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGc                                                                                       <  1:3390132/93‑1 (MQ=255)

  aGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGc                                                                                       >  2:1437186/1‑93 (MQ=255)

   gcgcCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCa                                                                                      <  1:4268554/93‑1 (MQ=255)

   gcgcCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCa                                                                                      >  1:192677/1‑93 (MQ=255)

   gcgcCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCa                                                                                      >  1:5973670/1‑93 (MQ=255)

      ccAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGa                                                                                   >  2:5789749/1‑93 (MQ=255)

       cAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGaa                                                                                  >  1:3036017/1‑93 (MQ=255)

        aGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAAc                                                                                 <  2:704199/93‑1 (MQ=255)

        aGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAAc                                                                                 >  1:127341/1‑93 (MQ=255)

         gACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACt                                                                                >  1:4356615/1‑93 (MQ=255)

         gACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACt                                                                                >  1:5386044/1‑93 (MQ=255)

         gACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACt                                                                                <  1:1223310/93‑1 (MQ=255)

          aCGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTc                                                                               <  1:4969370/93‑1 (MQ=255)

           cGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCa                                                                              <  2:3035326/93‑1 (MQ=255)

             gCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCAcc                                                                            >  2:2072368/1‑93 (MQ=255)

             gCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCAcc                                                                            >  1:45107/1‑93 (MQ=255)

             gCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCAcc                                                                            >  2:4471449/1‑93 (MQ=255)

               tAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCgg                                                                          >  1:5066552/1‑93 (MQ=255)

               tAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCgg                                                                          >  2:5320307/1‑93 (MQ=255)

                aaaaaTTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGt                                                                         <  2:2780901/93‑1 (MQ=255)

                aaaaaTTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGt                                                                         >  1:954421/1‑93 (MQ=255)

                   aaTTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGCTGGa                                                                      >  2:337520/1‑93 (MQ=255)

                    aTTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGaa                                                                     <  1:4205007/93‑1 (MQ=255)

                     ttttACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAAt                                                                    <  1:4019046/93‑1 (MQ=255)

                       ttACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCg                                                                  <  2:1809777/93‑1 (MQ=255)

                            tACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATa                                                             <  2:1331694/93‑1 (MQ=255)

                            tACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATa                                                             <  1:675726/93‑1 (MQ=255)

                             aCCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATaa                                                            <  1:5216313/93‑1 (MQ=255)

                             aCCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATaa                                                            <  1:1758605/93‑1 (MQ=255)

                             aCCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATaa                                                            <  2:658663/93‑1 (MQ=255)

                              ccACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAg                                                           <  2:4422672/93‑1 (MQ=255)

                              ccACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAg                                                           <  1:2196757/93‑1 (MQ=255)

                               cACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGc                                                          >  1:2284557/1‑93 (MQ=255)

                                aCGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCg                                                         >  2:1577523/1‑93 (MQ=255)

                                  gAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGc                                                       >  2:4556350/1‑93 (MQ=255)

                                  gAAATAGACAATTAACGCATACAGATACCACAAGCAGGTGCTGGCGGCGATCATCCAGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGc                                                       >  1:6309473/1‑93 (MQ=255)

                                   aaaTAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCa                                                      <  2:1651688/93‑1 (MQ=255)

                                   aaaTAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGGGGAATCGGCATAAGCGGCa                                                      >  1:1516226/1‑93 (MQ=255)

                                    aaTAGACAATTAACGCATACAGATACCACAAGCTTGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCAt                                                     >  1:2411634/1‑93 (MQ=255)

                                    aaTAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCAt                                                     >  2:2202600/1‑93 (MQ=255)

                                    aaTAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCAt                                                     <  2:5746825/93‑1 (MQ=255)

                                    aaTAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCAt                                                     <  2:810977/93‑1 (MQ=255)

                                    aaTAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCAt                                                     <  1:4751333/93‑1 (MQ=255)

                                    aaTAGACAATAAACGCATACAGATAACACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCAt                                                     <  2:1093946/93‑1 (MQ=255)

                                     aTAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCAtt                                                    >  1:2830235/1‑93 (MQ=255)

                                      tAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTg                                                   <  1:1748343/93‑1 (MQ=255)

                                       aGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTgg                                                  <  1:5434421/93‑1 (MQ=255)

                                        gACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGa                                                 <  2:5440107/93‑1 (MQ=255)

                                         aCAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATCAGCGGCATTGGAg                                                <  2:6051659/93‑1 (MQ=255)

                                             ttAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCAt                                            >  1:1120990/1‑93 (MQ=255)

                                                  gCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCtt                                       >  2:3806888/1‑93 (MQ=255)

                                                   cATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTa                                      <  2:1163871/93‑1 (MQ=255)

                                                   cATACAGATACCACAAGCTGGTGCTGGCGGCGATCATCCTGTGCACGAACTCACCGGTGAAATCGGCATAAGACGCATTGGAGGCATTGCGTa                                      >  2:3499045/1‑93 (MQ=255)

                                                    aTACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTaa                                     >  1:2108782/1‑93 (MQ=255)

                                                     tACAGTTACCACAAGCTGGTGCTGTCGGTGATCATCCGGTGCAGGAACTCACCGTTGGAATCGTCATAAGCGGCATTGGATGCATTGCTTaaa                                    <  2:3345822/93‑1 (MQ=255)

                                                     tACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTaaa                                    >  1:4796560/1‑93 (MQ=255)

                                                     tACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTaaa                                    >  1:3653832/1‑93 (MQ=255)

                                                       cAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATc                                  <  2:5912756/93‑1 (MQ=255)

                                                       cAGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATc                                  <  1:3010358/93‑1 (MQ=255)

                                                        aGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATcg                                 >  2:1545444/1‑93 (MQ=255)

                                                        aGATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATcg                                 >  2:4919005/1‑93 (MQ=255)

                                                         gATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATcgg                                >  2:1849908/1‑92 (MQ=255)

                                                         gATACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATcgc                                >  1:3826789/1‑93 (MQ=255)

                                                          aTACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATcgcg                               <  2:2324821/93‑1 (MQ=255)

                                                           tACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCgcgc                              >  2:3643113/1‑93 (MQ=255)

                                                           tACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCgcgc                              <  2:4071781/93‑1 (MQ=255)

                                                           tACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCgcgc                              >  2:311842/1‑93 (MQ=255)

                                                           tACCACAAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCgcgc                              >  1:5752505/1‑93 (MQ=255)

                                                              cacaAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCa                           <  1:1562847/93‑1 (MQ=255)

                                                              cacaAGCGGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCa                           <  2:1154774/93‑1 (MQ=255)

                                                               acaAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAg                          >  1:6147975/1‑93 (MQ=255)

                                                               acaAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAg                          <  2:5852063/93‑1 (MQ=255)

                                                               acaAGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAg                          <  1:2730185/93‑1 (MQ=255)

                                                                 aaGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGt                        <  1:1768249/93‑1 (MQ=255)

                                                                  aGCTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTg                       >  1:4440945/1‑93 (MQ=255)

                                                                    cTGGTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTgcc                     >  2:4449304/1‑93 (MQ=255)

                                                                      ggTGCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTgccag                   >  1:2009266/1‑93 (MQ=255)

                                                                         gCTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTgccagcca                <  1:3067688/93‑1 (MQ=255)

                                                                          cTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCAc               <  2:2124102/93‑1 (MQ=255)

                                                                          cTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCAc               <  2:5105485/93‑1 (MQ=255)

                                                                          cTGGCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCAc               <  2:2357977/93‑1 (MQ=255)

                                                                            ggCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTg             <  2:3708965/93‑1 (MQ=255)

                                                                             gCGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTgg            <  2:2525811/93‑1 (MQ=255)

                                                                              cGGTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGt           >  2:5674528/1‑93 (MQ=255)

                                                                               ggTGATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGtt          >  2:4588282/1‑93 (MQ=255)

                                                                                  gATCATCCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCAACTGGTTCAg       >  1:1348717/1‑93 (MQ=255)

                                                                                   atcatcCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAgc      >  1:4386651/1‑93 (MQ=255)

                                                                                     catcCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAgcgc    <  1:3758516/93‑1 (MQ=255)

                                                                                     catcCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAgcgc    <  1:1438748/93‑1 (MQ=255)

                                                                                      atcCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGGCAGCCACTGGTTCAGCGCa   >  2:4181643/1‑93 (MQ=255)

                                                                                       tcCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCAc  <  1:4077313/93‑1 (MQ=255)

                                                                                       tcCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCAc  >  2:499209/1‑93 (MQ=255)

                                                                                       tcCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCAc  >  2:2886417/1‑93 (MQ=255)

                                                                                       tcCTGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCAc  >  1:6186022/1‑93 (MQ=255)

GCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCAC  >  NC_000913/1173893‑1174072

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