BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,374,901	A→T	5.7%	L11F (TTA→TTT)	dcuD →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,374,901	0	A	T	5.7%	359.3 / 8.1	117	L11F (TTA→TTT)	dcuD	putative transporter
Reads supporting (aligned to +/- strand): ref base A (59/51); new base T (3/4); total (62/55)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 134 total aligned reads displayed.

ACAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATA  >  NC_000913/3374755‑3375049

aCAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGTGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑attaat                                                                                                                                                    >  5:125712/1‑149 (MQ=255)

aCAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑attaat                                                                                                                                                    >  7:406860/1‑149 (MQ=255)

   aTTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATAAAACATCGTTGTTATCAATCTGCCGTTTATGGGCTTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTAc                                                                                                                                                 <  2:263285/149‑1 (MQ=255)

    ttAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGtaattaatac                                                                                                                                                  >  4:232668/1‑146 (MQ=255)

    ttAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTTATTAc                                                                                                                                                 >  4:40644/1‑148 (MQ=255)

    ttAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTAc                                                                                                                                                 >  7:127790/1‑148 (MQ=255)

       aaaaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATg                                                                                                                                             <  8:392736/149‑1 (MQ=255)

       aaaaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAAACTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATg                                                                                                                                             >  5:84520/1‑149 (MQ=255)

         aaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATggg                                                                                                                                           <  2:414946/149‑1 (MQ=255)

                  tCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGAt                                                                                                                                  >  8:62770/1‑149 (MQ=255)

                   ccGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATc                                                                                                                                 >  4:102296/1‑149 (MQ=255)

                             aaTAGCCTTTCAAAAAATATAACATAGTTGTTTTCAATCTGCCGTTTAGGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACt                                                                                                                       <  6:373547/149‑1 (MQ=255)

                                    ttGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAAcc                                                                                                                >  1:313720/1‑149 (MQ=255)

                                    ttGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACAGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAAcc                                                                                                                >  6:289103/1‑149 (MQ=255)

                                     tGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCt                                                                                                               >  8:268583/1‑149 (MQ=255)

                                        aaaaaaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAg                                                                                                            <  1:68232/149‑1 (MQ=255)

                                         aaaaaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGa                                                                                                                                               <  4:56726/113‑1 (MQ=255)

                                         aaaaaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAgg                                                                                                           >  5:106244/1‑149 (MQ=255)

                                          aaaaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCCTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCAACGT‑ATTAATTACGATGGGCAATTTGATCCTGAAAAACTACAAACCTCAggt                                                                                                          >  7:233367/1‑149 (MQ=255)

                                                  aCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctg                                                                                                   <  3:102296/148‑1 (MQ=255)

                                                  aCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgct                                                                                                    <  2:377742/147‑1 (MQ=255)

                                                    aTCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACaaa                                                                                                                  <  1:231119/131‑1 (MQ=255)

                                                    aTCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAggt                                                                                                          >  8:111017/1‑139 (MQ=255)

                                                    aTCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAggt                                                                                                          <  7:111017/139‑1 (MQ=255)

                                                       gttgttTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGc                                                                                              <  3:41741/148‑1 (MQ=255)

                                                        ttgttTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCg                                                                                            >  1:329001/1‑149 (MQ=255)

                                                           ttttCAATCTGCCGTGTCTGGGATTGACCGGTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAg                                                                                         <  6:125712/149‑1 (MQ=255)

                                                           tttgCAATCTGCCGTGTATGGGATTGACCGTTTTCTTTTGACCCGGAGATCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCa                                                                                          <  5:165830/148‑1 (MQ=255)

                                                            tttCAATCTGCCGTTTATTGTATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAg                                                                                         <  3:43303/148‑1 (MQ=255)

                                                                aaTCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATc                                                                                    >  1:107362/1‑149 (MQ=255)

                                                                      ccGTTTATGGGATTGACCGTTTTCTTTTGACACTGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTcc                                                                                <  1:277123/147‑1 (MQ=255)

                                                                      ccGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTcc                                                                                <  1:406423/147‑1 (MQ=255)

                                                                      ccGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTc                                                                                 <  8:32571/146‑1 (MQ=255)

                                                                      ccGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCtg                                                                              <  8:122411/149‑1 (MQ=255)

                                                                      ccGGTTATTGGATTGTCCGTCTTCTTTTGACACGGCGTTCAACAATGCTCGGAATAATTCTATCTGTCATCGT‑ATTAATTACTATGGGCTATTTGATACTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCtg                                                                              <  6:227154/149‑1 (MQ=255)

                                                                        gTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTc                                                                                 <  6:8233/144‑1 (MQ=255)

                                                                        gTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCtgat                                                                            >  1:393708/1‑149 (MQ=255)

                                                                           tATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGACCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCtgatga                                                                          <  5:91379/148‑1 (MQ=255)

                                                                            aTGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTgatgat                                                                         <  2:269208/148‑1 (MQ=255)

                                                                             tGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGAtgt                                                                       >  6:412158/1‑149 (MQ=255)

                                                                                gATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTAt                                                                                     <  1:46743/132‑1 (MQ=255)

                                                                                gATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTAt                                                                                     <  3:389171/132‑1 (MQ=255)

                                                                                gATTGACAGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTAt                                                                                     >  2:46743/1‑132 (MQ=255)

                                                                                 aTTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCgg                                                                   <  7:268583/149‑1 (MQ=255)

                                                                                 aTTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCgg                                                                   >  8:159552/1‑149 (MQ=255)

                                                                                  ttGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGgt                                                                  <  3:82040/149‑1 (MQ=255)

                                                                                    gACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggc                                                                                                 >  6:53046/1‑116 (MQ=255)

                                                                                     aCCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTTTCATCGT‑ATTAATTACGa                                                                                                                                               <  1:255153/69‑1 (MQ=255)

                                                                                     aCCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGgtgt                                                                <  6:84520/148‑1 (MQ=255)

                                                                                     aCCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTc                                                               >  7:271855/1‑149 (MQ=255)

                                                                                     aCCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTc                                                               >  2:60481/1‑149 (MQ=255)

                                                                                     aCCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTc                                                               >  2:124533/1‑149 (MQ=255)

                                                                                     aCCGTTTTCTTTTGACACGGAGTTCAACAATGTTAGGCATAATTATATCTGTCATCTT‑ATTAATTACGa                                                                                                                                               >  2:255153/1‑69 (MQ=255)

                                                                                      ccGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCt                                                              >  5:198658/1‑149 (MQ=255)

                                                                                                ttGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAggg                                                      <  5:289103/147‑1 (MQ=255)

                                                                                                  gACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCg                                                  >  3:50172/1‑149 (MQ=255)

                                                                                                   acacGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCgg                                                 >  2:302215/1‑149 (MQ=255)

                                                                                                     acGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑TTTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCggtg                                               <  5:89272/149‑1 (MQ=255)

                                                                                                     acGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCggtg                                               <  8:127790/149‑1 (MQ=255)

                                                                                                     acGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCggtg                                               <  6:106244/149‑1 (MQ=255)

                                                                                                      cGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCggtgg                                              >  4:19378/1‑149 (MQ=255)

                                                                                                       ggAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAgg                                                       >  8:425531/1‑139 (MQ=255)

                                                                                                       ggAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAgg                                                       <  7:425531/139‑1 (MQ=255)

                                                                                                        gAGTTCAACAATGTTCGGCATAATTATATCTGTCATCTT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGGATCTTCCTTATGATGTGCGGTGTCTGGTTAGGGTTCggtgggg                                            >  6:57716/1‑147 (MQ=255)

                                                                                                           ttCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCg                                                                    >  4:185805/1‑122 (MQ=255)

                                                                                                           ttCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCg                                                                    <  3:185805/122‑1 (MQ=255)

                                                                                                           ttCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGCTGATGGCTGCCGCAGGTATCTTCATGATTATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTAc                                         >  4:67277/1‑149 (MQ=255)

                                                                                                             caacaaTGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACt                                        <  5:412158/148‑1 (MQ=255)

                                                                                                             caacaaTGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACt                                        <  2:230012/148‑1 (MQ=255)

                                                                                                                  ttGTTCTGCATAATTATATCTGTCATCGT‑ATTTATTTCGATGGGCTATTTGATGGTTTAAAACTACTAACCGAATGTGGTTCTTGCTGCCGCAGGTTTCTTCCTGATGATTTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATcc                                  <  3:67277/148‑1 (MQ=255)

                                                                                                                        ggCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTTGTTAGGGTTCGGTGGTGTACTCGATCacacca                             >  5:60419/1‑148 (MQ=255)

                                                                                                                         gCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGATAGGGTTCGGTGGTGTACTCGATCCCACCaa                            >  1:281366/1‑148 (MQ=255)

                                                                                                                          gaTAATTATATCTGTCATCGT‑ATTAATGTCGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCaa                            >  1:89665/2‑147 (MQ=255)

                                                                                                                               ttATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGc                                                                     >  3:201739/1‑101 (MQ=255)

                                                                                                                               ttATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAagcagc                     <  2:393708/149‑1 (MQ=255)

                                                                                                                                 atatCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACTAACCGCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTTggg                                                      >  3:11946/1‑114 (MQ=255)

                                                                                                                                 atatCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAggg                                                      <  4:11946/114‑1 (MQ=255)

                                                                                                                                  tatCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCgg                   <  7:159552/148‑1 (MQ=255)

                                                                                                                                     cTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTAc               >  8:377264/1‑149 (MQ=255)

                                                                                                                                      tGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACt              >  5:315927/1‑149 (MQ=255)

                                                                                                                                      tGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACt              >  4:156456/1‑149 (MQ=255)

                                                                                                                                      tGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACt              >  7:266492/1‑149 (MQ=255)

                                                                                                                                      tGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGAGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACt              >  3:47968/1‑149 (MQ=255)

                                                                                                                                      gcccATCGT‑AATTAATACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCccacc                              >  5:247034/4‑133 (MQ=255)

                                                                                                                                      gcccATCGT‑AATTAATACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCccacc                              <  6:247034/130‑1 (MQ=255)

                                                                                                                                       cccATCGT‑AATTAATACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACtt             >  4:40001/3‑149 (MQ=255)

                                                                                                                                        tCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCggtggt                                             >  8:191046/1‑116 (MQ=255)

                                                                                                                                           tCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGAtc         >  2:160950/1‑149 (MQ=255)

                                                                                                                                            cGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCt                                                                               <  6:87527/78‑1 (MQ=255)

                                                                                                                                            cGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCt                                                                               >  5:87527/1‑78 (MQ=255)

                                                                                                                                              taattaa‑tacGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGAtcgtc      <  3:40001/147‑1 (MQ=255)

                                                                                                                                              t‑ATTAATTACGATGGGCTATTTGATCCTGAAACACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGAtcgtc      >  3:106818/1‑149 (MQ=255)

                                                                                                                                              t‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGAtcgtc      <  6:315927/149‑1 (MQ=255)

                                                                                                                                               ttttaattaTGATGGGATATTTGATCCTGAATAACTACAAACCTGAGGTGGTGCTGGCTGCCGCATGTATCTTAGTGTTGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCACACCAAAAGCAGCGGCTACTTGAtcgtc      <  4:48844/147‑1 (MQ=255)

                                                                                                                                                attaattaCGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGgt                                                                  <  3:235410/88‑1 (MQ=255)

                                                                                                                                                attaattaCGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGgt                                                                  >  4:235410/1‑88 (MQ=255)

                                                                                                                                                 ttaattaCGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGa    >  3:162989/1‑149 (MQ=255)

                                                                                                                                                 ttaattaCGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGa    <  4:162989/149‑1 (MQ=255)

                                                                                                                                                   aattaCGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGata  >  4:213179/1‑149 (MQ=255)

                                                                                                                                                   aattaCGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGGGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGat   <  6:198658/148‑1 (MQ=255)

ACAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGT‑ATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATA  >  NC_000913/3374755‑3375049

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