BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,385,206	T→A	5.5%	P456P (CCT→CCA)	yjeM →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,385,206	0	T	A	5.5%	491.9 / 9.8	140	P456P (CCT→CCA)	yjeM	putative transporter
Reads supporting (aligned to +/- strand): ref base T (67/65); new base A (3/5); total (70/70)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

TCCTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTA  >  NC_000913/4385058‑4385353

tCCTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCt                                                                                                                                                     <  4:65140/149‑1 (MQ=255)

  cTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCtgt                                                                                                                                                   <  1:458107/149‑1 (MQ=255)

    cgcgCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTgg                                                                                                                                                 <  2:329738/149‑1 (MQ=255)

      cgcTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGtt                                                                                                                                               <  3:105601/149‑1 (MQ=255)

      cgcTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGtt                                                                                                                                               <  4:528821/149‑1 (MQ=255)

         tGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGaa                                                                                                                                            <  4:843904/149‑1 (MQ=255)

         tGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGaa                                                                                                                                            >  3:592652/1‑149 (MQ=255)

         tGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGaa                                                                                                                                            >  3:737085/1‑149 (MQ=255)

         tGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCGCGACAGACGGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCACCCACAGGTTGaa                                                                                                                                            <  3:383587/149‑1 (MQ=255)

          ggCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaat                                                                                                                                           <  4:770891/149‑2 (MQ=255)

           gCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGc                                                                                                                                          >  3:673608/1‑149 (MQ=255)

             tttCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCg                                                                                                                                        >  4:173667/1‑149 (MQ=255)

             tttCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCg                                                                                                                                        >  4:691192/1‑149 (MQ=255)

               tCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGga                                                                                                                                      >  4:163908/1‑149 (MQ=255)

                  cGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGAc                                                                                                                                   >  3:823835/1‑149 (MQ=255)

                    ttCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTg                                                                                                                                 >  4:46564/1‑149 (MQ=255)

                    ttCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTg                                                                                                                                 <  4:466323/149‑1 (MQ=255)

                    ttCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTg                                                                                                                                 >  1:164243/1‑149 (MQ=255)

                    ttCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTg                                                                                                                                 >  1:111467/1‑149 (MQ=255)

                    ttCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTg                                                                                                                                 >  1:209502/1‑149 (MQ=255)

                     tCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTgg                                                                                                                                <  4:50059/149‑1 (MQ=255)

                       tttAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGa                                                                                                                              >  1:676068/1‑149 (MQ=255)

                       tttAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGa                                                                                                                              <  3:170233/149‑1 (MQ=255)

                       tttAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGa                                                                                                                              >  4:170233/1‑149 (MQ=255)

                         tAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACa                                                                                                                            >  2:690579/1‑149 (MQ=255)

                          aaaGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAg                                                                                                                           >  4:42336/1‑149 (MQ=255)

                          aaaGCACGTCAGGATCTCGACAGACCGGTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGCGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGAGTCAg                                                                                                                           >  4:396925/1‑149 (MQ=255)

                           aaGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGc                                                                                                                          >  3:364089/1‑149 (MQ=255)

                             gCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGgaga                                                                                                                                    <  3:866650/137‑1 (MQ=255)

                             gCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGcac                                                                                                                        >  3:372084/1‑149 (MQ=255)

                                cGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACAtt                                                                                                                     >  4:285554/1‑149 (MQ=255)

                                  tCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtgt                                                                                                                   <  2:209502/149‑1 (MQ=255)

                                  tCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtgt                                                                                                                   <  4:324858/149‑1 (MQ=255)

                                  tCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCAGTGGTTGAAGCCGGAGACTGGGACGGCACATtgt                                                                                                                   >  1:859842/1‑149 (MQ=255)

                                  tCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaat                                                                                                                                           >  2:172262/1‑124 (MQ=255)

                                   cAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAGTCATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAAAGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGAAAGCACATtgtg                                                                                                                  <  4:81893/149‑1 (MQ=255)

                                    aGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTgg                                                                                                                 <  4:372084/149‑1 (MQ=255)

                                    aGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTgg                                                                                                                 >  2:65027/1‑149 (MQ=255)

                                    aGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTgg                                                                                                                 <  2:673047/149‑1 (MQ=255)

                                     ggATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGga                                                                                                                <  2:527845/149‑1 (MQ=255)

                                      gATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGgat                                                                                                               >  3:716942/1‑149 (MQ=255)

                                       aTCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGgatg                                                                                                              >  3:727983/1‑149 (MQ=255)

                                        tctcGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGgatga                                                                                                             >  1:687455/1‑149 (MQ=255)

                                         ctcGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGgatgat                                                                                                            >  1:724477/1‑149 (MQ=255)

                                            gACAGACCGTTTGTGCTTTTTATACCGCATTTGTCGGCAATGATTGCGACTGTGGTTTTCGTACTGGTGGCGACCCTTGCGAACGTCTTCACCATCATTGAACCTGTGGTCGAAGCCGGCGACTGGGACAGCACATTGTGGATGATTgg                                                                                                         >  2:369005/1‑149 (MQ=255)

                                            gACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTgg                                                                                                         >  2:792026/1‑149 (MQ=255)

                                            gACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTgg                                                                                                         >  3:513940/1‑149 (MQ=255)

                                             aCAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGtt                                                                                                                                               <  1:585502/110‑1 (MQ=255)

                                             aCAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGtt                                                                                                                                               >  2:585502/1‑110 (MQ=255)

                                                 aCCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGcc                                                                                                    >  1:294485/1‑149 (MQ=255)

                                                 aCCGTTTGCGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGcc                                                                                                    >  3:403114/1‑149 (MQ=255)

                                                      ttGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtc                                                                                               <  1:651687/149‑1 (MQ=255)

                                                        gtgATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtctt                                                                                             <  3:390599/149‑1 (MQ=255)

                                                         tgATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtcttc                                                                                            <  2:592927/149‑1 (MQ=255)

                                                          gATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtcttct                                                                                           >  4:611616/1‑149 (MQ=255)

                                                                    aCGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGtt                                                                                 >  2:354949/1‑149 (MQ=255)

                                                                     cGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTa                                                                                <  4:755253/149‑1 (MQ=255)

                                                                       cATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGc                                                                              <  4:332417/149‑1 (MQ=255)

                                                                         tttGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGAtt                                                                                                           <  3:604008/118‑1 (MQ=255)

                                                                           tGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATg                                                                          <  2:883205/149‑1 (MQ=255)

                                                                                gCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAt                                                                     >  3:330245/1‑149 (MQ=255)

                                                                                gCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAt                                                                     >  1:201094/1‑149 (MQ=255)

                                                                                 cAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAtt                                                                    <  3:469418/149‑1 (MQ=255)

                                                                                    tGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCTGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTAc                                                                 <  4:727983/149‑1 (MQ=255)

                                                                                    tGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTAc                                                                 <  1:98878/149‑1 (MQ=255)

                                                                                      aTTGCGACAGTGTTTTGCGTACTGGTGGTCACATGTGCGAACGCCTTCACCATCTTTCAACCTGTGCTTAAAGACGGAGTCTGAGACAGCACATTGTGGATGATTTGCGGCCCGGTTTTCTTCTCGCTGTTAGCGATGGCGATTTACCa                                                               <  2:239418/149‑1 (MQ=255)

                                                                                          cGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAAc                                                           <  1:65027/149‑1 (MQ=255)

                                                                                           gACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATg                                                                          <  4:120728/133‑1 (MQ=255)

                                                                                           gACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATg                                                                          >  3:120728/1‑133 (MQ=255)

                                                                                             cAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTAt                                                        <  2:687455/149‑1 (MQ=255)

                                                                                               gTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTg                                                      <  4:513940/149‑1 (MQ=255)

                                                                                               gTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTg                                                      >  2:354007/1‑149 (MQ=255)

                                                                                                 ggTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCAGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgca                                                    <  2:859842/149‑1 (MQ=255)

                                                                                                  gTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCAACATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGaa                                                            <  1:263302/140‑1 (MQ=255)

                                                                                                      tCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCaga                                               <  2:22088/149‑1 (MQ=255)

                                                                                                               tggtgACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCaaaa                                      <  3:285554/149‑1 (MQ=255)

                                                                                                                  tgACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAAt                                   >  1:569935/1‑149 (MQ=255)

                                                                                                                       tttGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCa                              <  1:43093/149‑1 (MQ=255)

                                                                                                                         tGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGt                            <  4:276284/149‑1 (MQ=255)

                                                                                                                           cGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTgg                          <  2:799581/149‑1 (MQ=255)

                                                                                                                            gAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTggg                         >  3:251433/1‑149 (MQ=255)

                                                                                                                            gAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTggg                         >  1:868846/1‑149 (MQ=255)

                                                                                                                                 tCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATc                                  >  3:239091/1‑135 (MQ=255)

                                                                                                                                 tCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                    <  4:363719/149‑1 (MQ=255)

                                                                                                                                        catcatTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAAc             <  4:592652/149‑1 (MQ=255)

                                                                                                                                        catcatTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAAc             >  4:153066/1‑149 (MQ=255)

                                                                                                                                        catcatTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAAc             <  2:503257/149‑1 (MQ=255)

                                                                                                                                         atcatTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAAcc            >  4:594714/1‑149 (MQ=255)

                                                                                                                                         atcatTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAAcc            >  4:131370/1‑149 (MQ=255)

                                                                                                                                            atTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGa         <  4:251433/149‑1 (MQ=255)

                                                                                                                                             ttCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGAt        >  3:811926/1‑149 (MQ=255)

                                                                                                                                             ttCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGAt        >  2:412965/1‑149 (MQ=255)

                                                                                                                                             ttCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGAt        >  2:807455/1‑149 (MQ=255)

                                                                                                                                                aaCCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGcc     >  3:627810/1‑149 (MQ=255)

                                                                                                                                                 aCCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGccc    >  3:876663/1‑149 (MQ=255)

                                                                                                                                                  ccTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCt   <  4:354859/149‑1 (MQ=255)

                                                                                                                                                   cTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTa  <  3:173667/149‑1 (MQ=255)

                                                                                                                                                   cTGTGGTTGAAGCCGGAGACTGGACCAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCACTGTTAGCGATGCCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTa  >  4:410482/1‑149 (MQ=255)

                                                                                                                                                    tgtgGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGaa                 >  2:158046/1‑133 (MQ=255)

                                                                                                                                                    tgtgGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGaa                 <  1:158046/133‑1 (MQ=255)

TCCTCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTA  >  NC_000913/4385058‑4385353

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