BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	440,662	C→A	5.3%	S147I (AGC→ATC) ‡	ispA ←	geranyltranstransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	440,662	0	C	A	5.3%	322.8 / 5.3	103	S147I (AGC→ATC) ‡	ispA	geranyltranstransferase
Reads supporting (aligned to +/- strand): ref base C (46/51); new base A (3/3); total (49/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

CGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAA  >  NC_000913/440513‑440811

cgcTTAATGCCCCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGTCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGgcgc                                                                                                                                                       >  2:249657/1‑150 (MQ=255)

cgcTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGgcgc                                                                                                                                                       <  2:1219497/150‑1 (MQ=255)

 gcTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCt                                                                                                                                                      >  1:685318/1‑150 (MQ=255)

    tAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGc                                                                                                                                                   <  1:596483/150‑1 (MQ=255)

       tGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAg                                                                                                                                                >  1:1261833/1‑150 (MQ=255)

       tGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAg                                                                                                                                                <  2:595515/150‑1 (MQ=255)

        gCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGt                                                                                                                                               >  1:799324/1‑150 (MQ=255)

         cACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtg                                                                                                                                              <  2:1321411/150‑2 (MQ=255)

          aCCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTc                                                                                                                                             <  1:4084/150‑1 (MQ=255)

            cAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtggta                                                                                                                                           <  1:437963/150‑5 (MQ=255)

            cAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtggta                                                                                                                                           <  1:936471/150‑5 (MQ=255)

             aaGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGa                                                                                                                                          <  2:537050/150‑1 (MQ=255)

                gCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCCCCGCCCATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGCAat                                                                                                                                       >  1:1399643/1‑150 (MQ=255)

                    aCGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAatcatc                                                                                                                                   <  1:1726560/150‑1 (MQ=255)

                     cggcggCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCg                                                                                                                                  <  2:1115535/150‑1 (MQ=255)

                     cggcggCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCg                                                                                                                                  >  2:1687663/1‑150 (MQ=255)

                     cggcggCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCg                                                                                                                                  <  2:331306/150‑1 (MQ=255)

                      ggcggcGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGATATCATCGa                                                                                                                                 >  1:1009266/1‑150 (MQ=255)

                           cgcgAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTc                                                                                                                            <  1:1466967/150‑1 (MQ=255)

                           cgcgAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTGCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTc                                                                                                                            >  2:1372026/1‑150 (MQ=255)

                           cgcgAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGAGTCCAGAGGTACGTGATTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTc                                                                                                                            >  1:1801037/1‑150 (MQ=255)

                               aaTCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTc                                                                                                                        <  1:404756/150‑1 (MQ=255)

                                aTCAATGCGCCGGTATTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTcg                                                                                                                       >  1:1442306/1‑150 (MQ=255)

                                   aaTGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCgg                                                                                                                    >  1:596831/1‑150 (MQ=255)

                                    aTGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGt                                                                                                                   <  1:299121/150‑1 (MQ=255)

                                                tATGACGATGAATACGCTCAAGCGCGGCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCg                                                                                                       <  1:42989/150‑1 (MQ=255)

                                                 aTGACGATCAATACGCTCAAGCGCGTCCAGAGGAACGTGTTCGCCTGCCGCGTCTAAATCTACTGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCgg                                                                                                      <  2:1392369/150‑1 (MQ=255)

                                                   gACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCTGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGACa                                                                                                    <  1:241514/150‑1 (MQ=255)

                                                         gAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCgg                                                                                              >  1:958980/1‑150 (MQ=255)

                                                          aaTACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGc                                                                                             <  2:844296/150‑1 (MQ=255)

                                                          aaTACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGc                                                                                             <  2:837007/150‑1 (MQ=255)

                                                             aCGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATc                                                                                          <  2:559329/150‑1 (MQ=255)

                                                             aCGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATc                                                                                          >  1:1659917/1‑150 (MQ=255)

                                                             aCGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACAACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGGCCGACACTTCCGGCATATCGGCATc                                                                                          <  2:410923/150‑1 (MQ=255)

                                                                cTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCTATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCGGGCATATCGGCATCGCt                                                                                       >  2:1473380/1‑150 (MQ=255)

                                                                cTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCt                                                                                       >  2:1473379/1‑150 (MQ=255)

                                                                cTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCt                                                                                       >  2:1480626/1‑150 (MQ=255)

                                                                cTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCt                                                                                       <  2:1616587/150‑1 (MQ=255)

                                                                 tCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCtt                                                                                      <  2:1286969/150‑1 (MQ=255)

                                                                   aaGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTaa                                                                                    <  1:1743608/150‑1 (MQ=255)

                                                                     gcgcGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTaaaa                                                                                  >  2:1307642/1‑150 (MQ=255)

                                                                      cgcgTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCTAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAAt                                                                                 <  1:1355232/150‑1 (MQ=255)

                                                                        cgTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGGAATACCACTGGCGCTCGCCAGTTCTGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCAGCGCTTAAAATCg                                                                               <  1:1229477/150‑1 (MQ=255)

                                                                          tCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGCCCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAAGCgag                                                                             >  2:832022/1‑150 (MQ=255)

                                                                           ccAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCgaga                                                                            <  1:843627/150‑1 (MQ=255)

                                                                                ggTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCGAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGcc                                                                       <  1:156342/150‑1 (MQ=255)

                                                                                ggTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGcc                                                                       >  2:760946/1‑150 (MQ=255)

                                                                                  tACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAg                                                                     <  2:1678488/150‑1 (MQ=255)

                                                                                    cGTGTTTGCCTTCCGCGTCTAACTCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCg                                                                   <  2:1461999/150‑1 (MQ=255)

                                                                                    cGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCg                                                                   >  1:1660523/1‑150 (MQ=255)

                                                                                    cGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCAGTCTATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGCCACGTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCg                                                                   >  1:305865/1‑150 (MQ=255)

                                                                                       gtTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGttt                                                                <  2:1376515/150‑1 (MQ=255)

                                                                                        tttGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTg                                                               <  1:1691212/150‑1 (MQ=255)

                                                                                        tttGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCACATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTg                                                               >  1:1711706/1‑150 (MQ=255)

                                                                                         ttGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGt                                                              <  1:1397370/150‑1 (MQ=255)

                                                                                          tGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAATTCGAGAACGCCAGCGTTTGTa                                                             >  1:1677800/1‑150 (MQ=255)

                                                                                            ccTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACTCCCGCGTTTGTaaa                                                           >  1:1695155/1‑150 (MQ=255)

                                                                                             cTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAg                                                          >  2:1553723/1‑150 (MQ=255)

                                                                                                 cgcgTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTc                                                      >  2:1858248/1‑150 (MQ=255)

                                                                                                 cgcgTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTc                                                      >  1:670381/1‑150 (MQ=255)

                                                                                                   cgTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGc                                                    <  1:249657/150‑1 (MQ=255)

                                                                                                   cgTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCGCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGc                                                    <  2:809645/150‑1 (MQ=255)

                                                                                                      cTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGGCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAg                                                 <  2:685318/150‑1 (MQ=255)

                                                                                                       tAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGc                                                >  1:1250167/1‑150 (MQ=255)

                                                                                                        aaaTCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCg                                               >  1:444531/1‑150 (MQ=255)

                                                                                                          aTCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCgag                                             >  2:878010/1‑150 (MQ=255)

                                                                                                          aTCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCACCAGCgag                                             <  2:452971/150‑1 (MQ=255)

                                                                                                           tCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCGTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCgaga                                            >  1:1500945/1‑150 (MQ=255)

                                                                                                            cTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGaa                                           >  2:473658/1‑150 (MQ=255)

                                                                                                                tGCCTGACCCCCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCATAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTCAAGCGTCGCCAGCGAGAATcgc                                       >  2:1540273/1‑150 (MQ=255)

                                                                                                                tGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATcgc                                       >  1:1351344/1‑150 (MQ=255)

                                                                                                                tGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATcgc                                       >  2:1049628/1‑150 (MQ=255)

                                                                                                                tGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATcgc                                       >  1:1614361/1‑150 (MQ=255)

                                                                                                                     gACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTg                                  >  2:1514591/1‑150 (MQ=255)

                                                                                                                     gACCACCGCACATTCCGGCAATACCACTGGCGATCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGATCCGCCACTTCCGGCATATCAGCATAGCTTAAAATCGAGAACGCCAGCGTTTCTAAAGCGTCGCCAGCGAGAATCGCGTATg                                  >  1:325316/1‑150 (MQ=255)

                                                                                                                      accaccGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGc                                 <  2:720094/150‑1 (MQ=255)

                                                                                                                         accGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTc                              <  2:1769992/150‑1 (MQ=255)

                                                                                                                         accGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTc                              >  2:1733855/1‑150 (MQ=255)

                                                                                                                          ccGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCg                             <  1:460529/150‑1 (MQ=255)

                                                                                                                          ccGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCg                             <  1:1858248/150‑1 (MQ=255)

                                                                                                                          ccGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGAATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCAACAGAGAGAATCGCGTTTGCTTCg                             >  1:929306/1‑150 (MQ=255)

                                                                                                                            gCACATTCCGGCAATACCTCTGGCGCTCGCCAGTTCAGACATCATCGACATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCGTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGcc                           <  2:307002/150‑1 (MQ=255)

                                                                                                                            gCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGcc                           <  2:817125/150‑1 (MQ=255)

                                                                                                                            gCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGcc                           >  1:558365/1‑150 (MQ=255)

                                                                                                                             cacaTTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGccc                          >  2:1851297/1‑149 (MQ=255)

                                                                                                                               caTTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGCCACTGCCGCCATATCGGCATCGCTTAAAAGAGAGACCGCCAGCGGTTGTAAAGCGTCGCCCGCGAGAATCGCGTTTGCTTCGACata                        >  1:239243/1‑148 (MQ=255)

                                                                                                                                 ttCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGTCCTATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCTCCCTCGAGAATCGCGTTAGCTTCGCCAAACt                      >  1:278122/1‑150 (MQ=255)

                                                                                                                                 ttCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACt                      >  2:1572364/1‑150 (MQ=255)

                                                                                                                                 ttCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACt                      <  2:597526/150‑1 (MQ=255)

                                                                                                                                   ccGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCCGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCACACtgc                    >  1:1724588/1‑148 (MQ=255)

                                                                                                                                   ccGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTc                    <  2:870188/150‑1 (MQ=255)

                                                                                                                                   ccGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTc                    <  2:877809/150‑1 (MQ=255)

                                                                                                                                    cGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTca                   <  1:182190/150‑1 (MQ=255)

                                                                                                                                     ggCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTcac                  >  2:1900419/1‑150 (MQ=255)

                                                                                                                                          tACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGc             <  2:1085292/150‑1 (MQ=255)

                                                                                                                                           tgaaCTGGCGAGCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCa            >  2:1880018/4‑150 (MQ=255)

                                                                                                                                            ccACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAg           >  2:153954/1‑150 (MQ=255)

                                                                                                                                                  gcgcTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTGGCCTAACTTCACATGGCAGGTTGGc     <  2:459326/150‑1 (MQ=255)

                                                                                                                                                    gcTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCaa   <  1:1309103/150‑1 (MQ=255)

                                                                                                                                                     cTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCTCCAAACTTCACATGGCAGGTTGGCaaa  >  1:1882453/1‑150 (MQ=255)

CGCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCAAA  >  NC_000913/440513‑440811

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