BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	699,132	C→G	5.0%	E16Q (GAG→CAG)	asnB ←	asparagine synthetase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	699,132	0	C	G	5.0%	342.1 / 5.9	99	E16Q (GAG→CAG)	asnB	asparagine synthetase B
Reads supporting (aligned to +/- strand): ref base C (49/45); new base G (2/3); total (51/48)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

AGAGAGGTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTGAAAAAAT  >  NC_000913/698988‑699279

agagagGTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACt                                                                                                                                                 >  4:382390/1‑149 (MQ=255)

 gagagGTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACt                                                                                                                                                 <  7:1543/148‑1 (MQ=255)

  agagGTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCTTCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGc                                                                                                                                               <  8:205302/149‑1 (MQ=255)

     ggTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTc                                                                                                                                            <  4:231820/149‑1 (MQ=255)

     ggTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTc                                                                                                                                            <  6:144707/149‑1 (MQ=255)

      gTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGTTGAGCTGCGtaa                                                                                                                                           <  7:175869/149‑3 (MQ=255)

      gTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGTTGAGCTGCGtaa                                                                                                                                           <  2:142331/149‑3 (MQ=255)

      gTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCt                                                                                                                                           >  3:141548/1‑149 (MQ=255)

      gTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATAAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCt                                                                                                                                           >  4:315479/1‑149 (MQ=255)

       ttGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTg                                                                                                                                          >  4:117990/1‑149 (MQ=255)

           gcCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGTTGAGCTGCGtaa                                                                                                                                           <  3:268689/144‑3 (MQ=255)

           gcCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGTTGAGCTGCGtaa                                                                                                                                           >  4:268689/1‑142 (MQ=255)

           gcCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGGGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGtttt                                                                                                                                      >  6:271622/1‑149 (MQ=255)

            cccccGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTg                                                                                                                                     <  4:141548/149‑1 (MQ=255)

             ccccGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGa                                                                                                                                    <  3:152428/149‑1 (MQ=255)

             ccccGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCAGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGa                                                                                                                                    >  3:76144/1‑149 (MQ=255)

                 gcgTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATc                                                                                                                                <  4:45940/149‑1 (MQ=255)

                 gcgTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATc                                                                                                                                <  5:271622/149‑1 (MQ=255)

                 gcgTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGACACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATc                                                                                                                                <  5:125956/149‑1 (MQ=255)

                   gTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCg                                                                                                                               >  7:201403/1‑148 (MQ=255)

                   gTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCg                                                                                                                               >  6:255780/1‑148 (MQ=255)

                   gTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGAATCAGGCGTGACAGCTCGAGGGCGTTCTTACGCAGCTCACATGCGGCTGTTTTGATATCg                                                                                                                               >  6:251230/1‑148 (MQ=255)

                    ttAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGaa                                                                                                                             <  7:174961/149‑1 (MQ=255)

                    ttAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGaa                                                                                                                             >  3:289215/1‑149 (MQ=255)

                       aCGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATa                                                                                                                           <  5:410630/148‑1 (MQ=255)

                       aCGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATAc                                                                                                                          <  3:171045/149‑1 (MQ=255)

                       aCGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATAc                                                                                                                          >  1:333050/1‑149 (MQ=255)

                         gTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGc                                                                                                                        >  8:206731/1‑149 (MQ=255)

                         gTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGc                                                                                                                        >  8:87834/1‑149 (MQ=255)

                            aacaaTTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCa                                                                                                                      >  5:268992/1‑148 (MQ=255)

                            aacaaTTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCa                                                                                                                      >  8:163241/1‑148 (MQ=255)

                               aaTTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCaaaaa                                                                                                                  <  3:117990/149‑1 (MQ=255)

                                 ttGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAAtt                                                                                                                <  7:206731/149‑1 (MQ=255)

                                 ttGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAAtt                                                                                                                <  4:76144/149‑1 (MQ=255)

                                   gACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTg                                                                                                               >  4:95432/1‑148 (MQ=255)

                                    aCAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGaa                                                                                                             <  7:163241/149‑1 (MQ=255)

                                    aCAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGaa                                                                                                             >  5:389807/1‑149 (MQ=255)

                                    aCAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGa                                                                                                              <  1:93634/148‑1 (MQ=255)

                                       gACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAacac                                                                                                          <  8:292344/149‑1 (MQ=255)

                                              gTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCt                                                                                                   >  7:271439/1‑149 (MQ=255)

                                              gTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCt                                                                                                   >  3:203464/1‑149 (MQ=255)

                                                gggCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGa                                                                                                                                    >  1:40103/1‑114 (MQ=255)

                                                gggCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGa                                                                                                                                    <  2:40103/114‑1 (MQ=255)

                                                      gaATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAAc                                                                                                     <  3:166883/139‑1 (MQ=255)

                                                      gaATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAAc                                                                                                     >  4:166883/1‑139 (MQ=255)

                                                            cGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTg                                                                                     <  4:203464/149‑1 (MQ=255)

                                                              ttATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTgg                                                                                    >  4:104940/1‑148 (MQ=255)

                                                                aTCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTAt                                                                                 >  1:90370/1‑149 (MQ=255)

                                                                 tCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTAtt                                                                                <  6:268992/149‑1 (MQ=255)

                                                                   gCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATttg                                                                              >  5:172819/1‑149 (MQ=255)

                                                                    cTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATttgt                                                                             <  2:90370/149‑1 (MQ=255)

                                                                    cTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATttgt                                                                             <  3:382390/149‑1 (MQ=255)

                                                                    cTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATttgt                                                                             >  5:394536/1‑149 (MQ=255)

                                                                         aTAAATACCTGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCACACCTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAAAACATAACCTTCTCCGTTAACCTGGTATttgttaga                                                                         >  5:389578/1‑145 (MQ=255)

                                                                         aTAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATttgt                                                                             <  6:389578/144‑1 (MQ=255)

                                                                         aTAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCtt                                                                        >  8:167367/1‑149 (MQ=255)

                                                                         aTAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCtt                                                                        <  6:204078/149‑1 (MQ=255)

                                                                         aTAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCtt                                                                        <  5:255780/149‑1 (MQ=255)

                                                                         aTAAATACCGGACCAGTCCGGGCCACGATGACGCATAAGGCGTGACAGCTAGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCCTTAACCTGGTATTTGTTGCtt                                                                        >  1:385765/1‑149 (MQ=255)

                                                                            aaTACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAAc                                                                                                     <  7:300168/117‑1 (MQ=255)

                                                                            aaTACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAAc                                                                                                     >  8:300168/1‑117 (MQ=255)

                                                                              tACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCGTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTtgttgt                                                                   <  8:313733/149‑1 (MQ=255)

                                                                                ccGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTgtgt                                                                 <  7:87834/149‑1 (MQ=255)

                                                                                 cGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGtt                                                                >  2:73343/1‑149 (MQ=255)

                                                                                      cAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCtt                                                           >  8:107269/1‑149 (MQ=255)

                                                                                        gTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACCTAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGt                                                         <  2:162104/149‑1 (MQ=255)

                                                                                           cGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTATCCTGGTATTTGTTGCTTGTTGTGTTTGCTTGttt                                                       <  6:389807/148‑1 (MQ=255)

                                                                                               ccACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTaaaa                                                   >  4:306038/1‑148 (MQ=255)

                                                                                                 aCGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTaaaaaaa                                                <  6:102145/149‑1 (MQ=255)

                                                                                                   gATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATg                                              <  1:73343/149‑1 (MQ=255)

                                                                                                   gATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATg                                              <  4:289215/149‑1 (MQ=255)

                                                                                                    aTGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGc                                             <  3:95432/149‑1 (MQ=255)

                                                                                                              aGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAagca                                    <  1:303433/148‑1 (MQ=255)

                                                                                                               ggCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCttgtt                                                                     <  7:63333/114‑1 (MQ=255)

                                                                                                               ggCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCttgtt                                                                     >  8:63333/1‑114 (MQ=255)

                                                                                                               ggCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAagcag                                   >  4:284192/1‑148 (MQ=255)

                                                                                                               ggCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAagcag                                   >  3:371258/1‑148 (MQ=255)

                                                                                                                gCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAgcagca                                 <  3:315479/149‑1 (MQ=255)

                                                                                                                     aCAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGcc                                            <  1:126227/133‑1 (MQ=255)

                                                                                                                     aCAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATAGCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCCGCACtgtg                            >  5:257129/1‑149 (MQ=255)

                                                                                                                     aCAGCACGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTTAACACATAACCATCTCCGTTAACCAAGAATTTGTTGCTTGTTGTGTATGCATGTTTAAAAAAATgac                                            >  2:126227/1‑131 (MQ=255)

                                                                                                                       aGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAAAACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTAGTTGTGTTTGCTAGTTTAAAAAAATGCCGCAAAGCAGCCCTGTgcg                          >  4:208511/1‑149 (MQ=255)

                                                                                                                          tCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTgcgc                         >  6:183740/1‑147 (MQ=255)

                                                                                                                            gAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGc                     <  3:58791/149‑1 (MQ=255)

                                                                                                                            gAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGc                     >  6:111266/1‑149 (MQ=255)

                                                                                                                            gAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGc                     <  5:111266/149‑1 (MQ=255)

                                                                                                                            gAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACACCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGc                     >  4:65985/1‑149 (MQ=255)

                                                                                                                             aGGGCTTTCTTACGCAGCTCAACTGCGTCTGGTTTGATATCGAATACGCCAAAAATTGCACCCATCACCTGCTCCGGTAACCTGGGATTTGTTGCTTGGTGTGTTTGCTTGTTGAAAAAAATGCCGCACAGCAGCACTGTGCGCACGc                     >  5:161501/1‑148 (MQ=255)

                                                                                                                                  tttCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTg               >  4:127602/1‑149 (MQ=255)

                                                                                                                                    tCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGc             >  6:220518/1‑149 (MQ=255)

                                                                                                                                    tCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGc             >  6:372585/1‑149 (MQ=255)

                                                                                                                                     cTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGc             <  3:306038/148‑1 (MQ=255)

                                                                                                                                      ttACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCgg           >  1:104089/1‑149 (MQ=255)

                                                                                                                                      agaCGCAGTTGAGCTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCgg           >  4:165317/3‑149 (MQ=255)

                                                                                                                                         cGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTg         >  1:346847/1‑148 (MQ=255)

                                                                                                                                            aGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTgc                           >  5:185410/1‑127 (MQ=255)

                                                                                                                                            aGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTgc                           >  5:191524/1‑127 (MQ=255)

                                                                                                                                            aGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTgc                           <  6:185410/127‑1 (MQ=255)

                                                                                                                                            aGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTgc                           <  6:191524/127‑1 (MQ=255)

                                                                                                                                                cAACTGCGTCTGTTTTGATATCGAATACACCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTGAAAAAAt  >  7:3887/1‑148 (MQ=255)

AGAGAGGTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTGAAAAAAT  >  NC_000913/698988‑699279

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