BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,986	T→C	T199A (ACG→GCG)	yidA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,986	0	T	C	82.4%	801.8 / 168.6	375	T199A (ACG→GCG)	yidA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (31/35); new base C (159/150); total (190/185)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

Only 100 of 377 total aligned reads displayed.

TTCAATCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCG  >  NC_000913/3873894‑3874078

  cAATCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCAc                                                                                            >  2:4246156/1‑93 (MQ=255)

   aaTCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCAcc                                                                                           <  1:939230/93‑1 (MQ=255)

   aaTCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCAcc                                                                                           >  2:6004955/1‑93 (MQ=255)

   aaTCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCAcc                                                                                           >  1:4595602/1‑93 (MQ=255)

    aTCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCt                                                                                          >  1:1183612/1‑93 (MQ=255)

     tCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCAGCCAGTGATTTCACCCCCGCACCtt                                                                                         >  2:5824861/1‑93 (MQ=255)

      cATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCttt                                                                                        >  1:4973086/1‑93 (MQ=255)

        ttGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGt                                                                                      >  1:357689/1‑93 (MQ=255)

        ttGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGt                                                                                      <  1:3785843/93‑1 (MQ=255)

           nnatatCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAACACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTaa                                                                                   <  1:1407051/91‑1 (MQ=255)

              tatCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAAcnn                                                                                >  2:2764542/1‑91 (MQ=255)

              tatCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAAcgc                                                                                >  1:3633642/1‑93 (MQ=255)

              tatCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAAcgc                                                                                >  2:1613983/1‑93 (MQ=255)

                  gTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGttt                                                                            >  1:3371774/1‑93 (MQ=255)

                   ttttCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTa                                                                           >  1:5879247/1‑93 (MQ=255)

                   ttttCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTa                                                                           >  2:6487343/1‑93 (MQ=255)

                    tttCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTAt                                                                          <  1:7432400/93‑1 (MQ=255)

                      tCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCg                                                                        >  1:571414/1‑93 (MQ=255)

                       ccTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGa                                                                       <  2:1847107/93‑1 (MQ=255)

                       ccTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGa                                                                       <  1:3279686/93‑1 (MQ=255)

                        cTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAg                                                                      <  2:734117/93‑1 (MQ=255)

                         tGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAgg                                                                     >  2:2706465/1‑93 (MQ=255)

                         tGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAgg                                                                     >  2:3054313/1‑93 (MQ=255)

                         tGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAgg                                                                     >  1:7344292/1‑93 (MQ=255)

                         tGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAgg                                                                     >  2:6267854/1‑93 (MQ=255)

                         tGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAgg                                                                     >  2:7437080/1‑93 (MQ=255)

                          gATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGa                                                                    <  2:2045284/93‑1 (MQ=255)

                          gATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGa                                                                    <  2:3920217/93‑1 (MQ=255)

                          gATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGa                                                                    >  1:1755541/1‑93 (MQ=255)

                          gATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGa                                                                    >  1:6074513/1‑93 (MQ=255)

                            tcgccaatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGAtt                                                                  <  2:3463345/93‑1 (MQ=255)

                            tcgccaatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGAtt                                                                  <  2:5750309/93‑1 (MQ=255)

                            tcgccaatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGAtt                                                                  >  1:6183035/1‑93 (MQ=255)

                            tcgccaatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGAtt                                                                  <  1:1584955/93‑1 (MQ=255)

                              gccaatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTc                                                                >  2:5441825/1‑93 (MQ=255)

                              gccaatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTc                                                                <  1:6914088/93‑1 (MQ=255)

                                  atcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAgg                                                            >  1:1096159/1‑93 (MQ=255)

                                   tcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGa                                                           <  2:761227/93‑1 (MQ=255)

                                    cgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGaa                                                          <  1:4416589/93‑1 (MQ=255)

                                     gccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAg                                                         >  2:6910858/1‑93 (MQ=255)

                                     gccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAg                                                         >  2:4521166/1‑93 (MQ=255)

                                      ccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGt                                                        >  1:2738595/1‑93 (MQ=255)

                                      ccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGt                                                        >  2:7262409/1‑93 (MQ=255)

                                      ccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGt                                                        >  1:4111606/1‑93 (MQ=255)

                                      ccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGt                                                        <  2:1097130/93‑1 (MQ=255)

                                       caTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTa                                                       >  1:5317488/1‑93 (MQ=255)

                                       caTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTa                                                       <  1:5955124/93‑1 (MQ=255)

                                       caTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTa                                                       <  1:236539/93‑1 (MQ=255)

                                       caTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTa                                                       >  2:4361610/1‑93 (MQ=255)

                                        aTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTAc                                                      <  2:4694232/93‑1 (MQ=255)

                                         tGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACg                                                     >  2:80694/1‑93 (MQ=255)

                                          gATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACgg                                                    >  2:507998/1‑93 (MQ=255)

                                           aTTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGgc                                                   >  1:2040591/1‑93 (MQ=255)

                                            tttCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGgcg                                                  >  1:7142394/1‑93 (MQ=255)

                                            tttCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGgcg                                                  <  1:3156303/93‑1 (MQ=255)

                                             ttcttcCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGgcgc                                                 >  1:57534/1‑93 (MQ=255)

                                               cttcCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCAc                                               <  2:1334791/93‑1 (MQ=255)

                                                 tcCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACtt                                             >  1:5096271/1‑93 (MQ=255)

                                                    ggTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTc                                          <  1:3412457/93‑1 (MQ=255)

                                                    ggTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTc                                          <  2:2668262/93‑1 (MQ=255)

                                                     gTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCa                                         <  2:6735290/93‑1 (MQ=255)

                                                     gTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCa                                         >  2:518926/1‑93 (MQ=255)

                                                      tttAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCAAGGAAGTACGGCGCACTTTTCAg                                        >  1:6661060/1‑93 (MQ=255)

                                                       ttAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGc                                       >  2:910325/1‑93 (MQ=255)

                                                       ttAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGc                                       >  1:5624631/1‑93 (MQ=255)

                                                         aaTACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGGTTTCGAGGAAGTACGGCGCACTTTTCAGCAc                                     <  1:4489229/93‑1 (MQ=255)

                                                          aTACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACg                                    >  2:5670507/1‑93 (MQ=255)

                                                          aTACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACg                                    <  2:3270714/93‑1 (MQ=255)

                                                           tACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACgg                                   >  1:1788164/1‑93 (MQ=255)

                                                           tACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTCATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACgg                                   <  1:5249426/93‑1 (MQ=255)

                                                            aCCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGt                                  >  2:2183938/1‑93 (MQ=255)

                                                            aCCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGt                                  <  1:7030115/93‑1 (MQ=255)

                                                             ccTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGta                                 >  1:1270634/1‑93 (MQ=255)

                                                             ccTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGta                                 >  2:1012266/1‑93 (MQ=255)

                                                             ccTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGta                                 <  1:3695054/93‑1 (MQ=255)

                                                                 acacGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATAtt                             >  2:3773989/1‑93 (MQ=255)

                                                                 acacGTCGGCCAGTGATTTAACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATAtt                             >  1:6229072/1‑93 (MQ=255)

                                                                   acGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTtc                           >  2:2384695/1‑93 (MQ=255)

                                                                     gTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTtctc                         <  2:6602057/93‑1 (MQ=255)

                                                                      tCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTctct                        >  2:6192444/1‑93 (MQ=255)

                                                                      tCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTctct                        <  1:2942451/93‑1 (MQ=255)

                                                                      tCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTctct                        <  2:4924543/93‑1 (MQ=255)

                                                                       cGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCtt                       <  1:5388137/93‑1 (MQ=255)

                                                                        ggCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCttt                      >  1:1352919/1‑93 (MQ=255)

                                                                          ccAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCa                    <  2:3129544/93‑1 (MQ=255)

                                                                           cAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCAc                   >  1:423563/1‑93 (MQ=255)

                                                                            aGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACt                  >  2:7126524/1‑93 (MQ=255)

                                                                             gTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACtt                 <  1:4688200/93‑1 (MQ=255)

                                                                             gTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACtt                 >  2:4046875/1‑93 (MQ=255)

                                                                             gTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACtt                 >  1:2671159/1‑93 (MQ=255)

                                                                             gTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACtt                 >  1:2196887/1‑93 (MQ=255)

                                                                               gATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTcc               <  2:1955523/93‑1 (MQ=255)

                                                                                 tttCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTg             >  1:6963655/1‑93 (MQ=255)

                                                                                   tCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCg           >  1:1674938/1‑93 (MQ=255)

                                                                                      cccccGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGaa        <  1:1012409/93‑1 (MQ=255)

                                                                                      cccccGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGaa        >  2:6786665/1‑93 (MQ=255)

                                                                                       ccccGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAAt       >  2:1516930/1‑93 (MQ=255)

                                                                                       ccccGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAAt       >  1:2605362/1‑93 (MQ=255)

                                                                                        cccGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATa      <  2:1679517/93‑1 (MQ=255)

                                                                                         ccGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATAc     <  1:2268460/93‑1 (MQ=255)

TTCAATCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCG  >  NC_000913/3873894‑3874078

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