BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,027,332:1	+C	5.7%	coding (594/993 nt)	ocd →	ornithine cyclodeaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,027,332	1	.	C	5.7%	435.3 / 18.1	105	coding (594/993 nt)	ocd	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand): ref base . (47/52); new base C (3/3); total (50/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 110 total aligned reads displayed.

AAAGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGATCACTACGT  >  NC_002947/5027183‑5027481

aaaGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCt                                                                                                                                                        >  2:2049112/1‑150 (MQ=255)

         cttcttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCa                                                                                                                                              <  1:2065057/150‑1 (MQ=255)

         cttctgGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCa                                                                                                                                              <  1:1542581/150‑1 (MQ=255)

          ttcttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑Ga                                                                                                                                                     >  1:2275714/1‑142 (MQ=255)

          ttcttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑Ga                                                                                                                                                     <  2:2275714/142‑1 (MQ=255)

          ttcttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAg                                                                                                                                             <  2:1598754/150‑1 (MQ=255)

          ttcttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGACACTGCCACAGAGCCGGTCGATGAAGCAAGGTAGAGCGCT‑GACATGGCAg                                                                                                                                             >  2:2323390/1‑150 (MQ=255)

             ttGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCa                                                                                                                                          >  2:1949329/1‑150 (MQ=255)

                           aaGCGAGTTCGGGTATGGGCTAGGCGTTCAATACAAGCACTGCAACGTGGCTTTCATCTCTCGAGCGTGTTCGGCTTTGGCCGCAATACCGGCAAATCGGGCGATGAAGCAATGTCGAGCGCT‑GACATTGACGTCAAATCTACGCCAAGc                                                                                                                            >  1:1151486/1‑150 (MQ=255)

                           aaGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGc                                                                                                                            >  2:1816610/1‑150 (MQ=255)

                            aGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCa                                                                                                                           >  1:104785/1‑150 (MQ=255)

                            aGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCa                                                                                                                           >  2:535072/1‑150 (MQ=255)

                                 gTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCaa                                                                                                                      <  1:1816610/150‑1 (MQ=255)

                                 gTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCaa                                                                                                                      >  1:757776/1‑150 (MQ=255)

                                 gTTCGGGTATGGGCAAGGAGCTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGGCGATGAAGAAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCaa                                                                                                                      >  2:450075/1‑150 (MQ=255)

                                  ttCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAAc                                                                                                                     <  1:808749/150‑1 (MQ=255)

                                     gggTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCaa                                                                                                                  >  1:757465/1‑150 (MQ=255)

                                       gTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTGGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAAtt                                                                                                                >  2:662504/1‑150 (MQ=255)

                                          tGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTa                                                                                                             <  2:104785/150‑1 (MQ=255)

                                           ggggTATGAGTTCCGCAAAAGCGCAGCAACGGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGGCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTAc                                                                                                            <  2:546406/150‑1 (MQ=255)

                                            ggCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACa                                                                                                           <  1:324834/150‑1 (MQ=255)

                                             gCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACAt                                                                                                          >  1:236514/1‑150 (MQ=255)

                                             gCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACAt                                                                                                          >  1:439058/1‑150 (MQ=255)

                                              cTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATc                                                                                                         >  1:959795/1‑150 (MQ=255)

                                                  gAGTTCCGCAAAAGCTCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAg                                                                                                     <  2:1978202/150‑1 (MQ=255)

                                                     ttCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCa                                                                                                  >  1:2266232/1‑150 (MQ=255)

                                                     ttCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCa                                                                                                  <  2:757465/150‑1 (MQ=255)

                                                         gCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCAc                                                                                                 <  2:612365/147‑1 (MQ=255)

                                                         gCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCAc                                                                                                 >  1:612365/1‑147 (MQ=255)

                                                         gCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTC‑AGCGCTCGACATTGct                                                                                                                                              >  2:401304/1‑101 (MQ=255)

                                                         gCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTC‑AGCGCTCGACATTGct                                                                                                                                              <  1:401304/102‑2 (MQ=255)

                                                          cAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGc                                                                                             <  2:439058/150‑1 (MQ=255)

                                                           aaaaGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCa                                                                                                      >  1:223638/1‑140 (MQ=255)

                                                           aaaaGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCa                                                                                                      <  2:223638/140‑1 (MQ=255)

                                                           aaaaGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCg                                                                                            <  2:1521391/150‑1 (MQ=255)

                                                             aaGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGc                                                                                          <  2:505780/150‑1 (MQ=255)

                                                              aGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGcc                                                                                         >  1:1766476/1‑150 (MQ=255)

                                                               gcgcAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGAAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGccc                                                                                        <  1:978582/150‑1 (MQ=255)

                                                                        cTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCat                                                                               >  1:338682/1‑150 (MQ=255)

                                                                          ggCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCatat                                                                             <  1:1661270/150‑1 (MQ=255)

                                                                              ttGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAAACTGATATCGCAGGCACTTGCGGCCCGGTCTCCATATCACg                                                                         >  1:737907/1‑150 (MQ=255)

                                                                                    cTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAgg                                                                                                    <  2:379938/117‑1 (MQ=255)

                                                                                    cTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAgg                                                                                                    >  1:379938/1‑117 (MQ=255)

                                                                                    cTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATg                                                                   >  1:1436700/1‑150 (MQ=255)

                                                                                      aTCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATggg                                                                 >  2:160135/1‑150 (MQ=255)

                                                                                          aGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCa                                                             >  2:755938/1‑150 (MQ=255)

                                                                                          aGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAAACAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCa                                                             <  2:1557906/150‑1 (MQ=255)

                                                                                          aGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCCGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCa                                                             >  2:203834/1‑150 (MQ=255)

                                                                                             ctctACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGac                                                          <  2:1766476/150‑1 (MQ=255)

                                                                                                tACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGc                                                                                          >  2:1880540/1‑115 (MQ=255)

                                                                                                tACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGc                                                                                          <  1:1880540/115‑1 (MQ=255)

                                                                                                tACGGCCTAGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACt                                                       <  2:1472274/150‑1 (MQ=255)

                                                                                                   ggCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGaa                                                    <  1:2049112/150‑1 (MQ=255)

                                                                                                   ggCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGaa                                                    <  1:252375/150‑1 (MQ=255)

                                                                                                    gCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAac                                                   >  2:682781/1‑150 (MQ=255)

                                                                                                      cTGGAAGCCACTGCCGCATAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAacac                                                 >  1:290861/1‑150 (MQ=255)

                                                                                                      cTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTg                                                      >  2:2390769/1‑145 (MQ=255)

                                                                                                      cTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTg                                                      <  1:2390769/145‑1 (MQ=255)

                                                                                                      cTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCGATGGGGGCAGACACTGAacac                                                 >  2:560141/1‑150 (MQ=255)

                                                                                                       tGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAAcaca                                                >  2:826295/1‑150 (MQ=255)

                                                                                                       tGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAAcaca                                                >  2:2217756/1‑150 (MQ=255)

                                                                                                               aCTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCACGCACTCTATCTCGTTTACAGCGCAGGCTCTTGTGGCCCGTTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAg                                        <  1:2098231/150‑1 (MQ=255)

                                                                                                               aCTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCa                                                                                                      >  2:781352/1‑88 (MQ=255)

                                                                                                               aCTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCa                                                                                                      <  1:781352/88‑1 (MQ=255)

                                                                                                               aCTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAg                                        >  1:2389833/1‑150 (MQ=255)

                                                                                                                cTGCCTCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAgg                                       <  1:2323390/150‑1 (MQ=255)

                                                                                                                 tGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAg                                        >  1:2165122/1‑148 (MQ=255)

                                                                                                                    cGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGc                                   <  1:511240/150‑1 (MQ=255)

                                                                                                                      cAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTc                                 <  1:826295/150‑1 (MQ=255)

                                                                                                                        gagTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCt                               <  1:160135/150‑1 (MQ=255)

                                                                                                                        gagTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCt                               <  1:1657183/150‑1 (MQ=255)

                                                                                                                          gtcggtcgATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGa                             <  2:236514/150‑1 (MQ=255)

                                                                                                                          gtcggtcgATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGa                             <  2:1436700/150‑1 (MQ=255)

                                                                                                                          gtcggtcgATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCCGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGa                             <  1:203834/150‑1 (MQ=255)

                                                                                                                           tcggtcgATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAg                            <  2:1151486/150‑1 (MQ=255)

                                                                                                                            cggtcgATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAgg                           <  2:608168/150‑1 (MQ=255)

                                                                                                                             ggtcgATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGt                          >  2:665624/1‑150 (MQ=255)

                                                                                                                              gtcgATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTa                         >  1:1276438/1‑150 (MQ=255)

                                                                                                                               tcgATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCCGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTaa                        <  1:1949329/150‑1 (MQ=255)

                                                                                                                               tcgATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTaa                        <  2:338682/150‑1 (MQ=255)

                                                                                                                                 gtgactgCAATGTC‑AGCGCTCGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAaca                                                  >  2:1830442/7‑122 (MQ=255)

                                                                                                                                 gtgactgCAATGTC‑AGCGCTCGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAaca                                                  <  1:1830442/116‑1 (MQ=255)

                                                                                                                                 gATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATa                      >  2:915692/1‑150 (MQ=255)

                                                                                                                                 gATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATa                      >  2:2348839/1‑150 (MQ=255)

                                                                                                                                    gAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATagca                   >  1:639051/1‑150 (MQ=255)

                                                                                                                                      aGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATagcagc                 <  1:2363616/150‑1 (MQ=255)

                                                                                                                                        cAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGcagcag               <  2:290861/150‑1 (MQ=255)

                                                                                                                                               gAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACgg                                                                        >  1:1604489/1‑86 (MQ=255)

                                                                                                                                               gAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACgg                                                                        <  2:1604489/86‑1 (MQ=255)

                                                                                                                                               gAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGATCa        <  2:1904790/150‑1 (MQ=255)

                                                                                                                                                aGCGCTCGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAg                            >  2:2361722/1‑130 (MQ=255)

                                                                                                                                                aGCGCTCGACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAg                            <  1:2361722/130‑1 (MQ=255)

                                                                                                                                                 gcgcT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGa                             <  2:250717/127‑1 (MQ=255)

                                                                                                                                                 gcgcT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGa                             >  1:250717/1‑127 (MQ=255)

                                                                                                                                                 gcgcT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGATCACt      <  1:535072/150‑1 (MQ=255)

                                                                                                                                                  cgcT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGATCACTa     <  2:2389833/150‑1 (MQ=255)

                                                                                                                                                     t‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAgg                                       <  2:2004705/113‑1 (MQ=255)

                                                                                                                                                     t‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAgg                                       >  1:2004705/1‑113 (MQ=255)

                                                                                                                                                     t‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGATCACTACGt  >  1:2252114/1‑150 (MQ=255)

                                                                                                                                                     t‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGATCACTACGt  >  1:814228/1‑150 (MQ=255)

AAAGCGCTCCTTCTTGTGCGGGATATCAAGCGAGTTCGGGTATGGGCTAGGAGTTCCGCAAAAGCGCAGCAACTGGCCTTGCAGCTATCGAGCCTCTACGGCCTGGAAGCCACTGCCGCAGAGTCGGTCGATGAAGCAATGTCGAGCGCT‑GACATTGCAGTCACCTGCACGCCAAGCAATCAACCAATTTTACATCGCAGGCACTTGCGGCCCGGTCTCCATATCACGGCTATGGGGGCAGACACTGAACACAAAAACGAGGTTGCTCCTGAGGTAATAGCAGCAGTCGATCACTACGT  >  NC_002947/5027183‑5027481

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