BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,388,669	C→A	5.5%	noncoding (99/130 nt)	PP_mr31 ←	rmf regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,388,669	0	C	A	5.5%	335.0 / 4.7	104	noncoding (99/130 nt)	PP_mr31	rmf regulatory region
Reads supporting (aligned to +/- strand): ref base C (53/45); new base A (3/3); total (56/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

GGAAGATTGAAGGGGCAAAGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTCATAAGCACGTAACCT  >  NC_002947/2388522‑2388816

ggaagATTGAAGGGGCAAAGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTGCg                                                                                                                                                    >  8:212280/1‑149 (MQ=255)

 gaagATTGAAGGGGCAAAGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTcc                                                                                                                                                     >  4:279392/1‑147 (MQ=255)

 gaagATTGAAGGGGCAAAGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTcc                                                                                                                                                     >  5:204476/1‑147 (MQ=255)

  aagATTGAAGGGGCAAAGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAg                                                                                                                                                  <  3:120509/149‑1 (MQ=255)

         aaGGGGCAAAGTTCGCGGGATTTGCCGGTGACCCCGTTTTGGTAGCCACGGGAATATGCTCTTTCGACCGGCTCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCGGAGAGGCCgac                                                                                                                                            <  3:68351/148‑3 (MQ=255)

         aaGGGGCAAAGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCTACGGATCAC‑GCTTATGTCTTCTCATTGGGTGTTGCCCTCACTTGTTGACTGTAATGTCGCGTCGGCCTCTCCGAGGCCggg                                                                                                                                            <  6:337172/148‑1 (MQ=255)

         aaGGGGCAAAGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCggg                                                                                                                                            >  8:196159/1‑148 (MQ=255)

         aaGGGGCAAAGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCggg                                                                                                                                            <  3:279392/148‑1 (MQ=255)

         aaGGGGCAAAGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCggg                                                                                                                                            >  7:331206/1‑148 (MQ=255)

         aaGGCGCAAAGTTCGCGGGATTTGCCGGTCACCCCGTATTGTCAGCCACGTTAATATGCTCTTTCCAACGCATCAC‑GCTTAAGTCTTCTCATCGGGTGTTGCCCTCACTTGTTCACTGTAATGCCCCGGCGGCCTCGGAGAGGCCgacg                                                                                                                                           <  4:284732/149‑4 (MQ=255)

            gggCAAAGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCCGGCag                                                                                                                                         >  2:133245/1‑148 (MQ=255)

              gCAAAGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGATTACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGt                                                                                                                                      >  1:302336/1‑149 (MQ=255)

                aaaGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCg                                                                                                                                    >  3:190914/1‑149 (MQ=255)

                aaaGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGCCGGCCTCGGAGAGGCCgacggg                                                                                                                                         >  1:213970/1‑139 (MQ=255)

                aaaGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGCCGGCCTCGGAGAGGCCgacggg                                                                                                                                         <  2:213970/144‑6 (MQ=255)

                    ttCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCACGTCGGCCACTCCGAGGCCGGGCACAGTCGTTCt                                                                                                                                >  4:18297/1‑149 (MQ=255)

                            aTTTGCCGGTGACCCCGTATTGGTAGCCACTTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGttt                                                                                                                        <  3:28233/149‑1 (MQ=255)

                               tGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTgc                                                                                                                      <  4:190914/148‑1 (MQ=255)

                                   ggTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGt                                                                                                                          >  4:256437/1‑140 (MQ=255)

                                   ggTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGt                                                                                                                          <  3:256437/140‑1 (MQ=255)

                                         cccGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCg                                                                                                           >  3:80146/1‑149 (MQ=255)

                                         cccGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCg                                                                                                           >  3:189358/1‑149 (MQ=255)

                                          ccGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTc                                                                                                                                     <  6:83419/122‑1 (MQ=255)

                                          ccGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTAGCTGCCgg                                                                                                          >  5:25869/1‑149 (MQ=255)

                                           cGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTgcgc                                                                                                                    <  6:263158/138‑1 (MQ=255)

                                           cGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTgcgc                                                                                                                    >  5:263158/1‑138 (MQ=255)

                                            gTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCg                                                                                                        >  8:192021/1‑149 (MQ=255)

                                                 ggTAGCCACGTGAATATGCTCTTTCCAACGGATCACTGCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGc                                                                                                    >  1:59170/1‑149 (MQ=255)

                                                 ggTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGc                                                                                                    <  6:204476/148‑1 (MQ=255)

                                                      ccACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCt                                                                                               >  2:225852/1‑148 (MQ=255)

                                                          gTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTTAAAg                                                                                          >  2:295989/1‑149 (MQ=255)

                                                           tGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGt                                                                                         <  4:80146/149‑1 (MQ=255)

                                                              atatGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCt                                                                                               >  5:263857/1‑140 (MQ=255)

                                                              atatGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCt                                                                                               <  6:263857/140‑1 (MQ=255)

                                                               tatGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCGCTACCAGGCCGGGGAGAGTCGTTGTGGCGGGTTTCGCGAGCGTCCGGCGGGGCGAGCAGaaaattccc                                                                                     >  6:235253/1‑143 (MQ=255)

                                                                 tGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTgccgcc                                                                                   >  5:204438/1‑149 (MQ=255)

                                                                  gCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCt                                                                                  >  1:133257/1‑149 (MQ=255)

                                                                   ctctTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTc                                                                                 >  2:49352/1‑149 (MQ=255)

                                                                   ctctTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTc                                                                                 >  7:176626/1‑149 (MQ=255)

                                                                    tctTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCg                                                                                >  2:250281/1‑149 (MQ=255)

                                                                         ccAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGc                                                                            >  1:136673/1‑148 (MQ=255)

                                                                         ccAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGc                                                                            <  1:225852/148‑1 (MQ=255)

                                                                          cAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCg                                                                           >  6:68080/1‑148 (MQ=255)

                                                                          cAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCg                                                                           >  6:283684/1‑148 (MQ=255)

                                                                          cAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCg                                                                           >  6:261658/1‑148 (MQ=255)

                                                                          cAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCACTCCGAGGCCGGGCAGAGTGGTTGTGCCGGTTTTCGCTCGCTGCCCGCGTGGGGGGCTGAAAGTGCCGCCTCGTTggg                                                                           >  6:35590/1‑146 (MQ=255)

                                                                           aaCGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGGCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCg                                                                           >  1:295202/1‑147 (MQ=255)

                                                                           aaCGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCg                                                                           >  3:122237/1‑147 (MQ=255)

                                                                            aCGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCGCGTTGcata                                                                         >  2:161632/1‑145 (MQ=255)

                                                                             cGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGc                                                                       <  2:59170/149‑1 (MQ=255)

                                                                              ggATCAC‑GCTTAAGTCTTCTCCTAGGCGGGTGCCCTCACTTGTTGACTGTACTGTCCCGTCGGCCTCTCCGCGGCCGGGCAGAGTCGTTCTACCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCg                                                                      <  8:124338/149‑1 (MQ=255)

                                                                                aTCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGc                                                                    <  6:282844/149‑1 (MQ=255)

                                                                                 tCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCTTCGCAACGAggcggca                                                                   <  3:171806/149‑2 (MQ=255)

                                                                                  cAC‑GCTTAAGGCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGCCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGGTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAACTGCCGCTTCGCAACGAggcggcac                                                                  <  7:55106/149‑3 (MQ=255)

                                                                                   aC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGa                                                                          <  2:224661/140‑1 (MQ=255)

                                                                                   aC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGa                                                                          >  1:224661/1‑140 (MQ=255)

                                                                                    c‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTga                                                                <  2:136673/149‑1 (MQ=255)

                                                                                      gCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTgag                                                               <  2:295202/149‑1 (MQ=255)

                                                                                       cTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTgaga                                                              >  7:142422/1‑149 (MQ=255)

                                                                                        ttAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTGTCCGAGGCCGGGCAGAGTCGTAATGCGGCTTTTCCCTACCTGCCGGCGTGGCGAGCTGAAACTGCGGCCTCGTTGGGAAGCGGCCTga                                                                >  8:105923/1‑146 (MQ=255)

                                                                                        ttAAGTCATCTCATAGGGTGTTGCCCTGACTTGTTGACTGTAATCTGCCGTCGGCCTCTCCCAGGCCGGGCAGAGTCGTTCTGCCGGATTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTCCCGCCTCGTTGCGAAGCGGCCTga                                                                <  7:105923/146‑1 (MQ=255)

                                                                                            gTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGt                                                         <  5:235253/149‑1 (MQ=255)

                                                                                             tcttctCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGtt                                                        <  1:250281/149‑1 (MQ=255)

                                                                                                tctcATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCt                                                      <  4:122237/148‑1 (MQ=255)

                                                                                                  tcATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGtt                                                        <  3:325947/144‑1 (MQ=255)

                                                                                                  tcATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGtt                                                        >  4:325947/1‑144 (MQ=255)

                                                                                                  tcATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATc                                                   >  1:127516/1‑149 (MQ=255)

                                                                                                  tcATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGATGCGAAGCGGCCTGAGACCAGTTCTATc                                                   >  7:169436/1‑149 (MQ=255)

                                                                                                    aTAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATc                                                   <  7:192021/147‑1 (MQ=255)

                                                                                                    aTAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCg                                                  >  5:132708/1‑148 (MQ=255)

                                                                                                     tAGGGTGTTGCCCTCTCTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCCCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGaa                                                <  3:18297/149‑1 (MQ=255)

                                                                                                          tgtTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGa                                           <  8:142422/149‑1 (MQ=255)

                                                                                                            ttGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGt                                         <  2:133257/149‑1 (MQ=255)

                                                                                                                ccTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACg                                     <  4:189358/149‑1 (MQ=255)

                                                                                                                 cTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACgg                                    <  7:14760/149‑1 (MQ=255)

                                                                                                                 cTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACgg                                    >  8:14760/1‑149 (MQ=255)

                                                                                                                    aCTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGaa                                                                         >  2:243490/1‑109 (MQ=255)

                                                                                                                    aCTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGaa                                                                         <  1:243490/109‑1 (MQ=255)

                                                                                                                      ttgttgACTGTAATGTCCCGTCGGCCGCTCCGCGGCCCGGCAGAGTGGTTGTGGCGGGTTTCGCGAGCGGCCCGCGGGGGGAGCGGAAAATTCCCCCGCGGTTCGAAGCGGCCGGAGACCACGTCTATAGAATGCGAGTCACCGtgtg                                >  6:250138/1‑148 (MQ=255)

                                                                                                                        gttgACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGa                               <  6:204438/147‑1 (MQ=255)

                                                                                                                        gttgACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAAt                             >  7:192912/1‑149 (MQ=255)

                                                                                                                         ttgACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAAt                             <  7:12834/148‑1 (MQ=255)

                                                                                                                         ttgACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAAt                             <  5:261658/148‑1 (MQ=255)

                                                                                                                           gACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGGGCGCTGAAAGTGCCGCCGCGTTGCGCAGCGGCCTGCGACCAGtt                                                        >  6:225740/1‑119 (MQ=255)

                                                                                                                           gACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGa                           <  5:250138/148‑1 (MQ=255)

                                                                                                                             cTGTAATGTCCCGTCGGCCTCTCCGAGGCCGCGCAGAGTCGTTCTTCCGGTTTGCCCTCGCTGACCGCGTGGCGACCTGAAAGTGCCGCCACGTTGCGAAGCGGCCGGAGACCAGTTCGATAGAATGCGAGTGCCGCTGTGAATGAcc                         >  4:281308/1‑148 (MQ=255)

                                                                                                                                     tCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTc                 <  8:176626/148‑1 (MQ=255)

                                                                                                                                     tCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTc                 >  8:17371/1‑148 (MQ=255)

                                                                                                                                     tCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTCa                >  5:201167/1‑149 (MQ=255)

                                                                                                                                       ccGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTCAt               >  1:126746/1‑148 (MQ=255)

                                                                                                                                       ccGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTCAt               >  4:153754/1‑148 (MQ=255)

                                                                                                                                         gTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTCATAAg            >  6:258168/1‑149 (MQ=255)

                                                                                                                                         cccGGCCTCGGAGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTCATAAg            >  7:219995/3‑149 (MQ=255)

                                                                                                                                          ccGGCCTCGGAGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTCATAAg            >  6:331557/2‑148 (MQ=255)

                                                                                                                                            ggCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTCATAAGCa          <  5:35590/148‑1 (MQ=255)

                                                                                                                                             gCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTg                                                                             >  5:260445/1‑80 (MQ=255)

                                                                                                                                             gCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTg                                                                             <  6:260445/80‑1 (MQ=255)

                                                                                                                                               ctctCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTCATAAGCACGt       >  2:130289/1‑148 (MQ=255)

                                                                                                                                                   ccGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTCATAAGCACGTAACCt  >  3:163793/1‑149 (MQ=255)

                                                                                                                                                    cGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTCATAAGCACGTAAc    <  1:130289/146‑1 (MQ=255)

GGAAGATTGAAGGGGCAAAGTTCGCGGGATTTGCCGGTGACCCCGTATTGGTAGCCACGTGAATATGCTCTTTCCAACGGATCAC‑GCTTAAGTCTTCTCATAGGGTGTTGCCCTCACTTGTTGACTGTAATGTCCCGTCGGCCTCTCCGAGGCCGGGCAGAGTCGTTCTGCCGGTTTGCGCTCGCTGCCGGCGTGGCGAGCTGAAAGTGCCGCCTCGTTGCGAAGCGGCCTGAGACCAGTTCTATCGAATGCGAGTCACGGTGTGAATGACCGATTTGTCATAAGCACGTAACCT  >  NC_002947/2388522‑2388816

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