BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,643,234	C→T	6.9%	G650D (GGC→GAC) ‡	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,643,234	0	C	T	6.9%	309.7 / 8.8	101	G650D (GGC→GAC) ‡	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand): ref base C (43/50); new base T (2/5); total (45/56)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

Only 100 of 104 total aligned reads displayed.

TAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATG  >  NC_000913/3643085‑3643382

tAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGc                                                                                                                                                      >  2:1021238/1‑150 (MQ=255)

 aaGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGcc                                                                                                                                                     >  2:1690323/1‑150 (MQ=255)

   ggCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCCGAACCGCCAc                                                                                                                                                   >  2:907023/1‑150 (MQ=255)

         ggagacGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCa                                                                                                                                             >  1:686518/6‑150 (MQ=255)

          tctcCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAg                                                                                                                                            >  1:1687488/1‑150 (MQ=255)

           ctcCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGa                                                                                                                                           <  1:1137058/150‑1 (MQ=255)

           ctcCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCACCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCAGGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGCACCGCCACGGCTCAGa                                                                                                                                           >  1:230885/1‑150 (MQ=255)

             ccGGTGCCTGTTTCATCAATTAATCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAGCCGTGGCGGTTCAGaag                                                                                                                                         <  2:1653989/150‑2 (MQ=255)

             ccGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGTCTCTTCTGAACCGCCACGGCTCAGAAt                                                                                                                                         <  1:711306/150‑1 (MQ=255)

               gtaggcTGTTTCATCAATTAAGCAGATTTTCACTGATTGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAAtgt                                                                                                                                       <  2:124774/145‑1 (MQ=255)

                gTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAGCCGTGGCGGTTCAGaagagc                                                                                                                                      <  1:449/150‑5 (MQ=255)

                 tGCCTGTTTCATCTATTAAGCAGATTTTCACTGAATGAGCAGCCCTTAATGCCGTAATGCTCCAGCATCTCATCCAGCTGCGGTTCACGACCACGCAAGCGTTTGAACAGATCCATCGGCTCTTCTGAGCCGTGGCGGTTCAGaagagcc                                                                                                                                     <  2:1506532/150‑6 (MQ=255)

                 tGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAGCCGTGGCGGTTCAGaagagcc                                                                                                                                     <  2:910639/150‑6 (MQ=255)

                 tGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAGCCGTGGCGGTTCAGaagagcc                                                                                                                                     <  1:1628859/150‑6 (MQ=255)

                 tGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGAGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCAACGGCTCAGAAttttt                                                                                                                                     >  2:1497405/1‑149 (MQ=255)

                      gTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAg                                                                                                                                >  1:525754/1‑150 (MQ=255)

                          catcaATTAAGCAGATTTTCACTGAATGATCCGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCAAGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGacc                                                                                                                            >  2:1258737/1‑148 (MQ=255)

                          catcaATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAAc                                                                                                                            >  1:1171566/1‑150 (MQ=255)

                           atcaATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACg                                                                                                                           <  2:1473951/150‑1 (MQ=255)

                             caATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGGAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCCCGGCTCAGAATGTTGTCGAGGAACGAc                                                                                                                         >  1:234612/1‑150 (MQ=255)

                              aaTTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGCACAGATCCATCGGCTCTTCTGAAACGCCACGGATCAGAATGTTGTCGAGGAACGCCt                                                                                                                        >  2:152447/1‑150 (MQ=255)

                              aaTTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACt                                                                                                                        >  2:907054/1‑150 (MQ=255)

                              aaTTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACt                                                                                                                        <  1:18132/150‑1 (MQ=255)

                              aaTTAACCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACt                                                                                                                        <  1:24156/150‑1 (MQ=255)

                               aTTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGGCGAGGAACGACTg                                                                                                                       <  2:49520/150‑1 (MQ=255)

                                 tAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACAACTGcc                                                                                                                     >  2:1381437/1‑150 (MQ=255)

                                   aGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCg                                                                                                                   <  1:675718/150‑1 (MQ=255)

                                    gCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTATGAACAGATCCATCGGCTCTACGGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCgg                                                                                                                  >  2:593419/1‑150 (MQ=255)

                                     cAGATTTTCACTTAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTATGAACCGCCACGTCTCAGAATGTTTTAGAGTCACGACTGCCCGGt                                                                                                                 >  1:1272845/1‑150 (MQ=255)

                                        aTTTTCACTGAATGCGCAGCCCTTAATGCCGTAATGCTCCAGCATCGCCTCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGACACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTc                                                                                                              <  1:471248/150‑1 (MQ=255)

                                        aTTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGGTCACGACCACGGAAGCGTTTGAACAGCTCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTc                                                                                                              <  1:471247/150‑1 (MQ=255)

                                          tttCACTGAATGATCCGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCAc                                                                                                            <  2:880658/150‑1 (MQ=255)

                                             cACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTTAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGt                                                                                                         >  2:1204370/1‑150 (MQ=255)

                                             cACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGt                                                                                                         >  2:642227/1‑150 (MQ=255)

                                              aCTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGtt                                                                                                        >  1:1348875/1‑150 (MQ=255)

                                                tGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGCACGACTGCCCGGTTTCACGGTTGa                                                                                                      >  1:1404803/1‑150 (MQ=255)

                                                tGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGa                                                                                                      <  2:1666372/150‑1 (MQ=255)

                                                  aaTGATCAGCCCTAAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTACAAGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGATCAGAATGATGTCGAGGAACGACTGCCCGGTTTCACGGTTGaaa                                                                                                    <  1:876068/150‑1 (MQ=255)

                                                      aTCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGCATGTTGTCGAGGAGCGACTGCCCGGTTAAGCGGTGGAAAATGc                                                                                                >  2:653150/1‑150 (MQ=255)

                                                       tCAGCCCTTACTGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGcc                                                                                               <  1:347752/150‑1 (MQ=255)

                                                       tCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGcc                                                                                               <  2:609885/150‑1 (MQ=255)

                                                       tCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGAGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGcc                                                                                               <  2:1056168/150‑1 (MQ=255)

                                                        cAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACAACTGCCCGGTTTCACGGTTGAAAATGccc                                                                                              <  1:1381437/150‑1 (MQ=255)

                                                         aGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCct                                                                                             <  2:1315746/150‑1 (MQ=255)

                                                            ccTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGCCTGCCCGGTTTCACGGTTGAAAATGCCCTCtt                                                                                          >  1:1551821/1‑150 (MQ=255)

                                                               tAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCt                                                                                       <  2:529687/150‑1 (MQ=255)

                                                               tAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCt                                                                                       <  2:460882/150‑1 (MQ=255)

                                                                    ccGTAATGCTCCCGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCATAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAg                                                                                  >  2:824813/1‑150 (MQ=255)

                                                                    ccGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCCTCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAATGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAg                                                                                  >  1:1078017/1‑150 (MQ=255)

                                                                    ccGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCATAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAg                                                                                  >  2:484874/1‑150 (MQ=255)

                                                                    ccGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAg                                                                                  <  1:510819/150‑1 (MQ=255)

                                                                      gTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAgcg                                                                                <  2:686518/150‑1 (MQ=255)

                                                                        aaTGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTGCACGGTTGAAAATGCCCTCTTCCTCAAAGcgcg                                                                              <  2:539391/150‑1 (MQ=255)

                                                                            cTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGaa                                                                          >  2:942043/1‑150 (MQ=255)

                                                                              ccAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGCACCGCCACGGCTCACAATGAAGTCGAGCAGCGACTGCCCGGATTCTCGGTTGAAAATGCCCACTACCTCAAAGCGCGAGAAAg                                                                        >  1:129903/1‑150 (MQ=255)

                                                                              ccAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAg                                                                        <  2:242815/150‑1 (MQ=255)

                                                                                   atcgcatcCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCt                                                                   <  1:1497405/150‑1 (MQ=255)

                                                                                    tcgcatcCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATTCACTCTTCCTCAAAGCGAGAGAAAGCAACTg                                                                  >  2:268260/1‑150 (MQ=255)

                                                                                    tcgcatcCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATATg                                                                  <  2:1153204/150‑1 (MQ=255)

                                                                                     cgcatcCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAACCGCGAGAAAGCATCTgc                                                                 >  2:83137/1‑150 (MQ=255)

                                                                                       catcCAGCTGCGGTTCACGACCACTGAAGAGTTTTAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGTTTTCACGGTTGAAATTGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTgccg                                                               <  1:833733/150‑1 (MQ=255)

                                                                                       catcCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTgccg                                                               <  1:1669656/150‑1 (MQ=255)

                                                                                       catcCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTgccg                                                               <  1:531870/150‑1 (MQ=255)

                                                                                                    ttCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGc                                                  >  1:772259/1‑150 (MQ=255)

                                                                                                    ttCACGACCACGGAAGCGGTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGc                                                  <  1:269086/150‑1 (MQ=255)

                                                                                                     tCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGcc                                                 <  1:580294/150‑1 (MQ=255)

                                                                                                      cACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCCCGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGccc                                                >  1:1329597/1‑150 (MQ=255)

                                                                                                        cGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCcac                                              >  2:1667505/1‑150 (MQ=255)

                                                                                                         gACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACTGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAACCcaca                                             <  2:1362314/150‑1 (MQ=255)

                                                                                                         gACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCcaca                                             >  1:1553807/1‑150 (MQ=255)

                                                                                                         gACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCcaca                                             <  1:1522943/150‑1 (MQ=255)

                                                                                                           ccACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAgg                                           <  1:1690323/150‑1 (MQ=255)

                                                                                                            cACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGt                                          <  1:441491/150‑1 (MQ=255)

                                                                                                                gAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCt                                      >  1:316226/1‑150 (MQ=255)

                                                                                                                 aaGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTg                                     <  2:1496480/150‑1 (MQ=255)

                                                                                                                  aGCGTTTGAACTGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTGGGCGAGGACCGACTGCCCGGTTTCACGGTTGAAAATGCCCTCATCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTgt                                    <  2:842131/150‑1 (MQ=255)

                                                                                                                     gTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGCGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCCCAGGTAGCTgtagt                                 >  1:1216710/1‑150 (MQ=255)

                                                                                                                        tGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAAc                              <  1:1297563/150‑1 (MQ=255)

                                                                                                                         gaccAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTTTCTAGTTACGACTGCCAGGTTACACGGTTTAATATGCCCTCTTCCTCAAAGCGGGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGAAGCTGTAGTAAcc                             <  1:1068374/147‑1 (MQ=255)

                                                                                                                         gAACAGATCCATCGGCTCTTCTGAACCGCCCCGGCTCAGAATGTTGTCGGGGCAGGACTGCCGGGTATCCGGGTTGAAAATGCCCACTTCCTCAAAGCGCGAGCAAGCATCTGCCGCCAGTACGACAGCCCCCAGGTAGCAGTATTAccc                             >  2:242801/1‑150 (MQ=255)

                                                                                                                          aaCAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCt                            <  2:1551821/150‑1 (MQ=255)

                                                                                                                          aaCAGATCCAACGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAAccg                            >  1:1443213/1‑149 (MQ=255)

                                                                                                                            cAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGc                          >  2:1527546/1‑150 (MQ=255)

                                                                                                                             aGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGCCTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCg                         >  1:514346/1‑150 (MQ=255)

                                                                                                                               aaccATCGGCTCTTCTGAACCGCCACGGCTCAGAATTTTGTCGACGAACGACTGCCCGGTTTCACGGGTGAAAATGCCCTCTTCCTCAAAGCGCGAGATCGCATCTGCCGCCATTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGc                       <  2:234612/148‑1 (MQ=255)

                                                                                                                                tCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCa                      <  2:906951/150‑1 (MQ=255)

                                                                                                                                 ccATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCTCCAGTACGTCAGCCCACAGGTATCTGTTGTAACCTCCGGCAt                     >  2:1052067/1‑150 (MQ=255)

                                                                                                                                 ccATCGGCTCTTCCGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCAt                     <  1:907023/150‑1 (MQ=255)

                                                                                                                                   aTCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATaa                   <  1:679139/150‑1 (MQ=255)

                                                                                                                                   aTCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATaa                   <  2:1288589/150‑1 (MQ=255)

                                                                                                                                    tCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAac                  <  2:1652523/150‑1 (MQ=255)

                                                                                                                                    tCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAac                  <  1:907054/150‑1 (MQ=255)

                                                                                                                                     cGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGATAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAacc                 >  2:352308/1‑150 (MQ=255)

                                                                                                                                     cGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAacc                 <  1:973317/150‑1 (MQ=255)

                                                                                                                                         tcttctGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCg             <  2:1348875/150‑1 (MQ=255)

                                                                                                                                          attctGAGCCGTGGCGGTTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCgg            >  1:1412425/2‑150 (MQ=255)

                                                                                                                                           ttctGAGCCGTGGCGGTTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCCCAGGTAGCTGTACTAACCTGCGGCATAACCCCCGGc           >  1:837760/1‑150 (MQ=255)

                                                                                                                                              tGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCCCAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGaa        >  1:1557860/1‑150 (MQ=255)

                                                                                                                                                aaCCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGaaaa      <  2:1687488/150‑1 (MQ=255)

                                                                                                                                                    gCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCAGAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATg  >  2:213194/1‑150 (MQ=255)

TAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATG  >  NC_000913/3643085‑3643382

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`