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	98.6%	457.4 / ‑2.5	146	E223G (GAA→GGA)	qor	quinone oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (86/58); total (88/58)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-01

Only 100 of 160 total aligned reads displayed.

AGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTC  >  NC_000913/4261053‑4261231

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

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

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

 gaATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTcccc                                                                                       >  1:3340514/1‑93 (MQ=255)

 gaATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTcccc                                                                                       >  1:3807129/1‑93 (MQ=255)

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

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

    tGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGAAGGCAATCCAGCGACCGTCCCCAgg                                                                                    >  1:4045423/1‑93 (MQ=255)

      ccTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGgtg                                                                                  >  1:1543763/1‑93 (MQ=255)

      ccTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGgtg                                                                                  >  1:2182957/1‑93 (MQ=255)

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

        tAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGtc                                                                                <  1:4414774/93‑1 (MQ=255)

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

         aaGTTCACACCGGTAACCGCACCTGATGAGTTGCTAAAACTGACCATTACGCAGCGGCGTTGCGGGCGGTCCAGCGACCGTCCCCAGGTGtct                                                                               >  2:3788109/1‑93 (MQ=255)

         aaGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGtct                                                                               <  2:1884850/93‑1 (MQ=255)

          aGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGtctc                                                                              >  1:1993126/1‑93 (MQ=255)

          aGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGtctc                                                                              <  1:652116/93‑1 (MQ=255)

           gTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTctct                                                                             <  1:367474/93‑1 (MQ=255)

           gTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGAGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTctct                                                                             >  2:4405201/1‑93 (MQ=255)

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

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

              cacaCCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGcc                                                                          >  1:4157454/1‑93 (MQ=255)

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

                 acCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCAc                                                                       <  1:985527/93‑1 (MQ=255)

                 acCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCAc                                                                       >  2:3924141/1‑93 (MQ=255)

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

                  ccGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACg                                                                      >  2:3023604/1‑93 (MQ=255)

                  ccGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACg                                                                      >  1:1289816/1‑93 (MQ=255)

                   cGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACgg                                                                     >  2:1816906/1‑93 (MQ=255)

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

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

                      tAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAAt                                                                  <  2:4487278/93‑1 (MQ=255)

                      tAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAAt                                                                  <  2:2889691/93‑1 (MQ=255)

                       aaCCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATc                                                                 >  2:4157962/1‑93 (MQ=255)

                        aCCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCg                                                                >  2:1202747/1‑93 (MQ=255)

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

                          cGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCgca                                                              >  1:344995/1‑91 (MQ=255)

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

                              ccTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTAcacc                                                          >  2:1408003/1‑93 (MQ=255)

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

                               cTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTAcacca                                                         <  2:2285542/93‑1 (MQ=255)

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

                                   tgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCAAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCAcgcg                                                     >  2:2890782/1‑93 (MQ=255)

                                      gTTGCCAAAACTGATCATTACGCCGCGGCGTTGCAGGCAATACAGCGACCGTCCACAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCAc                                                  >  2:4636974/1‑93 (MQ=255)

                                      gTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCATGCGCAc                                                  >  2:2153388/1‑93 (MQ=255)

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

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

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

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

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

                                              aaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTAc                                          <  1:2231419/93‑1 (MQ=255)

                                               aCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTAcc                                         >  1:1065286/1‑93 (MQ=255)

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

                                                cTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTAccg                                        >  1:2614732/1‑93 (MQ=255)

                                                     cATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGt                                   >  2:2783659/1‑93 (MQ=255)

                                                     cATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGt                                   <  2:1730508/93‑1 (MQ=255)

                                                      aTTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTg                                  >  1:3750051/1‑93 (MQ=255)

                                                       ttAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGa                                 >  1:2962845/1‑93 (MQ=255)

                                                       ttAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGa                                 <  2:3576097/93‑1 (MQ=255)

                                                       agaaGCCGTGGCGTTGCGGGCAATCCAGCGACCGTCCCCAGGTGTCTCTCCCCACGAAAACGTACACCACGCGCACTTTCTTTCCGCCGGTGa                                 <  2:3576098/91‑1 (MQ=255)

                                                        tAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGAt                                >  2:3890146/1‑93 (MQ=255)

                                                         acgCCGCGGCGTTGCACGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGAtc                               <  2:3872375/91‑1 (MQ=255)

                                                         aaGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGAtc                               >  1:1998315/1‑93 (MQ=255)

                                                          aGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGAtct                              <  2:4679809/93‑1 (MQ=255)

                                                             cgcgGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCtt                           >  1:4048925/1‑93 (MQ=255)

                                                             cgcgGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCtt                           >  2:641681/1‑93 (MQ=255)

                                                             cgcgGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCtt                           >  2:2028840/1‑93 (MQ=255)

                                                                 gcgTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTaa                       >  2:4520186/1‑93 (MQ=255)

                                                                 gcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTac                       >  1:3661233/1‑92 (MQ=255)

                                                                 gcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTaa                       >  2:1761766/1‑93 (MQ=255)

                                                                 gcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTaa                       >  1:4337130/1‑93 (MQ=255)

                                                                 gcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTaa                       <  1:4599390/93‑1 (MQ=255)

                                                                  cgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAAc                      <  2:1327894/93‑1 (MQ=255)

                                                                   gTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAAcc                     >  1:1085796/1‑93 (MQ=255)

                                                                    ttGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCg                    >  2:3637838/1‑93 (MQ=255)

                                                                     tGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGc                   <  1:1726355/93‑1 (MQ=255)

                                                                       cAGGCAATCCAGCTACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTc                 <  2:736111/93‑1 (MQ=255)

                                                                       cAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTc                 <  1:611465/93‑1 (MQ=255)

                                                                        aGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCg                >  1:415457/1‑93 (MQ=255)

                                                                        aGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCg                >  2:3861640/1‑93 (MQ=255)

                                                                            aaTCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCa            >  1:1577191/1‑93 (MQ=255)

                                                                            aaTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCa            <  1:3415305/93‑1 (MQ=255)

                                                                            aaTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCa            <  1:3283425/93‑1 (MQ=255)

                                                                             aTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAg           >  2:1325930/1‑93 (MQ=255)

                                                                             aTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAg           <  1:1914770/93‑1 (MQ=255)

                                                                              tCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGa          >  2:538650/1‑93 (MQ=255)

                                                                              tCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGa          <  2:2814785/93‑1 (MQ=255)

                                                                              tCCAGCGACCGTCCCCAGGAGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTTACCACTCGACCAGa          >  1:1695203/1‑93 (MQ=255)

                                                                               ccAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGAt         <  1:704163/93‑1 (MQ=255)

                                                                                cAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATc        >  1:3077458/1‑93 (MQ=255)

                                                                                 aGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATcc       >  1:2340626/1‑93 (MQ=255)

                                                                                  gCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCct      >  2:127713/1‑93 (MQ=255)

                                                                                  gCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCct      >  2:3323780/1‑93 (MQ=255)

                                                                                  gCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCct      <  1:758519/93‑1 (MQ=255)

                                                                                   cGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCctc     >  1:2939338/1‑93 (MQ=255)

                                                                                    gACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCctct    >  2:3695249/1‑93 (MQ=255)

                                                                                    gACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCctct    <  2:1470116/93‑1 (MQ=255)

                                                                                     aCCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATACTCtt   >  2:1055851/1‑93 (MQ=255)

                                                                                      ccGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTc  <  1:2317714/93‑1 (MQ=255)

                                                                                      ccGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTc  >  1:1323439/1‑93 (MQ=255)

AGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTC  >  NC_000913/4261053‑4261231

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