BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,261,143	T→C	E223G (GAA→GGA)	qor ←	quinone oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,261,143	0	T	C	99.4%	1038.8 / ‑3.8	338	E223G (GAA→GGA)	qor	quinone oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (188/148); total (188/150)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

Only 100 of 362 total aligned reads displayed.

AGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCA  >  NC_000913/4261053‑4261232

agaATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTccc                                                                                         >  2:846291/1‑93 (MQ=255)

agaATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTccc                                                                                         <  2:5486900/93‑1 (MQ=255)

agaATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTccc                                                                                         <  1:7421351/93‑1 (MQ=255)

 gaATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTcccc                                                                                        <  2:3395673/93‑1 (MQ=255)

  aaTGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCa                                                                                       <  1:485298/93‑1 (MQ=255)

   aTGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAg                                                                                      >  2:6870/1‑93 (MQ=255)

    tGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAgg                                                                                     >  2:3322124/1‑93 (MQ=255)

    tGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAgg                                                                                     >  1:7849425/1‑93 (MQ=255)

     gCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGgt                                                                                    <  2:2700399/93‑1 (MQ=255)

     gCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGgt                                                                                    >  1:1001574/1‑93 (MQ=255)

      ccTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGgtg                                                                                   <  1:4431913/93‑1 (MQ=255)

       cTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGgtgt                                                                                  >  1:1371926/1‑93 (MQ=255)

        tAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGtc                                                                                 >  2:1410883/1‑93 (MQ=255)

           gTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTctct                                                                              >  2:665121/1‑93 (MQ=255)

           gTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATGAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTctct                                                                              >  1:4008427/1‑93 (MQ=255)

            ttCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTg                                                                             <  2:4631434/93‑1 (MQ=255)

            ttCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTg                                                                             >  2:5809999/1‑93 (MQ=255)

             tCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGc                                                                            >  2:7882246/1‑93 (MQ=255)

               acacCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGccc                                                                          >  1:4120064/1‑93 (MQ=255)

                  ccGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACg                                                                       <  1:2103400/93‑1 (MQ=255)

                    ggTAGCCGCACCGGATGAGTTGCCAAAACTGACCATGAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCGGCCCACGGa                                                                     <  2:3790644/93‑1 (MQ=255)

                    ggTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGa                                                                     >  2:3086767/1‑93 (MQ=255)

                    ggTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGa                                                                     >  2:1193556/1‑93 (MQ=255)

                     gTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGaa                                                                    >  1:7722172/1‑93 (MQ=255)

                      tAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAAt                                                                   >  2:6566651/1‑93 (MQ=255)

                        aCCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCg                                                                 <  1:2842804/93‑1 (MQ=255)

                         ccGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGt                                                                >  2:116844/1‑93 (MQ=255)

                          cGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTa                                                               <  2:6013918/93‑1 (MQ=255)

                          cGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTa                                                               <  2:4911837/93‑1 (MQ=255)

                          cGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTa                                                               <  1:6482880/93‑1 (MQ=255)

                           gCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTac                                                              >  1:1062331/1‑93 (MQ=255)

                            cACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTaca                                                             <  1:5235289/93‑1 (MQ=255)

                            cACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTaca                                                             >  1:1640126/1‑93 (MQ=255)

                            cACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTaca                                                             <  2:1958591/93‑1 (MQ=255)

                            cACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTaca                                                             >  2:161273/1‑93 (MQ=255)

                             aCCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTacac                                                            >  1:1499911/1‑93 (MQ=255)

                             aCCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTacac                                                            <  1:7658749/93‑1 (MQ=255)

                               cTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTAcacca                                                          >  1:3320212/1‑93 (MQ=255)

                                tgatgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTAcaccac                                                         <  2:3927871/93‑1 (MQ=255)

                                tgatgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTAcaccac                                                         <  1:130085/93‑1 (MQ=255)

                                tgatgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTAcaccac                                                         <  1:86748/93‑1 (MQ=255)

                                 gatgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCAcg                                                        >  1:3829876/1‑93 (MQ=255)

                                  atgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCAcgc                                                       >  2:1090341/1‑93 (MQ=255)

                                  atgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCAcgc                                                       <  1:2191007/93‑1 (MQ=255)

                                   tgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACAAcgcc                                                      >  1:1168053/1‑92 (MQ=255)

                                    gaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACgcgc                                                     >  2:6606963/1‑93 (MQ=255)

                                     aGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCa                                                    >  1:4552119/1‑93 (MQ=255)

                                        tGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACtt                                                 <  1:5055561/93‑1 (MQ=255)

                                         gCCAAAACTGACCATTAAGCCGCGGCTTTGCAGGCTACGCAGCGACCGTACCTATGCGTTTCTGCCCACCGCGTCGTACTCCCCGCGCTCtta                                                <  2:4827509/93‑2 (MQ=255)

                                         gCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACttt                                                >  2:5177511/1‑93 (MQ=255)

                                          ccAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTc                                               <  2:7417046/93‑1 (MQ=255)

                                           cAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCt                                              <  1:5475360/93‑1 (MQ=255)

                                           cAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCt                                              >  2:6414610/1‑93 (MQ=255)

                                            aaaaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCtt                                             >  2:7720411/1‑93 (MQ=255)

                                             aaaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTa                                            >  1:3138885/1‑93 (MQ=255)

                                               cccgACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTAcc                                          <  1:702428/90‑1 (MQ=255)

                                               aCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTAcc                                          <  1:7306155/93‑1 (MQ=255)

                                                 tGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTAccgc                                        >  1:5615093/1‑93 (MQ=255)

                                                   aCCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTAccgccg                                      >  2:7467663/1‑93 (MQ=255)

                                                     cATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCggc                                    >  2:2239211/1‑92 (MQ=255)

                                                        tAAGCCGCGGCGTTGCAGTCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGAt                                 <  2:5026011/93‑1 (MQ=255)

                                                        tAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGAt                                 >  2:6900837/1‑93 (MQ=255)

                                                        tAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGAt                                 <  1:5938555/93‑1 (MQ=255)

                                                           gCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGAtctc                              >  1:2741240/1‑93 (MQ=255)

                                                           gCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGAtctc                              >  1:3390606/1‑93 (MQ=255)

                                                             cgcgGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCtt                            >  2:7287101/1‑93 (MQ=255)

                                                             cgcgGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCtt                            >  2:5547890/1‑93 (MQ=255)

                                                              gcggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTTCCGCCGGCGATCTCttt                           <  2:2149361/93‑1 (MQ=255)

                                                               cggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCtttt                          <  1:247854/93‑1 (MQ=255)

                                                               cggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCtttt                          >  2:4775260/1‑93 (MQ=255)

                                                                ggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTa                         >  2:3641633/1‑93 (MQ=255)

                                                                 gcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTaa                        <  2:627099/93‑1 (MQ=255)

                                                                   gTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAAcc                      >  2:2554892/1‑93 (MQ=255)

                                                                   gTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAAcc                      <  1:3623597/93‑1 (MQ=255)

                                                                   gTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAAcc                      >  2:7092191/1‑93 (MQ=255)

                                                                   gTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCACACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAAcc                      <  1:2973129/93‑1 (MQ=255)

                                                                      gCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCt                   <  2:6996188/93‑1 (MQ=255)

                                                                      gCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCt                   >  1:6742077/1‑93 (MQ=255)

                                                                         ggCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGa                <  2:3038314/93‑1 (MQ=255)

                                                                         ggCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGa                >  2:2490217/1‑93 (MQ=255)

                                                                          gCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGAc               <  1:5382130/93‑1 (MQ=255)

                                                                            aaTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCa             <  1:6049228/93‑1 (MQ=255)

                                                                             aTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAg            >  2:4232570/1‑93 (MQ=255)

                                                                             aTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAg            >  2:5429620/1‑93 (MQ=255)

                                                                             aTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAg            >  2:4072160/1‑93 (MQ=255)

                                                                             aTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAg            >  2:3456325/1‑93 (MQ=255)

                                                                             aTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAg            >  2:4490964/1‑93 (MQ=255)

                                                                             aTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAg            >  1:3494039/1‑93 (MQ=255)

                                                                              tCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGa           <  1:4656588/93‑1 (MQ=255)

                                                                               ccAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGAt          <  1:305859/93‑1 (MQ=255)

                                                                               ccAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGAt          >  1:3566867/1‑93 (MQ=255)

                                                                               ccAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGAt          >  2:1371007/1‑93 (MQ=255)

                                                                               ccAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGAt          <  2:5882373/93‑1 (MQ=255)

                                                                                 aGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATcc        <  2:3268199/93‑1 (MQ=255)

                                                                                  gCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCct       <  1:2054985/93‑1 (MQ=255)

                                                                                  gCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCct       <  2:6163650/93‑1 (MQ=255)

                                                                                    gACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCAATCTCTTTTAACCGCTCGACCAGATCctct     <  1:2501969/93‑1 (MQ=255)

                                                                                     aCCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCctctc    >  1:617686/1‑92 (MQ=255)

                                                                                     aCCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCtt    >  2:2395595/1‑93 (MQ=255)

                                                                                       cGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCa  <  2:1816423/93‑1 (MQ=255)

AGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCA  >  NC_000913/4261053‑4261232

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