BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,160,251	A→G	T145T (ACT→ACC)	sufI ←	suppressor of ftsI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,160,251	0	A	G	99.7%	911.4 / NA	289	T145T (ACT→ACC)	sufI	suppressor of ftsI
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (150/138); total (150/139)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 297 total aligned reads displayed.

GTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGG  >  NC_000913/3160161‑3160340

 ttGGGGATAGTCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTa                                                                                        <  1:6991172/93‑1 (MQ=255)

 ttGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTa                                                                                        <  2:29286/93‑1 (MQ=255)

 ttGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTa                                                                                        <  2:4925252/93‑1 (MQ=255)

 ttGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGGTGGGGGTa                                                                                        >  1:1693575/1‑93 (MQ=255)

 gtGGGGATAGGCAGCGACTTGCGGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTa                                                                                        <  1:5375316/92‑1 (MQ=255)

  tGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTAt                                                                                       >  2:2840756/1‑93 (MQ=255)

      gATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGc                                                                                   <  1:5566278/93‑1 (MQ=255)

       aTAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCg                                                                                  >  1:2865594/1‑93 (MQ=255)

        tAGGCAGCGCCTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGt                                                                                 <  2:1963185/93‑1 (MQ=255)

        tAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGt                                                                                 >  1:4827869/1‑93 (MQ=255)

        tAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGt                                                                                 >  2:7198002/1‑93 (MQ=255)

        tAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGt                                                                                 <  2:7874489/93‑1 (MQ=255)

         aGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGGATTGGCGTg                                                                                >  2:6536056/1‑93 (MQ=255)

         aGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGGTGGGGGGATTGGCGTg                                                                                >  1:7636653/1‑93 (MQ=255)

          ggCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGa                                                                               <  2:5597765/93‑1 (MQ=255)

          ggCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGa                                                                               >  2:2635040/1‑93 (MQ=255)

            cAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATa                                                                             >  1:3095527/1‑93 (MQ=255)

            cAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATa                                                                             <  1:1393217/93‑1 (MQ=255)

                gACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACcac                                                                         <  2:5772661/93‑1 (MQ=255)

                gACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACcac                                                                         >  2:5816411/1‑93 (MQ=255)

                gACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACcac                                                                         <  2:1749941/93‑1 (MQ=255)

                 cattGCGGACTTCCTCTACCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACcaca                                                                        <  1:799326/91‑1 (MQ=255)

                 aCTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACcaca                                                                        >  1:1802712/1‑93 (MQ=255)

                 aCTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACcaca                                                                        >  2:4603982/1‑93 (MQ=255)

                  cTTGCTGACTTCATCTTCCCCCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACag                                                                       <  1:2123506/93‑1 (MQ=255)

                  cTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACag                                                                       >  2:1452501/1‑93 (MQ=255)

                  cTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACag                                                                       >  1:105782/1‑93 (MQ=255)

                  cTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGGATTGGCGTGATACCACag                                                                       >  1:6576854/1‑93 (MQ=255)

                  cTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGCTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACag                                                                       >  2:6300825/1‑93 (MQ=255)

                   ttGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACaga                                                                      >  2:5920877/1‑93 (MQ=255)

                    tGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACagag                                                                     >  1:7767563/1‑93 (MQ=255)

                     gCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGCCCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGt                                                                    <  1:5219065/93‑1 (MQ=255)

                       tGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAg                                                                  >  2:3459649/1‑93 (MQ=255)

                         aCTTCATCTCCCCCCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCt                                                                <  2:4364631/93‑1 (MQ=255)

                               tCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTc                                                          >  1:665852/1‑93 (MQ=255)

                               tCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTc                                                          >  2:2191116/1‑93 (MQ=255)

                               tCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTc                                                          >  1:7216848/1‑93 (MQ=255)

                                cTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCt                                                         <  1:191484/93‑1 (MQ=255)

                                 ttCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTg                                                        >  2:766579/1‑93 (MQ=255)

                                 ttCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTg                                                        >  1:2463757/1‑93 (MQ=255)

                                  tCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTgg                                                       <  1:6778353/93‑1 (MQ=255)

                                      ccaGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGaa                                                   >  1:2769062/1‑93 (MQ=255)

                                      ccaGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGaa                                                   <  1:3317062/93‑1 (MQ=255)

                                       caGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCAAAGAGTAGCTGCGTTCTGGCGAAt                                                  <  1:4132683/93‑1 (MQ=255)

                                        aGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATg                                                 <  2:6603157/93‑1 (MQ=255)

                                         gCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATgg                                                <  2:1692741/93‑1 (MQ=255)

                                          ccACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATggg                                               >  1:5720889/1‑93 (MQ=255)

                                           cacaTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGc                                              <  1:2624075/93‑1 (MQ=255)

                                           cacaTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGc                                              >  2:965658/1‑93 (MQ=255)

                                             caTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAg                                            >  2:7504715/1‑93 (MQ=255)

                                               ttCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTa                                          >  1:4248621/1‑93 (MQ=255)

                                               ttCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTa                                          <  1:7475123/93‑1 (MQ=255)

                                               ttCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTa                                          >  2:4167298/1‑93 (MQ=255)

                                               ttCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTa                                          <  1:6210967/93‑1 (MQ=255)

                                               ttCCGGCAAGGCCGTTATAGACCCGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTa                                          >  1:2591484/1‑93 (MQ=255)

                                               ttCCGGCAAGCCCGTTATAGCCCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTa                                          <  2:6226865/93‑1 (MQ=255)

                                                   ggCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACggg                                      >  2:6688790/1‑93 (MQ=255)

                                                   ggCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACggg                                      <  1:6417759/93‑1 (MQ=255)

                                                   ggCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACggg                                      >  2:5636832/1‑93 (MQ=255)

                                                   ggCAAGGCCGTGATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACggg                                      <  1:2009893/93‑1 (MQ=255)

                                                    gCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGt                                     <  2:68499/93‑1 (MQ=255)

                                                    gCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGt                                     <  2:7405227/93‑1 (MQ=255)

                                                    gCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGt                                     <  2:2304523/93‑1 (MQ=255)

                                                    gCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGt                                     <  2:3841109/93‑1 (MQ=255)

                                                      aaGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGc                                   >  2:780224/1‑93 (MQ=255)

                                                       aGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGcc                                  >  2:7058320/1‑93 (MQ=255)

                                                         gCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCa                                <  2:3685345/93‑1 (MQ=255)

                                                           cGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGt                              >  2:5022935/1‑93 (MQ=255)

                                                            gTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTc                             <  2:3122901/93‑1 (MQ=255)

                                                            gTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTc                             <  1:2065787/93‑1 (MQ=255)

                                                             ttATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCa                            >  1:4273025/1‑93 (MQ=255)

                                                             ttATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCa                            >  2:2350773/1‑93 (MQ=255)

                                                             ttATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCa                            >  1:3769641/1‑93 (MQ=255)

                                                             ttATAGACCTGCTGAGCAGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCa                            >  2:845544/1‑93 (MQ=255)

                                                              tataGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAg                           >  2:1224964/1‑93 (MQ=255)

                                                                taGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCg                         <  1:5882296/93‑1 (MQ=255)

                                                                  gACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGTCGAATGGGCAGTACGGGTGCCCAGTCAGCGtt                       >  1:4981621/1‑93 (MQ=255)

                                                                  gACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGtt                       <  2:3402110/93‑1 (MQ=255)

                                                                   aCCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGttt                      <  1:95684/93‑1 (MQ=255)

                                                                   aCCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGttt                      <  2:3309034/93‑1 (MQ=255)

                                                                     ctgctgAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCAAATGGGCAGTACGGGTGCCCAGTCAGCGTTTgg                    <  2:2746910/93‑1 (MQ=255)

                                                                      tgctgAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGc                   <  2:3165178/93‑1 (MQ=255)

                                                                      tgctgAGCCGTGCGGTTGGGGGGATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGc                   <  2:4881824/93‑1 (MQ=255)

                                                                        ctgAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGa                 <  1:4515298/93‑1 (MQ=255)

                                                                         tgAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAc                >  1:1606653/1‑93 (MQ=255)

                                                                          gAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAca               <  1:4625022/93‑1 (MQ=255)

                                                                            gCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAcatc             <  2:5239587/93‑1 (MQ=255)

                                                                            gCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAcatc             >  2:5133066/1‑93 (MQ=255)

                                                                            gCCGTACGGTTGGGGGTATTGGCGTGATACCACGGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAcatc             <  2:6657139/93‑1 (MQ=255)

                                                                               gTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACatcatc          >  1:3557101/1‑93 (MQ=255)

                                                                               gTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGCGCCCAGTCAGCGTTTGGCGACatcatc          >  2:2546176/1‑93 (MQ=255)

                                                                                tGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATcc         <  1:7072569/93‑1 (MQ=255)

                                                                                tGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATcc         <  2:2102397/93‑1 (MQ=255)

                                                                                  cGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGt       >  2:4057695/1‑93 (MQ=255)

                                                                                  cGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGt       <  1:7401278/93‑1 (MQ=255)

                                                                                  cGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGt       <  2:6979119/93‑1 (MQ=255)

                                                                                   ggTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTg      <  2:6439117/93‑1 (MQ=255)

                                                                                    gTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGc     <  1:237069/93‑1 (MQ=255)

                                                                                     ttGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGcc    <  1:3940109/93‑1 (MQ=255)

                                                                                      tGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCg   >  1:5622775/1‑93 (MQ=255)

GTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGG  >  NC_000913/3160161‑3160340

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