Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,938,823	C→G	7.4%	intergenic (‑71/‑1726)	ddlA ← / → PP_4348	D‑alanine‑‑D‑alanine ligase A/cystathionine beta‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,938,823	0	C	G	7.4%	196.3 / 4.5	65	intergenic (‑71/‑1726)	ddlA/PP_4348	D‑alanine‑‑D‑alanine ligase A/cystathionine beta‑lyase
Reads supporting (aligned to +/- strand):  ref base C (33/27);  new base G (3/2);  total (36/29)
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

GTAGCGGCGGAGCGTAGCGCCACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATCGCGGTCTCTGCGCCGGAATGCATCAGCGGCCACAAA  >  NC_002947/4938676‑4938971                                                                                                                                                      |                                                                                                                                                     gTAGCGGCGGAGCGTAGCGCCAATGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCTAATAAAAATAACAACGGCA‑TTTGc                                                                                                                                                      <  7:241893/148‑1 (MQ=255)   aGCGGCGGAGCGTAGCGCCACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGATAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGa                                                                                                                                                   >  8:68845/1‑149 (MQ=255)    gcggcgGAGCGTAGCGCCACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGa                                                                                                                                                   >  6:136421/1‑148 (MQ=255)     cggcggAGCGTAGCGCCACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGat                                                                                                                                                  <  4:98317/148‑1 (MQ=255)      ggcggATCGTAGCGCCACTGAATGATCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGTGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGata                                                                                                                                                 <  8:146151/148‑1 (MQ=255)          gAGCGTAGCGCCACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTCCCGATATTCAg                                                                                                                                            <  7:283061/149‑1 (MQ=255)              gTAGCGCCACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAAt                                                                                                                                        <  4:199406/149‑1 (MQ=255)                 gcgcCACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGc                                                                                                                                                      <  7:272523/132‑1 (MQ=255)                 gcgcCACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGc                                                                                                                                                      >  8:272523/1‑132 (MQ=255)                   gcCACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATgg                                                                                                                                      >  2:176200/1‑146 (MQ=255)                    ccACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATATAAATAACAACGGCA‑TTTGCCGATTTTCAGGAATGGaa                                                                                                                                    >  1:50935/1‑147 (MQ=255)                    ccACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGaa                                                                                                                                    >  1:30748/1‑147 (MQ=255)                    ccACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAaga                                                                                                                                  <  2:30748/149‑1 (MQ=255)                     cACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGAAATACAGGAAAGGAagac                                                                                                                                 >  2:120794/1‑148 (MQ=255)                      aCTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGc                                                                                                                                <  7:103893/149‑1 (MQ=255)                               tctGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCATGTTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCttt                                                                                                                        >  1:283755/1‑149 (MQ=255)                                 tGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGGCGATATTCAg                                                                                                                                            >  5:229140/1‑126 (MQ=255)                                 tGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGGCGATATTCAg                                                                                                                                            >  5:223759/1‑126 (MQ=255)                                 tGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGGCGATATTCAg                                                                                                                                            <  6:229140/126‑1 (MQ=255)                                 tGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGGCGATATTCAg                                                                                                                                            <  6:223759/126‑1 (MQ=255)                                 tGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGc                                                                                                                      >  3:309960/1‑148 (MQ=255)                                   aGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATAAAATAAAAATAACAACGGCA‑TTTGGCGATATTCAGGAATGGAAGAGCATGCGTATGCaaa                                                                                                                   >  8:119999/1‑149 (MQ=255)                                        ctccccccAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCa                                                                                                                <  7:119999/145‑1 (MQ=255)                                            cccaAAAAAACAGAGCTATCATATTTTTTTACAT‑TTAGACTGGTCTCTGTTCTCATATTTTGGCATCTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAgg                                                                                                          <  1:324071/149‑1 (MQ=255)                                              ccATAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAgggg                                                                                                        <  1:120794/149‑1 (MQ=255)                                              ccAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAAGCACTAAgggg                                                                                                        >  1:4091/1‑149 (MQ=255)                                               cAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCATTTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAgggg                                                                                                        <  1:97797/149‑1 (MQ=255)                                                aaaaaCCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGc                                                                                                      <  5:156520/149‑1 (MQ=255)                                                       aGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTg                                                                                               <  2:283755/149‑1 (MQ=255)                                                           cGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTggg                                                                                             <  1:216692/147‑1 (MQ=255)                                                               cacaTTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCg                                                                                       >  6:170314/1‑149 (MQ=255)                                                                   ttttGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGgtcg                                                                                   >  2:258348/1‑149 (MQ=255)                                                                        tttCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGc                                                                              >  5:164599/1‑149 (MQ=255)                                                                        tttCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGc                                                                              >  3:309254/1‑149 (MQ=255)                                                                          taaT‑TTCGACAGGCCTCAGTCAGGAGATTTTGCCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGcc                                                                             <  4:116071/146‑1 (MQ=255)                                                                          tCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGcc                                                                             <  8:180554/148‑1 (MQ=255)                                                                             t‑ttCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGccccc                                                                          >  7:89630/1‑148 (MQ=255)                                                                               ttCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTg                                                                                                                       >  6:106698/1‑102 (MQ=255)                                                                                  gACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAACCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCaaaa                                                                      <  2:4091/148‑1 (MQ=255)                                                                                     aGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGt                                                                  <  3:3777/149‑1 (MQ=255)                                                                                        ccTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTcac                                                               <  1:176200/149‑1 (MQ=255)                                                                                         cTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCa                                                                                                                >  1:270324/1‑99 (MQ=255)                                                                                         cTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCa                                                                                                                <  2:270324/99‑1 (MQ=255)                                                                                         cTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTcaca                                                              >  6:325819/1‑149 (MQ=255)                                                                                              tcagGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGACCCCAAAATCGTCACACCTCa                                                         >  3:16680/1‑149 (MQ=255)                                                                                                agGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCAt                                                        <  4:309254/148‑1 (MQ=255)                                                                                                agGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCAt                                                        <  7:68845/148‑1 (MQ=255)                                                                                                agGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATc                                                       <  6:318166/149‑1 (MQ=255)                                                                                                 ggAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCT‑GGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTa                                                     >  4:18823/1‑149 (MQ=255)                                                                                                  gagaTTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGc                                                                                                                      >  1:166370/1‑84 (MQ=255)                                                                                                  gagaTTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGc                                                                                                                      <  2:166370/84‑1 (MQ=255)                                                                                                  gagaTTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTa                                                     >  6:240137/1‑149 (MQ=255)                                                                                                    gaTTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCAGCCCCCAAAATCGTCACACCTCATCTa                                                     >  5:181399/1‑147 (MQ=255)                                                                                                             tACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGc                                          >  8:372561/1‑149 (MQ=255)                                                                                                             tACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGc                                          >  1:321961/1‑149 (MQ=255)                                                                                                               cTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGa                                        >  8:302847/1‑149 (MQ=255)                                                                                                                  tatccaatatccaaTAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATcg                                     >  1:191175/1‑149 (MQ=255)                                                                                                                   atccaatatccaatAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATcgc                                    >  7:292289/1‑149 (MQ=255)                                                                                                                    tccaatatccaataAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATcgcg                                   <  2:321961/149‑1 (MQ=255)                                                                                                                    tccaatatccaataAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATcgcg                                   >  7:214938/1‑149 (MQ=255)                                                                                                                    tccaatatccaataAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATcgc                                    <  8:292289/148‑1 (MQ=255)                                                                                                                     ccaatatccaataAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATCGCgg                                  >  7:287593/1‑149 (MQ=255)                                                                                                                      caatatccaataAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATCGCgg                                  >  3:352841/1‑148 (MQ=255)                                                                                                                        atatccaataAAAATAACAACGGCT‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATCGCGGtct                               >  6:289073/1‑149 (MQ=255)                                                                                                                              aataAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATCGCGGTCTCTgcgc                         >  8:186816/1‑149 (MQ=255)                                                                                                                                        acaacGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATCGCGGTCTCTGCGCCGGAATGCAt               <  1:18267/149‑1 (MQ=255)                                                                                                                                          aacGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATCGCGGTCTCTGCGCCGGAATGCAt               >  5:181577/1‑147 (MQ=255)                                                                                                                                             ggCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATCGCGGTCTCTGCGCCGGAATGCATCAGCg          <  5:170314/149‑1 (MQ=255)                                                                                                                                                      ccGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATCGCGGTCTCTGCGCCGGAATTCATCAGCGGCCACaaa  >  5:400/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GTAGCGGCGGAGCGTAGCGCCACTGAATGCTCTGAGCTCCGCCCCCCAAAAACCAGTGCGATCACATTTTGTTTCAT‑TTCGACAGGCCTCAGTCAGGAGATTTTGGCTACTGTATCCAATATCCAATAAAAATAACAACGGCA‑TTTGCCGATATTCAGGAATGGAAGAGCATGCCTTTGCAAATCACTAAGGGGGCGACCCTGGGCGAGCGGTCGTCCGCCCCCAAAATCGTCACACCTCATCTACAGATTCGAGCGATCGCGGTCTCTGCGCCGGAATGCATCAGCGGCCACAAA  >  NC_002947/4938676‑4938971

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