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	98.4%	572.9 / ‑1.2	183	A110A (GCA→GCG)	b2972	putative peptidase
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (114/66); total (116/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

Only 100 of 189 total aligned reads displayed.

ACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAG  >  NC_000913/3111950‑3112129

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

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

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

      gccAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCaa                                                                                   <  2:563803/91‑1 (MQ=255)

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

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

           tCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACt                                                                              >  1:1393186/1‑93 (MQ=255)

           tCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACt                                                                              <  2:3270997/93‑1 (MQ=255)

             aGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTcc                                                                            >  2:366945/1‑93 (MQ=255)

              gCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTcca                                                                           >  2:1175360/1‑93 (MQ=255)

               ccAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTccac                                                                          <  2:1266774/93‑1 (MQ=255)

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

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

                 aGGCGGATAATATCACCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTccacca                                                                        <  2:2386786/93‑1 (MQ=255)

                  ggCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAg                                                                       >  1:2990307/1‑93 (MQ=255)

                  ggCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAg                                                                       >  2:4452671/1‑93 (MQ=255)

                     ggATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCgg                                                                    >  1:3433361/1‑93 (MQ=255)

                      gATAATATCATCCCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGa                                                                   <  2:934160/93‑1 (MQ=255)

                      gATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGa                                                                   <  2:3156431/93‑1 (MQ=255)

                      gATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGa                                                                   <  1:739731/93‑1 (MQ=255)

                        taataTCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATa                                                                 >  1:2936217/1‑93 (MQ=255)

                         aataTCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATaa                                                                >  1:1565690/1‑93 (MQ=255)

                         aataTCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATaa                                                                >  1:4168932/1‑93 (MQ=255)

                             tcatcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCt                                                            <  2:2726531/93‑1 (MQ=255)

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

                             tcatcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCt                                                            <  1:1286516/93‑1 (MQ=255)

                               atcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTg                                                          >  1:567586/1‑93 (MQ=255)

                                tcaCCGCCAGCCCCCATCCACTTTCTGGCCTGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGc                                                         >  1:2389696/1‑93 (MQ=255)

                                 caCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCt                                                        <  2:3314071/93‑1 (MQ=255)

                                   ccGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGa                                                      <  2:3853106/93‑1 (MQ=255)

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

                                     gCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAtt                                                    >  1:493937/1‑93 (MQ=255)

                                      ccAGCCCCCATCCCCTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAttt                                                   <  2:972082/93‑1 (MQ=255)

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

                                         gCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCg                                                >  2:1943591/1‑93 (MQ=255)

                                         gCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCg                                                >  1:967066/1‑93 (MQ=255)

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

                                           ccccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGcc                                              >  1:2267924/1‑93 (MQ=255)

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

                                            cccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCt                                             <  2:1572228/93‑1 (MQ=255)

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

                                             ccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTg                                            <  2:2495961/93‑1 (MQ=255)

                                             ccATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTg                                            >  2:2528972/1‑93 (MQ=255)

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

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

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

                                               aTCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTgga                                          >  2:1177845/1‑91 (MQ=255)

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

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

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

                                                 ccACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAg                                        <  1:2864141/93‑1 (MQ=255)

                                                 ccACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAg                                        >  1:2015024/1‑93 (MQ=255)

                                                 ccACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAg                                        <  1:1226368/93‑1 (MQ=255)

                                                   ccTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGTATAACGCTTGCTGATTTTCGCCTGACAGTc                                      <  1:458462/92‑1 (MQ=255)

                                                   aCTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTc                                      >  2:3399939/1‑93 (MQ=255)

                                                   aCTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTc                                      >  2:2273466/1‑93 (MQ=255)

                                                    cTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTcg                                     >  2:4090699/1‑93 (MQ=255)

                                                      ttCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTcgcg                                   >  2:3665916/1‑93 (MQ=255)

                                                      ttCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTcgcg                                   >  1:4180303/1‑93 (MQ=255)

                                                       tCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCgg                                  >  2:1275959/1‑93 (MQ=255)

                                                         tGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCa                                >  2:3897379/1‑93 (MQ=255)

                                                          ggCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAg                               >  2:4522253/1‑93 (MQ=255)

                                                          ggCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAg                               >  1:345337/1‑93 (MQ=255)

                                                          ggCCAGACCAGACTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAg                               >  2:613428/1‑93 (MQ=255)

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

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

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

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

                                                                 ccagacTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCtt                        >  2:3814646/1‑93 (MQ=255)

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

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

                                                                     acTCGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAg                    <  1:2344905/93‑1 (MQ=255)

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

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

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

                                                                        cGCCAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCAt                 >  1:4673101/1‑93 (MQ=255)

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

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

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

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

                                                                           cAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAg              <  2:1763544/93‑1 (MQ=255)

                                                                           cAGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAg              >  1:3161243/1‑93 (MQ=255)

                                                                            aGCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGc             >  1:1861812/1‑93 (MQ=255)

                                                                             gCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGcc            >  2:3009928/1‑93 (MQ=255)

                                                                             gCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGcc            <  1:3092899/93‑1 (MQ=255)

                                                                             gCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGcc            >  1:756478/1‑93 (MQ=255)

                                                                             gCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGcc            <  1:44215/93‑1 (MQ=255)

                                                                             gCAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGAATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGcc            <  1:1922377/93‑1 (MQ=255)

                                                                              cAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCa           >  1:616297/1‑93 (MQ=255)

                                                                              cAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCa           >  2:860549/1‑93 (MQ=255)

                                                                              cAAAAAAGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCa           <  2:4023066/93‑1 (MQ=255)

                                                                               aaaaaaGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAg          <  2:4454891/93‑1 (MQ=255)

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

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

                                                                                  aaaGCAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGgag       >  2:1344589/1‑93 (MQ=255)

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

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

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

                                                                                      cAAGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACa   <  2:1420224/93‑1 (MQ=255)

                                                                                       aaGCGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAg  >  2:197302/1‑93 (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` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`