BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,923,978	C→G	5.9%	E46Q (GAA→CAA)	yqcC ←	DUF446 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,923,978	0	C	G	5.9%	356.0 / 11.1	101	E46Q (GAA→CAA)	yqcC	DUF446 family protein
Reads supporting (aligned to +/- strand): ref base C (47/48); new base G (3/3); total (50/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

Only 100 of 102 total aligned reads displayed.

GCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCT  >  NC_000913/2923830‑2924123

gCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTc                                                                                                                                                   >  2:254741/1‑149 (MQ=255)

  cAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCa                                                                                                                                                 >  3:255631/1‑149 (MQ=255)

   aGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGTTGCag                                                                                                                                                <  1:260703/149‑2 (MQ=255)

   aGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTGGAGTGGTTGCag                                                                                                                                                <  4:11223/149‑2 (MQ=255)

   aGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCAt                                                                                                                                                <  1:20563/149‑1 (MQ=255)

    gAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATg                                                                                                                                               <  1:259699/149‑1 (MQ=255)

    gAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATg                                                                                                                                               <  2:508330/149‑1 (MQ=255)

         aGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCACAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTc                                                                                                                                          <  3:40201/149‑1 (MQ=255)

          gcgcgcGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTcc                                                                                                                                         <  2:693384/149‑1 (MQ=255)

            gcgcgTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCTCTCCGGTGGTTCCATGGTGTCCAt                                                                                                                                       <  2:550331/149‑1 (MQ=255)

                    cGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACg                                                                                                                               <  4:688404/149‑1 (MQ=255)

                    cGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACg                                                                                                                               <  2:671812/149‑1 (MQ=255)

                     ggATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACgg                                                                                                                              >  2:512980/1‑149 (MQ=255)

                       aTGGTCCGTCGCCAGCGCCATTTCATAGTAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAGCCACTCCAGTGGTTCCAGGGTGTCCATAAAGAACGGtt                                                                                                                            <  2:4220/149‑1 (MQ=255)

                       aTGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGtt                                                                                                                            <  3:371787/149‑1 (MQ=255)

                         ggTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTgg                                                                                                                          <  3:876964/149‑1 (MQ=255)

                           tCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGt                                                                                                                        >  4:313102/1‑149 (MQ=255)

                            ccGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTa                                                                                                                       >  4:194228/1‑149 (MQ=255)

                              gTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACt                                                                                                                     >  4:472735/1‑149 (MQ=255)

                               tCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTa                                                                                                                    <  4:464597/149‑1 (MQ=255)

                                  ccAGCGCCATTTCATAATAGGGTGCAACGACAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTa                                                                                                                 >  3:781555/1‑149 (MQ=255)

                                    aGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTaaa                                                                                                               <  2:777756/149‑1 (MQ=255)

                                    aGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTaaa                                                                                                               <  4:644643/149‑1 (MQ=255)

                                     gcgcCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAAt                                                                                                              <  4:255631/149‑1 (MQ=255)

                                     gcgcCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAAt                                                                                                              <  4:231268/149‑1 (MQ=255)

                                      cgcCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAAtt                                                                                                             >  1:452321/1‑149 (MQ=255)

                                      cgcCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAAtt                                                                                                             >  4:452693/1‑149 (MQ=255)

                                        ccATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGa                                                                                                           <  3:194228/149‑1 (MQ=255)

                                         cATTTCATAATAGGGTGCAACGGCAAAAGCGCCTGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGAt                                                                                                          >  3:672316/1‑149 (MQ=255)

                                         cATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGAt                                                                                                          <  2:711128/149‑1 (MQ=255)

                                         cATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGAt                                                                                                          <  3:338750/149‑1 (MQ=255)

                                           tttCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGc                                                                                                        >  2:305953/1‑149 (MQ=255)

                                                taataGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTg                                                                                                   >  2:894928/1‑149 (MQ=255)

                                                 aataGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTgg                                                                                                  <  1:305953/149‑1 (MQ=255)

                                                  ataGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTggg                                                                                                 <  2:58523/149‑1 (MQ=255)

                                                  ataGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTggg                                                                                                 <  3:358830/149‑1 (MQ=255)

                                                   taGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTgggg                                                                                                >  3:514456/1‑149 (MQ=255)

                                                           aaCGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGt                                                                                        <  2:751180/149‑1 (MQ=255)

                                                                 aaaaGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGcc                                                                                  >  1:411504/1‑149 (MQ=255)

                                                                 aaaaGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGcc                                                                                  >  3:441897/1‑149 (MQ=255)

                                                                    aGCGCCTGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGt                                                                               <  4:672316/149‑1 (MQ=255)

                                                                    aGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGt                                                                               <  3:622837/149‑1 (MQ=255)

                                                                     gcgcCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTg                                                                              <  3:313102/149‑1 (MQ=255)

                                                                     gcgcCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTg                                                                              <  1:512980/149‑1 (MQ=255)

                                                                       gcCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCt                                                                            <  1:254741/149‑1 (MQ=255)

                                                                        cccGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTg                                                                           >  1:651188/1‑149 (MQ=255)

                                                                        cccGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTg                                                                           >  2:303703/1‑149 (MQ=255)

                                                                          cGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGAt                                                                         <  2:441927/149‑1 (MQ=255)

                                                                            gTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGt                                                                       >  3:174198/1‑149 (MQ=255)

                                                                               aCGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGAGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCa                                                                    >  4:605434/1‑149 (MQ=255)

                                                                                    tgtttgttATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAg                                                               <  4:174198/149‑1 (MQ=255)

                                                                                        tgttATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAAc                                                           <  2:411504/149‑1 (MQ=255)

                                                                                         gttATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACg                                                          <  3:452693/149‑1 (MQ=255)

                                                                                          ttATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCGTGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGc                                                         >  3:729812/1‑149 (MQ=255)

                                                                                            aTCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCtt                                                       >  1:493668/1‑149 (MQ=255)

                                                                                               tAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCCTCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCaa                                                    <  4:536350/149‑1 (MQ=255)

                                                                                               tAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCaa                                                    >  4:585169/1‑149 (MQ=255)

                                                                                               tAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCaa                                                    >  2:292757/1‑149 (MQ=255)

                                                                                                aaCAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAg                                                   >  1:784605/1‑149 (MQ=255)

                                                                                                aaCAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAg                                                   >  4:555799/1‑149 (MQ=255)

                                                                                                aaCAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAg                                                   >  1:391004/1‑149 (MQ=255)

                                                                                                 aCAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAgc                                                  >  2:776219/1‑149 (MQ=255)

                                                                                                     aTCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCt                                              >  1:849356/1‑149 (MQ=255)

                                                                                                     aTCGTGCATTCGAGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCt                                              <  2:535169/149‑1 (MQ=255)

                                                                                                          gCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAAc                                         <  1:292757/149‑1 (MQ=255)

                                                                                                            aTTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTgg                                                                                                                          <  2:715350/66‑1 (MQ=255)

                                                                                                            aTTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTgg                                                                                                                          >  1:715350/1‑66 (MQ=255)

                                                                                                             ttCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTgg                                      >  3:654395/1‑149 (MQ=255)

                                                                                                              tCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGga                                     <  4:440681/149‑1 (MQ=255)

                                                                                                                  gggATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACg                                 >  2:413587/1‑149 (MQ=255)

                                                                                                                    gATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGtt                                                                                       <  3:790294/93‑1 (MQ=255)

                                                                                                                    gATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGtt                                                                                       >  4:790294/1‑93 (MQ=255)

                                                                                                                     aTCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGt                                                                                               >  3:17897/1‑84 (MQ=255)

                                                                                                                     aTCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGt                                                                                               <  4:17897/84‑1 (MQ=255)

                                                                                                                         gCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCgg                          >  2:930463/1‑149 (MQ=255)

                                                                                                                             ccACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCAt                      >  4:310440/1‑149 (MQ=255)

                                                                                                                             ccACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCAt                      >  2:190562/1‑149 (MQ=255)

                                                                                                                               catggaaCCACTGGAGTGGTTGCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGa                    >  1:128465/6‑149 (MQ=255)

                                                                                                                               aCTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGa                    <  2:784605/149‑1 (MQ=255)

                                                                                                                                atggaaCCACTGGAGTGGTTGCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTc                                                      >  4:508831/5‑114 (MQ=255)

                                                                                                                                atggaaCCACTGGAGTGGTTGCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTc                                                      <  3:508831/110‑1 (MQ=255)

                                                                                                                                atggaaCCACTGGAGTGGTTGCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAg                   >  1:281516/5‑149 (MQ=255)

                                                                                                                                atggaaCCACTGGAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGga                                     <  2:346264/127‑1 (MQ=255)

                                                                                                                                atggaaCCACTGGAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGga                                     >  1:346264/5‑131 (MQ=255)

                                                                                                                                 tGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGt                  >  1:598963/1‑149 (MQ=255)

                                                                                                                                  gCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTg                 <  1:894928/149‑1 (MQ=255)

                                                                                                                                   cAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTgg                >  1:607936/1‑149 (MQ=255)

                                                                                                                                    aaCCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGt               >  2:401551/1‑149 (MQ=255)

                                                                                                                                        aCTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATa           >  1:751290/1‑149 (MQ=255)

                                                                                                                                          tCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCAt            <  3:822518/146‑1 (MQ=255)

                                                                                                                                          tCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATaaa         <  3:28896/149‑1 (MQ=255)

                                                                                                                                            cAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGc       >  3:827987/1‑149 (MQ=255)

                                                                                                                                            cAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGc       >  1:356340/1‑149 (MQ=255)

                                                                                                                                               tGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAAc    >  2:424905/1‑149 (MQ=255)

                                                                                                                                               tGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAAc    <  2:607936/149‑1 (MQ=255)

                                                                                                                                                 gTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATaaa         <  3:364465/142‑1 (MQ=255)

                                                                                                                                                 gTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATaaa         <  3:385726/142‑1 (MQ=255)

                                                                                                                                                 gTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATaaa         >  4:385726/1‑142 (MQ=255)

                                                                                                                                                 gTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATaaa         >  4:364465/1‑142 (MQ=255)

                                                                                                                                                 gTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCt  <  1:424905/149‑1 (MQ=255)

GCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCT  >  NC_000913/2923830‑2924123

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