BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,112,040	T→C	A110A (GCA→GCG)	b2972 ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,112,040	0	T	C	99.7%	1095.7 / ‑5.3	352	A110A (GCA→GCG)	b2972	putative peptidase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (182/168); total (183/169)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

Only 100 of 369 total aligned reads displayed.

ACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAG  >  NC_000913/3111950‑3112129

 cTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCt                                                                                        >  2:5699158/1‑93 (MQ=255)

 cTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCt                                                                                        <  2:4727343/93‑1 (MQ=255)

    gcgGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTc                                                                                     <  1:536516/93‑1 (MQ=255)

     cgGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCa                                                                                    >  2:2003103/1‑93 (MQ=255)

      ggCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCaa                                                                                   >  1:4740147/1‑93 (MQ=255)

         aaTCAGCCAGGCGGATAATATCCTCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATaa                                                                                <  2:4592079/93‑1 (MQ=255)

         aaTCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATaa                                                                                >  1:926450/1‑93 (MQ=255)

         aaTCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATaa                                                                                <  2:4627850/93‑1 (MQ=255)

          aTCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAAc                                                                               >  2:121788/1‑93 (MQ=255)

            cAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTc                                                                             <  2:704675/93‑1 (MQ=255)

                cAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTccacc                                                                         <  2:3976160/93‑1 (MQ=255)

                cAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTccacc                                                                         <  2:1520168/93‑1 (MQ=255)

                cAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTcaacc                                                                         >  1:5665790/1‑93 (MQ=255)

                 aGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTccacca                                                                        <  1:5990456/93‑1 (MQ=255)

                 aGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTccacca                                                                        <  2:7223859/93‑1 (MQ=255)

                 aGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTccacca                                                                        >  1:5136028/1‑93 (MQ=255)

                 aGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTccacca                                                                        <  1:1395834/93‑1 (MQ=255)

                   gCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGc                                                                      <  1:2595927/93‑1 (MQ=255)

                      gATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGa                                                                   >  2:6588228/1‑93 (MQ=255)

                      gATAATACCATCGCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACGCCACCAGCGGa                                                                   <  1:7438778/93‑1 (MQ=255)

                       ataataTCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGAt                                                                  <  2:6471004/93‑1 (MQ=255)

                         aataTCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATaa                                                                <  2:16313/93‑1 (MQ=255)

                           tatCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACg                                                              >  2:367157/1‑93 (MQ=255)

                             tcatcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCt                                                            >  2:2496613/1‑93 (MQ=255)

                             tcatcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCt                                                            >  1:1147796/1‑93 (MQ=255)

                              catcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCtt                                                           <  2:1732182/93‑1 (MQ=255)

                               atcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTg                                                          <  1:2964585/93‑1 (MQ=255)

                                ccaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGc                                                         <  1:4591094/92‑1 (MQ=255)

                                  aCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTg                                                       >  1:1977023/1‑93 (MQ=255)

                                   ccGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGa                                                      >  1:3697150/1‑93 (MQ=255)

                                    cGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAt                                                     <  1:1475374/93‑1 (MQ=255)

                                    cGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAt                                                     <  2:7359634/93‑1 (MQ=255)

                                      ccAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAttt                                                   <  1:5244224/93‑1 (MQ=255)

                                      ccACCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAttt                                                   <  2:6957081/93‑1 (MQ=255)

                                       cAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAtttt                                                  >  1:3379730/1‑93 (MQ=255)

                                       cAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAtttt                                                  >  1:1021165/1‑93 (MQ=255)

                                       cAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAtttt                                                  >  1:6965015/1‑93 (MQ=255)

                                       cAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAtttt                                                  >  2:4173124/1‑93 (MQ=255)

                                       cAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAtttt                                                  <  2:2773167/93‑1 (MQ=255)

                                        aGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTc                                                 <  2:159354/93‑1 (MQ=255)

                                          cccccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGc                                               <  1:5729539/93‑1 (MQ=255)

                                           ccccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGcc                                              <  1:1256907/93‑1 (MQ=255)

                                           ccccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGcc                                              <  1:4209314/93‑1 (MQ=255)

                                            cccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCt                                             >  1:7117766/1‑93 (MQ=255)

                                            cccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCt                                             >  2:6254674/1‑93 (MQ=255)

                                              cATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGa                                           >  2:5858490/1‑93 (MQ=255)

                                              cATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGa                                           <  1:6254675/93‑1 (MQ=255)

                                               aTCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGAc                                          >  1:2440433/1‑93 (MQ=255)

                                               aTCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGAc                                          >  2:3079486/1‑93 (MQ=255)

                                                tCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACa                                         <  1:718832/93‑1 (MQ=255)

                                                tCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACa                                         >  2:1163854/1‑93 (MQ=255)

                                                tCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACa                                         <  1:3501114/93‑1 (MQ=255)

                                                tCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACa                                         >  1:5393742/1‑93 (MQ=255)

                                                tCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACa                                         >  1:6507261/1‑93 (MQ=255)

                                                tCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTACACCAGCGGATAACGCTTGCTGATTTTCGCCTGACa                                         >  1:4467152/1‑93 (MQ=255)

                                                  cACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGt                                       <  2:910319/93‑1 (MQ=255)

                                                  cACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGt                                       >  2:5629958/1‑93 (MQ=255)

                                                  cACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGt                                       >  1:2458944/1‑93 (MQ=255)

                                                     tttCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTcgc                                    <  2:2861131/93‑1 (MQ=255)

                                                      ttCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTcgcg                                   <  2:4197755/93‑1 (MQ=255)

                                                        cTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGc                                 <  1:1596557/93‑1 (MQ=255)

                                                           gCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAGCTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAgc                              <  2:7445043/93‑1 (MQ=255)

                                                           gCCAGACAAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAgc                              <  1:7290523/93‑1 (MQ=255)

                                                             cagaccagacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAgcgc                            >  2:6045577/1‑93 (MQ=255)

                                                               gaccagacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGccc                          >  2:2707607/1‑93 (MQ=255)

                                                               gaccagacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGccc                          <  2:5394529/93‑1 (MQ=255)

                                                               gaccagacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGccc                          >  2:4244750/1‑93 (MQ=255)

                                                                accagacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCt                         >  1:803759/1‑93 (MQ=255)

                                                                 ccagacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCtt                        <  1:3620809/93‑1 (MQ=255)

                                                                 ccagacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCtt                        <  2:7714083/93‑1 (MQ=255)

                                                                  cagacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCttt                       >  2:5080613/1‑93 (MQ=255)

                                                                   agacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTg                      <  2:7515384/93‑1 (MQ=255)

                                                                   agacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTg                      >  2:4781539/1‑93 (MQ=255)

                                                                    gacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGa                     <  1:3192920/93‑1 (MQ=255)

                                                                    gacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGa                     >  2:5181026/1‑93 (MQ=255)

                                                                      cTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGc                   <  1:6580777/93‑1 (MQ=255)

                                                                       tCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCa                  >  1:612455/1‑93 (MQ=255)

                                                                       ccGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCCCCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCa                  <  2:2245820/92‑1 (MQ=255)

                                                                         gCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATc                <  1:3970098/93‑1 (MQ=255)

                                                                         gCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATc                <  1:4867569/93‑1 (MQ=255)

                                                                         gCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATc                >  2:4449102/1‑93 (MQ=255)

                                                                         gCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATc                <  2:5898866/93‑1 (MQ=255)

                                                                         gCCAGCAAAAAAACAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATc                <  2:3248816/93‑1 (MQ=255)

                                                                          ccAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCa               <  2:1356080/93‑1 (MQ=255)

                                                                            aGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGc             <  2:4098703/93‑1 (MQ=255)

                                                                            aGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGc             <  1:1799357/93‑1 (MQ=255)

                                                                               aaaaaaTCAAGCGCTGTCACTACCTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCTGCAGCGCCCTTTGAGCATCAGCCAg          <  2:711320/93‑1 (MQ=255)

                                                                               aaaaaaGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAg          >  2:6736915/1‑93 (MQ=255)

                                                                                aaaaaGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAgg         >  2:2918491/1‑93 (MQ=255)

                                                                                aaaaaGCAAGCGCTGTCAATAACGCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAgg         <  1:1068756/93‑1 (MQ=255)

                                                                                 aaaaGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGga        >  2:840767/1‑93 (MQ=255)

                                                                                 aaaaGCAAGCGCTGTCAATAACGCCCCCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGga        <  2:3598846/93‑1 (MQ=255)

                                                                                   aagcaagcGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGgaga      <  1:7591706/93‑1 (MQ=255)

                                                                                    aGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGaa     >  2:5508603/1‑93 (MQ=255)

                                                                                    aGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGaa     >  1:4313256/1‑93 (MQ=255)

                                                                                    aGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGaa     <  1:40772/93‑1 (MQ=255)

                                                                                     gCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAAc    <  1:6804276/93‑1 (MQ=255)

                                                                                     gCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAAc    >  1:2766441/1‑93 (MQ=255)

                                                                                       aaGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAg  <  1:7815231/93‑1 (MQ=255)

                                                                                       aaGCGCGGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAg  <  2:2630391/93‑1 (MQ=255)

ACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAG  >  NC_000913/3111950‑3112129

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