BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,412,794	A→T	6.3%	L136F (TTA→TTT)	yjfM →	DUF1190 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,412,794	0	A	T	6.3%	841.5 / 5.6	262	L136F (TTA→TTT)	yjfM	DUF1190 family protein
Reads supporting (aligned to +/- strand): ref base A (82/159); new base T (18/3); total (100/162)
Fisher's exact test for biased strand distribution p-value = 4.43e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

GCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCG  >  NC_000913/4412646‑4412942

gcgcgcGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTaa                                                                                                                                                     <  2:80122/150‑1 (MQ=255)

  gcgcgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGc                                                                                                                                                   <  1:4925381/150‑1 (MQ=255)

  gcgcgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTTTTAAGc                                                                                                                                                   >  2:1170134/1‑150 (MQ=255)

    gcgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGTTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCg                                                                                                                                                 <  1:5253351/150‑1 (MQ=255)

    gcgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCg                                                                                                                                                 <  1:2059939/150‑1 (MQ=255)

     cgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCgg                                                                                                                                                <  2:2788833/150‑1 (MQ=255)

      gTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCggg                                                                                                                                               <  1:5622563/150‑1 (MQ=255)

      gTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGGGTGGGTTTTTGTTAAGCCggg                                                                                                                                               >  2:1610073/1‑150 (MQ=255)

        ggCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGtg                                                                                                                                             >  2:1759411/1‑149 (MQ=255)

        ggCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTGTTAATCCGGGtt                                                                                                                                             >  2:4577971/1‑150 (MQ=255)

         gCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGTTTTTTTTTAAGCCGGGTTa                                                                                                                                            >  2:2850728/1‑150 (MQ=255)

         gCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTGTTAAGGCggggtg                                                                                                                                            >  2:1259604/1‑147 (MQ=255)

          cTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGGTTTTTGTTAAGCCGGGGTAt                                                                                                                                           >  2:343896/1‑150 (MQ=255)

          cTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTGTTAATCCGGGTTAt                                                                                                                                           >  1:388188/1‑150 (MQ=255)

          cTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGGGTGGATTTTTGTTTATCCgggggat                                                                                                                                           >  1:2500027/1‑146 (MQ=255)

              aacaacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAAGCCCGGTTATTCGt                                                                                                                                       >  2:5179765/1‑150 (MQ=255)

              aacaacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGTTTTTTTTTAATCCCGGTGATTGGt                                                                                                                                       >  2:2472801/1‑150 (MQ=255)

               acaacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTa                                                                                                                                      <  1:4741244/150‑1 (MQ=255)

                caccGCCACAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCTGCAGACCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGTACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTaa                                                                                                                                     <  2:5387242/147‑1 (MQ=255)

                caacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTaa                                                                                                                                     <  1:2738003/150‑1 (MQ=255)

                    gCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGAGTTTTTTTTATACGGGGGATTAGTAAAAAt                                                                                                                                 >  1:4006498/1‑150 (MQ=255)

                      cAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATcg                                                                                                                               <  2:1274004/150‑1 (MQ=255)

                       aaaaCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATcgc                                                                                                                              <  1:4494271/150‑1 (MQ=255)

                           cGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAATCCGGGGTATTCGTAAAGATAGCGAtt                                                                                                                          >  2:3601881/1‑149 (MQ=255)

                               aTTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTTTTAAGCCCGGTTTTTAGTAAAGAACGCGATGAGcc                                                                                                                      >  2:2669776/1‑150 (MQ=255)

                               aTTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCGGGGTTTTTTTTTATACGGGGTATTAGTAAAGAAAGCGATGAGcc                                                                                                                      >  2:714172/1‑150 (MQ=255)

                                 tttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgt                                                                                                                    <  1:1985270/150‑1 (MQ=255)

                                  ttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTTTTCGTAAAGATCGCGATGAGCCgtt                                                                                                                   <  1:5692585/150‑1 (MQ=255)

                                     tGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgtttgt                                                                                                                <  2:5821714/150‑1 (MQ=255)

                                      gCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGGCTGGATTTTTGTTAAAACGGGGTGTTAGTAAAGAAAGCCCTGAGCCgagtgtt                                                                                                               >  1:5528165/1‑150 (MQ=255)

                                        cGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGATTTTTTTTAATCCGGGTTATTAGTAAAGATAGCGATGAGCCGTTTTTTta                                                                                                             >  2:892007/1‑150 (MQ=255)

                                         gATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTtat                                                                                                            <  1:4828991/150‑1 (MQ=255)

                                         gATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTtat                                                                                                            <  1:2010267/150‑1 (MQ=255)

                                           tGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATaa                                                                                                          <  2:3275095/150‑1 (MQ=255)

                                           tGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATaa                                                                                                          <  2:1784374/150‑1 (MQ=255)

                                           tGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATaa                                                                                                          <  2:1005764/150‑1 (MQ=255)

                                           tGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATaa                                                                                                          <  1:1920897/150‑1 (MQ=255)

                                             ttCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAATCCGGGTTATTAGTAAATATAGCGAAGAGCCGTTTGTTTATAACa                                                                                                        >  2:2348581/1‑150 (MQ=255)

                                                cGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTTTTTATACGGGTTATTAGTAAAGATAGCGAAGATCCGTTTTTTtattaaagac                                                                                                     >  1:73667/1‑143 (MQ=255)

                                                 gaagaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTGTTAAGCCGGGGTATTAGTTAAGATCGCGATGTTCCGGTTTTTTTTaaaaaccg                                                                                                    >  2:3974940/1‑144 (MQ=255)

                                                 gaagaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGTTTTTTTTTTAGGCCGGGTATTAGGAAAGATCGCGATGATCCGTTTGTTTTTAAAAacgg                                                                                                    >  2:624605/1‑150 (MQ=255)

                                                  aagaaTATGACTCAGCAGAACTGTCATTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAgcggc                                                                                                   <  1:2207046/150‑1 (MQ=255)

                                                  aagaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAATTTGCTACTACGACACTGTTGATCAGAGTTGGTTCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAgcggc                                                                                                   <  2:189008/150‑1 (MQ=255)

                                                     aaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGt                                                                                                <  2:3516137/150‑1 (MQ=255)

                                                     aaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGt                                                                                                <  1:3045735/150‑1 (MQ=255)

                                                      atatGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGtt                                                                                               <  2:2025231/150‑1 (MQ=255)

                                                       tatGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACATTGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTTTTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTc                                                                                              <  1:38495/150‑1 (MQ=255)

                                                          gACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCctc                                                                                           <  2:1983387/150‑1 (MQ=255)

                                                             tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAATACGGGGTATTAGTAAATATCGCGATGTGCCGTTTGTTTATAAAAGCGGCGGGTCATCttt                                                                                        >  2:4081030/1‑150 (MQ=255)

                                                             tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTATGGAATTTTGTTTAGCCGGGTTATTAGTAAAGATCGCCATGTTCCGTTTTTTTATAAAACCGGCGGGTCCTCttt                                                                                        >  2:3131299/1‑150 (MQ=255)

                                                                gcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGGGGGGGTTTTTTTTTATCAGGGTGATTAGTAAAGATAGCGATGAGCCGGTTGTTTTTAAAACCGGCGGGTCCACgttttg                                                                                     >  1:4925176/1‑148 (MQ=255)

                                                                  agAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTTAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCtt                                                                                   <  2:2068869/150‑1 (MQ=255)

                                                                   gAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTc                                                                                  <  2:3096109/150‑1 (MQ=255)

                                                                    aaCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcg                                                                                 <  1:3196763/150‑1 (MQ=255)

                                                                    aaCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcg                                                                                 <  1:2454303/150‑1 (MQ=255)

                                                                          cAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCa                                                                           <  1:5662112/150‑1 (MQ=255)

                                                                            gTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGAGGGATTTTTTTAAGCCCGGTGATTCGTAAAGATCGCGATGAGCCCGTTGTTTATAAAAACGGCGGGTCCTCTGTTGAGTCGCGcccccg                                                                         >  2:1526790/1‑147 (MQ=255)

                                                                               tAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGGTTTTTTTTTTGACGGGGTTTTAGTAAAGATCGCGATGAGCCGGTTGTTTTTTAAAGCGGCGGGTGCTCATTTTCTTAGCGgccccacag                                                                      >  1:4228300/1‑144 (MQ=255)

                                                                                 aGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGATCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCGTCGCGCCCAGTCTGgc                                                                    >  2:4235088/1‑150 (MQ=255)

                                                                                 aGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATAGCGATGAGCCGTTTGTTTATAAAAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgc                                                                    >  2:4702653/1‑150 (MQ=255)

                                                                                  gTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgcg                                                                   <  1:2661078/150‑1 (MQ=255)

                                                                                   tACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGAGGTTATTTTTTAAGCCGGGTTATTAGTAAATATGGCGATGATCCGTTTTTTTATAAAAGAGGCGGTTACACTTTTGCGTCGCGCCCCATCTGgcgc                                                                  >  2:3686197/1‑150 (MQ=255)

                                                                                    aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTTAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca                                                                 <  1:1652890/150‑1 (MQ=255)

                                                                                    aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca                                                                 <  1:2668015/150‑1 (MQ=255)

                                                                                    aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca                                                                 <  1:3584386/150‑1 (MQ=255)

                                                                                    aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca                                                                 <  1:892007/150‑1 (MQ=255)

                                                                                      gAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGTTTTTTGTTAAGCCCGGTTATTGGTAAAGATAGCGATGAGCCGGTTGTTTATAAAAGCGGCGGTTCCACTTTTTCTTCGCGCCCAGTCTGgcgcccc                                                               >  2:5023600/1‑147 (MQ=255)

                                                                                              ctactaCGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGAATTTTTTTTAGCCGGGTTATTGGTAAAGATAGCGGTGATCCGTTTGTTTATAAAAGAGGCGGGTCCTCTTTTGGTTCGCGCCCCCGATGGCGCAGCACATCTgg                                                       >  2:760546/1‑150 (MQ=255)

                                                                                              ctactaCGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTgg                                                       <  1:1531626/150‑1 (MQ=255)

                                                                                                 ctacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGa                                                    <  1:2924075/150‑1 (MQ=255)

                                                                                                  tacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTACTGGTGAt                                                   <  2:3860631/150‑1 (MQ=255)

                                                                                                  tacGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGAATTTTTTTAATCCGGGTTATTCGTAAAGATCCCGATGAGCCCTTTTTTTATAACACCCGCGGTTACTCTTTTGCTTATCGCCCACTCAGGCCCAGCAAATATGGTGAt                                                   >  2:258472/1‑150 (MQ=255)

                                                                                                    cGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTa                                                 <  2:4635587/150‑1 (MQ=255)

                                                                                                     gACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTTGCTTTTTTTTAAACCCGGTTATTCGTAAAGAGCGCGATGAGCCGTTTTTTTATAACAGCGGCGGTTACGCATTTTCTTCGCGCCCAGTCTGGCGCAGCACATATggggttccc                                                >  2:2327652/1‑143 (MQ=255)

                                                                                                         aTGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGTGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCt                                            <  1:2260345/150‑1 (MQ=255)

                                                                                                         aTGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCt                                            <  2:1562706/150‑1 (MQ=255)

                                                                                                         aTGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCt                                            <  1:4245993/150‑1 (MQ=255)

                                                                                                         aTGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCt                                            <  1:708630/150‑1 (MQ=255)

                                                                                                         aTGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCt                                            <  1:647876/150‑1 (MQ=255)

                                                                                                          tGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTg                                           <  2:1730415/150‑1 (MQ=255)

                                                                                                             tGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGgcg                                        <  2:3025719/150‑1 (MQ=255)

                                                                                                                 cAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTcc                                    <  2:3184761/150‑1 (MQ=255)

                                                                                                                       tGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCt                              <  2:388188/150‑1 (MQ=255)

                                                                                                                         gATCCCGGTCGTTTCTTGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTgg                            <  1:3341082/150‑1 (MQ=255)

                                                                                                                         gATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTgg                            <  2:2219432/150‑1 (MQ=255)

                                                                                                                                gTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaaaaaa                     <  1:4404425/150‑1 (MQ=255)

                                                                                                                                    tttCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGt                 >  2:5324712/1‑150 (MQ=255)

                                                                                                                                      tCTGGATTTTTGTTTAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCt               <  1:1097695/150‑1 (MQ=255)

                                                                                                                                      tCTGGATTTTTGTTAAGCCGGGTTATTCGTAACGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCt               <  2:5480089/150‑1 (MQ=255)

                                                                                                                                      tCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGTGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCt               <  2:4373187/150‑1 (MQ=255)

                                                                                                                                      tCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCt               <  2:1446340/150‑1 (MQ=255)

                                                                                                                                        tGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCCGCACTTCTGGTGATTACTCCAGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTa             >  1:4708909/1‑150 (MQ=255)

                                                                                                                                        tGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTa             >  2:4903152/1‑150 (MQ=255)

                                                                                                                                        tGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTa             >  1:621637/1‑150 (MQ=255)

                                                                                                                                        tGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTa             >  1:2994690/1‑150 (MQ=255)

                                                                                                                                              tttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTc       <  1:2503813/150‑1 (MQ=255)

                                                                                                                                              tttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTc       >  1:1325515/1‑150 (MQ=255)

                                                                                                                                                 gTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCggg    >  1:3709199/1‑150 (MQ=255)

                                                                                                                                                  ttAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggc   <  2:3224386/150‑1 (MQ=255)

                                                                                                                                                  ttAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggc   <  2:3224812/150‑1 (MQ=255)

GCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCG  >  NC_000913/4412646‑4412942

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