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.7%	1230.4 / ‑5.0	387	E223G (GAA→GGA)	qor	quinone oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (180/206); total (181/206)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 413 total aligned reads displayed.

AGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCA  >  NC_000913/4261053‑4261232

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

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

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

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

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

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

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

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

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

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

              cacaCCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGcc                                                                           <  2:2574760/93‑1 (MQ=255)

                cacCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCa                                                                         <  2:7566017/93‑1 (MQ=255)

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

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

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

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

                        aCCGCACCTCATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCg                                                                 <  1:8527093/93‑1 (MQ=255)

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

                          cGCACCTGATGAGTCGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTa                                                               <  1:5690046/93‑1 (MQ=255)

                           gCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTac                                                              <  1:7545097/93‑1 (MQ=255)

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

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

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

                              ccTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTAcacc                                                           <  2:4821441/93‑1 (MQ=255)

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

                              ccTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGGCTCTGCCCACGGAATCGTAcacc                                                           >  2:5563061/1‑93 (MQ=255)

                               cTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGCCCGTCCCCAGGTGTCTCTGCCCACGGAATCGTAcacca                                                          <  2:3993289/93‑1 (MQ=255)

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

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

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

                                tgatgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTAcaccac                                                         >  2:7819633/1‑93 (MQ=255)

                                tgatgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGGCTCTGCCCACGGAATCGTAcaccac                                                         >  1:9448585/1‑93 (MQ=255)

                                tgatgaGTTGCAAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTAcaccac                                                         <  1:7183929/93‑1 (MQ=255)

                                   tgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCAcgcg                                                      >  2:6673127/1‑93 (MQ=255)

                                   tgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCAcgcg                                                      >  1:4770821/1‑93 (MQ=255)

                                   tgaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCAcgcg                                                      >  1:8436632/1‑93 (MQ=255)

                                    gaGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACgcgc                                                     <  1:9756775/93‑1 (MQ=255)

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

                                      gTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCAc                                                   <  1:688282/93‑1 (MQ=255)

                                      gTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCAc                                                   <  1:2516417/93‑1 (MQ=255)

                                           cAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGTATCGTACACCACGCGCACTTTCt                                              >  2:7202557/1‑93 (MQ=255)

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

                                             aaaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTa                                            <  1:9166519/93‑1 (MQ=255)

                                              aaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTAc                                           >  1:5544618/1‑93 (MQ=255)

                                              aaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTAc                                           >  1:7315390/1‑93 (MQ=255)

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

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

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

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

                                                cTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTAccg                                         <  2:375531/93‑1 (MQ=255)

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

                                                 tGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTAccgc                                        <  2:2054176/93‑1 (MQ=255)

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

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

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

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

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

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

                                                         aaGCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGAtc                                <  1:602931/93‑1 (MQ=255)

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

                                                           gCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGAtctc                              >  2:1648273/1‑93 (MQ=255)

                                                           gCCGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGAtctc                              >  1:4569702/1‑93 (MQ=255)

                                                            ccGCGGCGTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATctct                             <  2:9002849/93‑1 (MQ=255)

                                                            cagcggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATctct                             <  1:6791573/91‑1 (MQ=255)

                                                              gcggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCttt                           >  1:6523671/1‑93 (MQ=255)

                                                              gcggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCttt                           <  2:8025300/93‑1 (MQ=255)

                                                              gcggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCttt                           >  1:8201119/1‑93 (MQ=255)

                                                                ggcgTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTa                         <  1:4143812/93‑1 (MQ=255)

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

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

                                                                   gTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAAcc                      <  1:1294569/93‑1 (MQ=255)

                                                                   gTTGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAAcc                      <  2:4410703/93‑1 (MQ=255)

                                                                   gTTGCAGGCAATCCAGCGACCGGCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAAcc                      >  1:5849129/1‑93 (MQ=255)

                                                                     tGCAGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGc                    >  1:4739343/1‑93 (MQ=255)

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

                                                                        aGGCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACAGCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCg                 >  1:6369251/1‑93 (MQ=255)

                                                                         ggCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCgc                >  2:2244484/1‑92 (MQ=255)

                                                                         ggCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGa                <  2:9266216/93‑1 (MQ=255)

                                                                         ggCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGa                >  1:9578793/1‑93 (MQ=255)

                                                                         ggCAATCCAGCGACCGTCCCCAGGAGTCTCTGCCCACGGAATCGTACACCACGCGAACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCgc                >  1:2087777/1‑92 (MQ=255)

                                                                         ggCAATCAAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGa                >  2:4657677/1‑93 (MQ=255)

                                                                          gCAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGAc               >  2:9402041/1‑93 (MQ=255)

                                                                           cAATCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGAcc              <  2:7968647/93‑1 (MQ=255)

                                                                            aaTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCa             >  1:9834469/1‑93 (MQ=255)

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

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

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

                                                                             aTCCAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTCACCGCCGGTGATCTCTTTTAACCGCTCGACCAg            >  1:2604201/1‑93 (MQ=255)

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

                                                                                cAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGAta         >  2:5837722/1‑92 (MQ=255)

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

                                                                                cAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATc         <  2:8854610/93‑1 (MQ=255)

                                                                                cAGCGACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATc         <  1:7381422/93‑1 (MQ=255)

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

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

                                                                                    gACCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGCTCctct     >  1:322356/1‑93 (MQ=255)

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

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

                                                                                     aCCGTCCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCtt    >  2:8593930/1‑93 (MQ=255)

                                                                                      ccGTCCCCAGGTGCCTCTGCCCACGGAATTGTACACCACGCGCACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGAGCctctcc   >  2:8238842/1‑91 (MQ=255)

AGAATGCCTAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCA  >  NC_000913/4261053‑4261232

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