BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,562,600	T→C	T274A (ACC→GCC)	yjiM ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,562,600	0	T	C	99.3%	456.7 / ‑5.5	146	T274A (ACC→GCC)	yjiM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (74/71); total (75/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 158 total aligned reads displayed.

GGAGCAGCCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGTGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTGA  >  NC_000913/4562511‑4562689

ggAGCAGCCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCTGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCa                                                                                        >  2:37912/1‑93 (MQ=255)

   gcagcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGtt                                                                                     >  2:1580680/1‑93 (MQ=255)

    cagcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGttt                                                                                    >  2:3041236/1‑93 (MQ=255)

    cagcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGttt                                                                                    >  2:284115/1‑93 (MQ=255)

     agcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGtttt                                                                                   >  2:3323462/1‑93 (MQ=255)

     agcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGtttt                                                                                   >  2:513061/1‑93 (MQ=255)

      gcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTc                                                                                  >  1:1417270/1‑93 (MQ=255)

      gcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTc                                                                                  >  2:1525179/1‑93 (MQ=255)

       ccAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCa                                                                                 <  2:1717804/93‑1 (MQ=255)

       ccAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCa                                                                                 >  1:3115734/1‑93 (MQ=255)

        cAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCAt                                                                                >  1:35742/1‑93 (MQ=255)

         aaTCGCCAGATTTTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCAtc                                                                               >  2:2309649/1‑92 (MQ=255)

          aTCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCACCGGCGCAGTTTTCATaa                                                                              >  2:121745/1‑93 (MQ=255)

           tCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAAc                                                                             <  1:3708417/93‑1 (MQ=255)

           tCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAAc                                                                             <  1:400470/93‑1 (MQ=255)

              ccAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGa                                                                          <  2:413735/93‑1 (MQ=255)

              ccAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGa                                                                          <  1:1563688/93‑1 (MQ=255)

              ccAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGa                                                                          <  2:3046774/93‑1 (MQ=255)

              ccAGATATTTATCCGCCAGCGCGACGTAGACATCGCCCGGTTCTGCCACGCATTGCTCTGTCGCTTTCGCCCCGGTGCAGTTTTCATAATCGa                                                                          >  1:4326356/1‑93 (MQ=255)

               cAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGAc                                                                         <  2:2784573/93‑1 (MQ=255)

                   tatTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAAcc                                                                     <  2:2937992/93‑1 (MQ=255)

                    atTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAAccc                                                                    <  1:3447333/93‑1 (MQ=255)

                    atTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCACGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAAccc                                                                    <  2:1897475/93‑1 (MQ=255)

                     tttATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCa                                                                   <  1:4633496/93‑1 (MQ=255)

                     tttATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCa                                                                   <  2:611273/93‑1 (MQ=255)

                       tATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAgc                                                                 >  2:1489372/1‑93 (MQ=255)

                       tATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAgc                                                                 <  2:3977918/93‑1 (MQ=255)

                       tATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAgc                                                                 >  1:828801/1‑93 (MQ=255)

                        agccGCCAGCGCGTCGTAGACATCGCTCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAgcc                                                                <  1:2481825/91‑1 (MQ=255)

                        aTCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAgcc                                                                <  1:3642730/93‑1 (MQ=255)

                        aTCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAgcc                                                                >  1:1463651/1‑93 (MQ=255)

                         tCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAgccg                                                               >  2:1639129/1‑93 (MQ=255)

                         tCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCCGCGCAGTGTTCATAACCGACAACCCAgccg                                                               >  2:3145603/1‑93 (MQ=255)

                          ccGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAgccgc                                                              >  1:4641174/1‑93 (MQ=255)

                            gCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCa                                                            >  1:725076/1‑93 (MQ=255)

                            gCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCa                                                            >  2:1064653/1‑93 (MQ=255)

                            gCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCa                                                            <  1:3972135/93‑1 (MQ=255)

                             ccAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCAt                                                           <  2:2743492/93‑1 (MQ=255)

                               aGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATtc                                                         >  2:2800174/1‑93 (MQ=255)

                                 cgcgTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACACCCCAGCCGCCATtctc                                                       <  1:312844/93‑1 (MQ=255)

                                  gcgTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTctct                                                      >  1:4529730/1‑93 (MQ=255)

                                     tcgtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCa                                                   >  1:329373/1‑93 (MQ=255)

                                       gtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAAt                                                 >  2:1605030/1‑93 (MQ=255)

                                        tAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATc                                                >  1:1025904/1‑93 (MQ=255)

                                          gACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATcgc                                              >  2:3446778/1‑93 (MQ=255)

                                           aCATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATcgcg                                             <  1:1036694/93‑1 (MQ=255)

                                           aCATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATcgcg                                             >  2:738507/1‑93 (MQ=255)

                                            cATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATcgcgc                                            >  2:2593924/1‑93 (MQ=255)

                                             aTCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATcgcgcg                                           >  1:4084625/1‑93 (MQ=255)

                                               cGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGca                                         >  1:2584512/1‑93 (MQ=255)

                                                gCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGCCAACCCAGCCGCCATTCTCTTCAATCGCGCGcac                                        >  1:501631/1‑93 (MQ=255)

                                                gCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGcac                                        >  1:921802/1‑93 (MQ=255)

                                                ccccGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGcac                                        <  2:2353603/92‑1 (MQ=255)

                                                 cccGTTTCCGCCACGCATTGCTCGGTCGCTTTCCCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGcacc                                       <  1:3604327/93‑1 (MQ=255)

                                                  ccGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGcacca                                      >  2:2728680/1‑93 (MQ=255)

                                                  ccGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGcacca                                      >  2:3261288/1‑93 (MQ=255)

                                                   cGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGcaccac                                     >  1:1501024/1‑93 (MQ=255)

                                                    gTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTAATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACt                                    <  2:2061129/93‑1 (MQ=255)

                                                     tttCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCTTAACCGCCAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACtt                                   >  2:3846926/1‑93 (MQ=255)

                                                     tttCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACtt                                   <  2:143779/93‑1 (MQ=255)

                                                     tttCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACtt                                   >  1:4612250/1‑93 (MQ=255)

                                                     tttCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTATCTTCAATCGCGCGcaccacgt                                   >  2:3752189/1‑91 (MQ=255)

                                                       tCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACtttt                                 >  2:122046/1‑93 (MQ=255)

                                                         tGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTc                               >  1:14309/1‑93 (MQ=255)

                                                         tGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTc                               <  1:2427148/93‑1 (MQ=255)

                                                         tGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTc                               <  1:4475874/93‑1 (MQ=255)

                                                          gCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTct                              <  2:476010/93‑1 (MQ=255)

                                                           ccACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTctg                             <  2:3644486/93‑1 (MQ=255)

                                                           ccACGCATTGCTCGGACGCTTTCGCCCCGGCGCAGTTTTCAGAACCGACAACCCAGCCGCCCTTCTCTTCAATCGCGCGCACCACTTTTTctg                             <  2:3644487/93‑1 (MQ=255)

                                                            cACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTctgc                            >  2:2842910/1‑93 (MQ=255)

                                                             aCGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTctgct                           <  2:1897009/93‑1 (MQ=255)

                                                             aCGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTctgct                           <  1:2029092/93‑1 (MQ=255)

                                                             aCGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTctgct                           >  2:1918182/1‑93 (MQ=255)

                                                             aCGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAAACGACAACTCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTctgct                           >  2:1448918/1‑93 (MQ=255)

                                                                 aTTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTgcgc                       >  2:3854187/1‑93 (MQ=255)

                                                                 aTTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTgcgc                       >  2:4314257/1‑93 (MQ=255)

                                                                 aTTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTgcgc                       >  2:2012264/1‑93 (MQ=255)

                                                                     cTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCgccgcc                   >  1:1122360/1‑93 (MQ=255)

                                                                      tCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCa                  <  1:2038442/93‑1 (MQ=255)

                                                                       cGGTCGCTTTCGCCCGGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCaa                 <  2:1056874/93‑1 (MQ=255)

                                                                       cGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCaa                 <  2:4673887/93‑1 (MQ=255)

                                                                         gTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATc               >  1:3336136/1‑93 (MQ=255)

                                                                          tCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCg              <  2:3276125/93‑1 (MQ=255)

                                                                           cGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCgg             <  2:3413555/93‑1 (MQ=255)

                                                                            gCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCggg            >  1:3018338/1‑93 (MQ=255)

                                                                            gCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCggg            >  2:1169287/1‑93 (MQ=255)

                                                                             cTTTCGCCCCGGCGCAGTTTTCATACCCGACACCCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGc           <  1:1862512/93‑1 (MQ=255)

                                                                             cTTTCGCCCCGGCGCAGTTTTCATACCCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGc           <  2:1793070/93‑1 (MQ=255)

                                                                             cTTTCGCCCCGGCGCAGTTTTCATAACCGACACCCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGc           <  1:2325857/93‑1 (MQ=255)

                                                                              tttCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCa          <  1:1713865/93‑1 (MQ=255)

                                                                              tttCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCa          <  1:261357/93‑1 (MQ=255)

                                                                              tttCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCa          >  2:3185795/1‑93 (MQ=255)

                                                                                tCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGc        <  1:3041236/93‑1 (MQ=255)

                                                                                tCGCCCCGGCGCAGTTTGCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGc        <  1:1879650/93‑1 (MQ=255)

                                                                                   ccccGGCGCAGTTTTCATAACCGACACCCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCgg     <  2:86425/93‑1 (MQ=255)

                                                                                   ccccGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCgg     <  2:1508364/93‑1 (MQ=255)

                                                                                     ccGGCGCAGTTTTCATAACCGACACCCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTg   <  1:4264135/93‑1 (MQ=255)

                                                                                      cGGCGCAGTTTTCATAACCGACACCCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTGa  <  2:1650348/93‑1 (MQ=255)

                                                                                      cGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTGa  >  1:2940218/1‑93 (MQ=255)

                                                                                      cGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTGa  >  2:1048934/1‑93 (MQ=255)

GGAGCAGCCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGTGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTGA  >  NC_000913/4562511‑4562689

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