BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,412,785	A→T	6.3%	G133G (GGA→GGT)	yjfM →	DUF1190 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,412,785	0	A	T	6.3%	649.2 / 2.1	220	G133G (GGA→GGT)	yjfM	DUF1190 family protein
Reads supporting (aligned to +/- strand): ref base A (75/120); new base T (12/2); total (98/122)
Fisher's exact test for biased strand distribution p-value = 1.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

CCGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACT  >  NC_000913/4412637‑4412933

ccGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGAt                                                                                                                                                     >  1:3076490/1‑150 (MQ=255)

 cGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGttt                                                                                                                                                    >  2:1318595/1‑150 (MQ=255)

    cATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGAttttt                                                                                                                                                 <  1:1949964/150‑1 (MQ=255)

      tttGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGt                                                                                                                                               <  2:1655256/150‑1 (MQ=255)

        tGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTa                                                                                                                                             <  2:36631/150‑1 (MQ=255)

         gcgcgcGTGGCTGGAACAACGCCAACACTGCTTTTTATGCCGATGTTCCGTAGATTATGACTCCGCATACCTGTCTGTCTAAGTACGATTATTGCTACTACGTCAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTaa                                                                                                                                            <  2:884339/150‑1 (MQ=255)

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

                 ggCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAATCCGGGtt                                                                                                                                    >  2:2427231/1‑150 (MQ=255)

                  gCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTTGAGTTTTTTTAATCCGGGttt                                                                                                                                   >  1:3687788/1‑149 (MQ=255)

                   cTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTATGGATTTTTTTTAATACGGGTTAt                                                                                                                                  >  2:2348994/1‑150 (MQ=255)

                    tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGTAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTAtt                                                                                                                                 <  1:3276342/150‑1 (MQ=255)

                    tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTAtt                                                                                                                                 <  1:3124720/150‑1 (MQ=255)

                       aacaacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGTTTTTTTTTTATCCGGGGTATTAGt                                                                                                                              >  1:910224/1‑150 (MQ=255)

                       aacaacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTGTTTATCCGGGTTTTTAGt                                                                                                                              >  1:556396/1‑150 (MQ=255)

                          aacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTaaa                                                                                                                           <  2:4100196/150‑1 (MQ=255)

                             gCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCGGGATTTTTTTTAATCCGGGTTATTCGTAAatat                                                                                                                        >  2:2726786/1‑150 (MQ=255)

                              ccAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGATTTTTTTTAATCCGGGTTATTAGTAAATATc                                                                                                                       >  2:4194422/1‑150 (MQ=255)

                              ccAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATc                                                                                                                       <  1:1030643/150‑1 (MQ=255)

                              ccAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATc                                                                                                                       <  1:1235213/150‑1 (MQ=255)

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

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

                               cAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGGGGGGGTTTTTGTTTATACGGGGTTTTCGTAAAGAtag                                                                                                                      >  2:1469805/1‑148 (MQ=255)

                                aaaaCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAAGCCGGGTTATTCGTAAAGATcgc                                                                                                                     >  2:3867082/1‑150 (MQ=255)

                                     ggCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTTGGTTTTTTTTATTGAGGGGTTTTAGTAAATATAGCGAttt                                                                                                                >  1:806826/1‑148 (MQ=255)

                                       cATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGc                                                                                                              <  1:109499/150‑1 (MQ=255)

                                         ttttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTTTTAAGCCGGGTTATTCTTAAAGATCGCGATGAGCCg                                                                                                            >  2:621206/1‑150 (MQ=255)

                                          tttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGTTTTTTGTTTATACGGGGTATTCGTAAAGATAGCGATGTGcccg                                                                                                           >  2:768431/1‑148 (MQ=255)

                                          tttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTATGGATTTTTTTTAAGCCCGGTTATTAGTAAATATAGAGATGAACCgt                                                                                                           >  2:2941275/1‑150 (MQ=255)

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

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

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

                                              tGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCGGGATTTTTTTTTATCCGGGTGTTTAGTAAAGATAGCGAAGAGCCggttgt                                                                                                       >  1:4069979/1‑150 (MQ=255)

                                               gCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGGGGGTTTTTTTTTAGACCGGTGATTGGTAAAGATCGCCATGAGACgtgtgtt                                                                                                      >  2:1208674/1‑150 (MQ=255)

                                                ccGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgtttgttt                                                                                                     <  1:1897837/150‑1 (MQ=255)

                                                 cGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAAGACGGGGTATTCGTAAACAAAGCCAAGAGCCGGTTTTTta                                                                                                    >  2:2323537/1‑150 (MQ=255)

                                                 cGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCAGGATTTTTGTTAAGACGGGTTATTCGTAAACATAGCGATCAACCCTTTGTTta                                                                                                    >  2:2478113/1‑150 (MQ=255)

                                                  gATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTGTAAGCCGGGGTATTCGTAAAGATCGCGATGAGCCGTTTGTTtat                                                                                                   >  2:2227941/1‑150 (MQ=255)

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

                                                          gaagaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAATCCGGGTTATTAGTAAAGATAGCGATGTGCCGTTTGTTTATAAAAgcgg                                                                                           >  2:3729893/1‑150 (MQ=255)

                                                          gaagaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGTTTTTTGTTTAGCCGGGGTAGTAGTAAAGATCGCGATGAGCCGTTTGTTTATAACAgcgg                                                                                           >  2:2145504/1‑150 (MQ=255)

                                                           aagaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAgcggc                                                                                          <  1:895221/150‑1 (MQ=255)

                                                                   gACTCAGCAGAACTGTCTGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCctc                                                                                  <  2:2564512/150‑1 (MQ=255)

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

                                                                   gACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGGGGATTTTTTTTTATACGGGGTATTAGTAAAGATAGGGATGTTCCGGTTGTTTATAAAAGcggcgggtacgc                                                                                  >  1:4248528/1‑144 (MQ=255)

                                                                      tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAAACCGGGGTAATCGTAAATATAGCGATGAGCCCGTTGTTTATAACAACCGCGGTTCCGCttt                                                                               >  2:1289517/1‑150 (MQ=255)

                                                                      tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTTAGCCGGGTTATTCGTAAAGATAGCGCTGAGCCGGTTGGTTTTAACAGCCGCCGTTGCTCttt                                                                               >  2:2655153/1‑150 (MQ=255)

                                                                      tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTAGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAGAAACGGTTATTGGTAAAGAAATCGAGGAGGCGATTTTTTTTAAAACAGGCGGttgcacgtt                                                                               >  1:1608942/1‑143 (MQ=255)

                                                                         gcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGGGGCATTTTGTTAATACGGGTGATTCGTAAATTTCGCGAAGAGCCCGTTTTTTATAAAAGCGGCGGTTACACATTTgg                                                                            >  2:523497/1‑149 (MQ=255)

                                                                          cagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGCTTTTTTTTTATCCGGGGTTTTCGGTAAGATAGCGGTGTGCCGGTTGTTTTTAAAAGCCGCGGGTGCTCATTTggt                                                                           >  1:632996/1‑148 (MQ=255)

                                                                           agAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGGTTTTTTTTTATCCGGGGTATTAGTAAAAGTCGCGATGAGCCGGTTGTTTTTTAAAGAGGCGGGTGCTCATTTGGtt                                                                          >  1:2688443/1‑150 (MQ=255)

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

                                                                                     gTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGGTTTTTGTTAATCCGGGGTATTCGTAAAGATCGCGATGAGCCGGTTTTTTATAAAAGCGGCGGTTCCTCTTTTTCTTGGCGCCCact                                                                >  2:3022666/1‑148 (MQ=255)

                                                                                         aaGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTgg                                                            <  1:1089697/150‑1 (MQ=255)

                                                                                          aGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgc                                                           <  1:4194422/150‑1 (MQ=255)

                                                                                             aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGATTTTTTTTAAGCCGGGTTATTAGTAAAGATCGCGATGAGCCGTTTTTTTATAACAGCGGCGGTTACACTTTTGCTTAGCGCCCAGTATGGCgca                                                        >  2:1784428/1‑150 (MQ=255)

                                                                                                aaaaTTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcagca                                                     <  2:3076490/150‑1 (MQ=255)

                                                                                                aaaaTTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcagca                                                     <  2:390743/150‑1 (MQ=255)

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

                                                                                                        tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                             <  1:2612232/150‑1 (MQ=255)

                                                                                                        tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                             <  1:461917/150‑1 (MQ=255)

                                                                                                        tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                             <  1:3986528/150‑1 (MQ=255)

                                                                                                             cGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGATTTTTTTTAAGCCGGGTTTTTAGTAAAGATCGCGATGAGCCGTTTTTTTATAAAACCGGCGGGTACACTTTTTCTTGGCGCCCAGTCATGCGCAGCACATAGGGTGATTa                                        >  2:2335782/1‑150 (MQ=255)

                                                                                                             cGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTTGTTTTTTGTTTATGCGGGGTTTTAGTTAAGTTGGCGATGTTGCGTTTGTTTTTTACAGGGGGGGTTACTATTTTTATTGGGGGCCAGTGTGGCGCAGCACTTGTGGTGATTa                                        >  2:1368705/1‑150 (MQ=255)

                                                                                                                cAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGTTTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTc                                     <  1:3883864/150‑1 (MQ=255)

                                                                                                                cAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTc                                     <  1:1770175/150‑1 (MQ=255)

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

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

                                                                                                                     ttGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAATCCCGGTTAATAGTAAAGATCGCGATGAGCCCTTTGTTTATAAAAGCCGCCGTTCCTCATTTTCATCGCGCCCACTCAGGCGCAGCACATCTGGTGATTACACCAGgc                                >  2:1773517/1‑150 (MQ=255)

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

                                                                                                                        aaCAGAGTTGGATCCCGGTCGTTTCTTGATTTTTTTTAATCCGGGTTATTCGTAAAGATCGCGATGAGCCGGTTTTTTATAAAACCGGCGGTTACTCTTTTTCTTCGCGCCCCCGATGGCGCAGCACATCGGGTGATTACACCAGgggcg                             >  2:811129/1‑149 (MQ=255)

                                                                                                                         aCAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTATGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATTTCAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTc                            <  1:44862/150‑1 (MQ=255)

                                                                                                                         aCAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTc                            <  1:79088/150‑1 (MQ=255)

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

                                                                                                                                 ggATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg                    <  1:1284251/150‑1 (MQ=255)

                                                                                                                                 ggATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg                    <  2:806826/150‑1 (MQ=255)

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

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

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

                                                                                                                                  gATCCCGGTCGTGTCTGGATTTTTGTTAAGCCGGGTTATTAGTAAAGATAGCCAAGATCCCTTTTTTTATAACAGCGGCGGTTACTCTTTTTCTTCGCGCCCAGTCTGGCGCAGCACATCTGGTGATTACTACAGGCGCTCCGGTTCTgg                   >  2:1425549/1‑150 (MQ=255)

                                                                                                                                   aTCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGc                  <  1:3791833/150‑1 (MQ=255)

                                                                                                                                       cGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaaaa              <  1:3639257/150‑1 (MQ=255)

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

                                                                                                                                           cGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAaga          <  2:3325870/150‑1 (MQ=255)

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

                                                                                                                                             tttCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGt        <  1:141321/150‑1 (MQ=255)

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

                                                                                                                                             tttCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGt        <  1:141235/150‑1 (MQ=255)

                                                                                                                                             tttCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAAAAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACACCTGGCGCTCAGGTTCTGGCAAAAAAGAGt        >  1:1836077/1‑150 (MQ=255)

                                                                                                                                              ttCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTc       <  2:3980353/150‑1 (MQ=255)

                                                                                                                                              gtCTGGATTTGTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTc       <  2:2955289/149‑1 (MQ=255)

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

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

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

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

                                                                                                                                                 tGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTa    <  2:446628/150‑1 (MQ=255)

                                                                                                                                                 tGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGGTCTGGCAAAAAAGAGTCTTa    >  1:1348606/1‑150 (MQ=255)

                                                                                                                                                  ggATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTAc   <  1:3305177/150‑1 (MQ=255)

                                                                                                                                                  ggATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTAc   <  1:1175544/150‑1 (MQ=255)

                                                                                                                                                   gATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTTCTTCGCGCCCAGTCTGGCGCAGCACGTCTGGTGATTAATCCAGGCGCTCCGGTTCTGGCAAAAAAGAGTATTAct  >  1:1876924/1‑150 (MQ=255)

                                                                                                                                                   gATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTAct  <  2:3437172/150‑1 (MQ=255)

CCGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACT  >  NC_000913/4412637‑4412933

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: `‑`