BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,366,854	T→G	8.2%	intergenic (‑81/+37)	dcuA ← / ← aspA	C4‑dicarboxylate antiporter/aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,366,854	0	T	G	8.2%	634.9 / 46.5	197	intergenic (‑81/+37)	dcuA/aspA	C4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand): ref base T (92/89); new base G (6/10); total (98/99)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

Only 100 of 198 total aligned reads displayed.

ATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTA  >  NC_000913/4366586‑4367088

agggcgaaaaTAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGGTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACg                                                                                                                                                                                                                                      <  4:387949/268‑1 (MQ=255)

   gCAGAAAGAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACTGGATGGTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCCCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCttt                                                                                                                                                                                                                               <  1:208880/279‑1 (MQ=255)

                                                           aCCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGGACCTTTTTACTTTGTACTACCCTGTCCGATTACTGTTCGCTTTCATCCTTAAAGCGtt                                                                                                                                                                           >  2:34825/1‑275 (MQ=255)

                                                                                                                         cTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATc                                                                                                                                                                                                                                           >  7:852335/1‑149 (MQ=255)

                                                                                                                         cTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATc                                                                                                                                                                                                                                           <  8:532819/149‑1 (MQ=255)

                                                                                                                            ttCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCt                                                                                                                                                                                                                                          >  2:153794/1‑147 (MQ=255)

                                                                                                                            ttCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGa                                                                                                                                                                                                                                        >  5:1307417/1‑149 (MQ=255)

                                                                                                                               cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGt                                                                                                                                                                                                                                     <  5:1163178/149‑1 (MQ=255)

                                                                                                                                ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                                                                                                                    >  7:724420/1‑149 (MQ=255)

                                                                                                                                 cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGc                                                                                                                                                                                                                                   >  6:1380662/1‑149 (MQ=255)

                                                                                                                                 cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGc                                                                                                                                                                                                                                   >  7:913132/1‑149 (MQ=255)

                                                                                                                                 cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCt                                                                                                                                                                                            >  4:25375/1‑188 (MQ=255)

                                                                                                                                 cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGc                                                                                                                                                                                                                                   <  6:1265252/149‑1 (MQ=255)

                                                                                                                                    aaTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCtt                                                                                                                                                                                                                                <  6:723509/149‑1 (MQ=255)

                                                                                                                                        tGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCt                                                                                                                                                                                                                                          >  6:797439/1‑135 (MQ=255)

                                                                                                                                        tGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCt                                                                                                                                                                                                                                          <  5:797439/135‑1 (MQ=255)

                                                                                                                                           cgcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtt                                                                                                                                                                                                                         <  8:63217/149‑1 (MQ=255)

                                                                                                                                            gcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAATAAAAAAGGCACGTCAGATGACGTGCCTTTTTTCtt                                                                                                                                                                                                                         <  8:1517993/149‑1 (MQ=38)

                                                                                                                                               cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGta                                                                                                                                                                                                                     >  7:540323/1‑149 (MQ=255)

                                                                                                                                                aagaagATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtac                                                                                                                                                                                                                    >  8:826865/1‑149 (MQ=255)

                                                                                                                                                aagaagATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtac                                                                                                                                                                                                                    >  5:510654/1‑149 (MQ=255)

                                                                                                                                                  gaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtacta                                                                                                                                                                                                                  >  8:885625/1‑149 (MQ=255)

                                                                                                                                                    agaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAcc                                                                                                                                                                                                                >  5:1423571/1‑149 (MQ=255)

                                                                                                                                                      aTCGCCAGCAAAACTATGATGAATTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTGTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGt                                                                            <  1:166288/279‑1 (MQ=255)

                                                                                                                                                         gCCAACAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTa                                                                                                                                                                                                           <  7:585561/149‑1 (MQ=255)

                                                                                                                                                             gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAt                                                                                                                                                                                                                          >  5:907244/1‑130 (MQ=255)

                                                                                                                                                             gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAt                                                                                                                                                                                                       <  7:407746/149‑1 (MQ=255)

                                                                                                                                                                   cTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGt                                                                                                                                                                                                 <  8:1494019/149‑1 (MQ=255)

                                                                                                                                                                    tATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCtttttt                                                                                                                                                                                                                            >  5:601016/1‑121 (MQ=255)

                                                                                                                                                                    tATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCtttttt                                                                                                                                                                                                                            <  6:601016/121‑1 (MQ=255)

                                                                                                                                                                       gatgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGc                                                                                                                                                                                             >  7:229283/1‑149 (MQ=255)

                                                                                                                                                                        atgaGTTCTACAACTAGCATCTTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTc                                                                                                                                                                                               <  8:1174250/146‑1 (MQ=255)

                                                                                                                                                                        atgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTGAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTc                                                                                                                                                                                               >  7:1174250/1‑146 (MQ=255)

                                                                                                                                                                        atgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCt                                                                                                                                                                                            <  5:928952/149‑1 (MQ=255)

                                                                                                                                                                          gaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAACCTGTACGATTACTGTTCGCttt                                                                                                                                                                                          >  5:1524327/1‑149 (MQ=255)

                                                                                                                                                                           aGTTCTACAACTTGCATATTAGCCGTCCTTGTTTTTTAACAAGTTGTTATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTGATTTGTACTACCCTGTACGATTGCTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTa                                                                                                                                          <  4:475712/196‑1 (MQ=255)

                                                                                                                                                                            gTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                                                                                                                        >  5:1661926/1‑149 (MQ=255)

                                                                                                                                                                             ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAttt                                                                                                                                                                                                                        >  8:722733/1‑116 (MQ=255)

                                                                                                                                                                             ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtt                                                                                                                                                                                                                         >  2:723051/1‑115 (MQ=255)

                                                                                                                                                                             ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAccc                                                                                                                                                                                                               >  1:37032/1‑125 (MQ=255)

                                                                                                                                                                             ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat                                                                                                                                                                                       >  6:88895/1‑149 (MQ=255)

                                                                                                                                                                             ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGaa                                                                                       >  4:202779/1‑245 (MQ=255)

                                                                                                                                                                                tACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAg                                                                                                                                                                                    >  7:144845/1‑149 (MQ=255)

                                                                                                                                                                                 aCAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGta                                                                                                                                                                                                                     <  4:764079/115‑1 (MQ=255)

                                                                                                                                                                                   aaCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGtat                                                                                                                                                                                 >  6:1190574/1‑149 (MQ=255)

                                                                                                                                                                                    aCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCtttt                                                                                                                                                                                                                              <  5:677027/103‑1 (MQ=255)

                                                                                                                                                                                    aCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCtttt                                                                                                                                                                                                                              >  6:677027/1‑103 (MQ=255)

                                                                                                                                                                                      tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGc                                                                                                                                                                              >  8:350959/1‑149 (MQ=255)

                                                                                                                                                                                      tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGc                                                                                                                                                                              >  5:827281/1‑149 (MQ=255)

                                                                                                                                                                                       aGCATGTTAGCCTTACTTGTGTTTTAACAAGTTGTTATTAGATTGTTACTTTTACGTTGCTGCTCAC‑AAGAAAAAAAGCGCGTCATGTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTGTGCTTTGTAAGACGGGTGCATCAGATTCTGTCCGGAGAAAATTTCGTCAAGTTCGGCTCCAGTCAACAGCCCGTGTTCCAGAACGATTTCCCTAACACTCTTACCGGTTACGGCACAGtttt                                                >  4:508318/1‑274 (MQ=255)

                                                                                                                                                                                        gCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGa                                                                                                                                                                                                                                        >  5:775437/1‑89 (MQ=255)

                                                                                                                                                                                           tatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGta                                                                                                                                                                                                                     <  4:740526/105‑1 (MQ=255)

                                                                                                                                                                                             ttAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGc                                                                                                                                                                       >  7:730429/1‑149 (MQ=255)

                                                                                                                                                                                                gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                                                                                                                    <  5:1351651/149‑1 (MQ=255)

                                                                                                                                                                                                gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                                                                                                                    >  8:1216601/1‑149 (MQ=255)

                                                                                                                                                                                                gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                                                                                                                    <  6:1311811/149‑1 (MQ=255)

                                                                                                                                                                                                       ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtac                                                                                                                                                                                                                    >  6:493574/1‑94 (MQ=255)

                                                                                                                                                                                                       ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtac                                                                                                                                                                                                                    <  5:493574/94‑1 (MQ=255)

                                                                                                                                                                                                       ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCTTCAGTATAGCGTTTTGCTTTGTAAGcc                                                                                                                                                             >  5:711355/1‑149 (MQ=255)

                                                                                                                                                                                                       ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGcc                                                                                                                                                             <  6:637051/149‑1 (MQ=255)

                                                                                                                                                                                                         gTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑CAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGtt                             <  2:51221/275‑1 (MQ=255)

                                                                                                                                                                                                          ttttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTAc                                                                                                                                         >  3:256614/1‑166 (MQ=255)

                                                                                                                                                                                                           tttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGt                                                                                                                                                         <  6:1554816/149‑1 (MQ=255)

                                                                                                                                                                                                            ttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTg                                                                                                                                                        >  5:120690/1‑149 (MQ=255)

                                                                                                                                                                                                               tAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATc                                                                                                                                                                                                                                           <  6:15388/63‑1 (MQ=255)

                                                                                                                                                                                                               tAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCtt                                                                                                                                                                                           >  4:127080/1‑111 (MQ=255)

                                                                                                                                                                                                                aaCAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATc                                                                                                                                                    >  5:184951/1‑149 (MQ=255)

                                                                                                                                                                                                                    aGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAatat                                                                                                                              >  2:758636/1‑167 (MQ=255)

                                                                                                                                                                                                                     gTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCAt                                                                                                                                                     <  7:845876/143‑1 (MQ=255)

                                                                                                                                                                                                                      ttGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTc                                                                                                                                              >  7:1310478/1‑149 (MQ=255)

                                                                                                                                                                                                                       tGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTTGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCt                                                                                                                                             <  6:510654/149‑1 (MQ=255)

                                                                                                                                                                                                                         atatTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGt                                                                                                                                           <  8:1310478/149‑1 (MQ=255)

                                                                                                                                                                                                                          tatTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatc                                                                                                                                                                                      <  6:937566/105‑1 (MQ=255)

                                                                                                                                                                                                                            ttAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                                                                                                                     >  5:286673/1‑104 (MQ=255)

                                                                                                                                                                                                                            ttAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                                                                                                                     <  6:286673/104‑1 (MQ=255)

                                                                                                                                                                                                                            ttAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACg                                                                                                                                        >  7:1505975/1‑149 (MQ=255)

                                                                                                                                                                                                                              aGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                                                                                                                           >  5:443539/1‑128 (MQ=255)

                                                                                                                                                                                                                               gATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGgag                                                                                                                                     <  6:839408/149‑1 (MQ=255)

                                                                                                                                                                                                                                aTTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGgaga                                                                                                                                    >  5:572685/1‑149 (MQ=255)

                                                                                                                                                                                                                                    ttATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGa                      >  2:230458/1‑255 (MQ=255)

                                                                                                                                                                                                                                    ttATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGgta  >  2:90011/1‑275 (MQ=255)

                                                                                                                                                                                                                                       tttttAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATc                                                                                                                             <  6:711355/149‑1 (MQ=255)

                                                                                                                                                                                                                                       tttttAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATc                                                                                                                             <  6:791389/149‑1 (MQ=255)

                                                                                                                                                                                                                                       tttttAAGTTACGGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATc                                                                                                                             >  8:84843/1‑149 (MQ=255)

                                                                                                                                                                                                                                          ttAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTc                                                                                                                          <  7:478969/149‑1 (MQ=255)

                                                                                                                                                                                                                                               ttACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGtt                                                                                                                     >  8:519048/1‑149 (MQ=255)

                                                                                                                                                                                                                                               ttACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGtt                                                                                                                     >  7:1052273/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   tGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGc                                                                                                                 <  6:1661926/149‑1 (MQ=255)

                                                                                                                                                                                                                                                        ac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCaa                                                                                                                        <  8:931616/137‑1 (MQ=255)

                                                                                                                                                                                                                                                        ac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCa                                                                                                             >  8:1630601/1‑148 (MQ=255)

                                                                                                                                                                                                                                                        ac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCa                                                                                                             <  7:1630601/148‑1 (MQ=255)

                                                                                                                                                                                                                                                         c‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGt                                                                                                           >  5:107163/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          aaataaaaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGt                                                                                                           >  7:690904/5‑149 (MQ=255)

                                                                                                                                                                                                                                                           aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                                                                                                          >  6:1450371/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 aaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGa                                                                                                    <  7:1542129/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                 aaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCg                                                                                                                  >  7:1067175/1‑135 (MQ=255)

                                                                                                                                                                                                                                                                   aGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc                                                                                                  >  8:1118590/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   aGGCACGTCAGATGACGTGCCTTTTTTATCTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                                                                                                          >  7:943828/1‑141 (MQ=255)

                                                                                                                                                                                                                                                                     gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGGGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc                                                                                                >  7:1244516/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                             tCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGa                                                                                        >  8:380443/1‑149 (MQ=255)

ATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTA  >  NC_000913/4366586‑4367088

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