Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,622,907	C→G	7.2%	G150G (GGG→GGC)	PP_1422 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,622,907	0	C	G	7.2%	220.9 / 5.7	69	G150G (GGG→GGC)	PP_1422	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base G (3/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

CGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATC  >  NC_002947/1622760‑1623054                                                                                                                                                     |                                                                                                                                                     cGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCTGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGc                                                                                                                                                      <  4:139531/148‑1 (MQ=255)        tCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGcccaccc                                                                                                                                                >  2:60596/1‑147 (MQ=255)          ccAGGCCTTGATTATCTTGACCCACCGAAACGGCTGGAGTTGCTCGATCACATAATGAAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTATGCGGCATAGAGGAGTACCTGCA‑GCGAGCCAAACTGGTTGACCAGcccacccccgc                                                                                                                                            >  8:262845/1‑145 (MQ=255)           cAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGt                                                                                                                                           >  4:18685/1‑149 (MQ=255)            aGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTg                                                                                                                                          <  3:167744/149‑1 (MQ=255)              gCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAACGTGGc                                                                                                                                        <  7:286539/149‑1 (MQ=255)                cTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCt                                                                                                                                      <  8:299979/149‑1 (MQ=255)                   gAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGt                                                                                                                                   >  2:33220/1‑149 (MQ=255)                           tGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGcc                                                                                                                                       <  6:69910/137‑1 (MQ=255)                           tGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGcc                                                                                                                                       >  5:69910/1‑137 (MQ=255)                            gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGGCCACCTTGGGTGGGCTGGTcaacc                                                                                                                              <  5:169254/145‑5 (MQ=255)                             tCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTga                                                                                                                          <  3:18685/148‑1 (MQ=255)                             tCCCCCAGCACCGGCGGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCTGCATAGATGAGTACCGGCA‑GCGCGCGAATCTGTTTGACCAGCCCACGCAAGGTGGCCTGGTCGGGTTTga                                                                                                                          <  8:229294/148‑1 (MQ=255)                               ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATATATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                       >  8:165073/1‑149 (MQ=255)                               ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                       <  1:311532/149‑1 (MQ=255)                                        cGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGc                                                                                                             >  6:293328/1‑149 (MQ=255)                                         ggCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCACTCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGc                                                                                                             <  1:253273/149‑1 (MQ=255)                                          gCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GTAGCGa                                                                                                           >  4:84708/1‑149 (MQ=255)                                                    cTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTg                                                                                                                           >  2:297642/1‑124 (MQ=255)                                                    cTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTg                                                                                                                           <  1:297642/124‑1 (MQ=255)                                                    cTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGAttt                                                                                                 <  5:301701/149‑1 (MQ=255)                                                       gATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCg                                                                                               >  6:291121/1‑148 (MQ=255)                                                         tCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAg                                                                                            <  3:210341/149‑1 (MQ=255)                                                          cacaTAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGc                                                                                           >  3:244302/1‑149 (MQ=255)                                                                aTGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGc                                                                                     >  5:121884/1‑149 (MQ=255)                                                                    aGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGTCAGATTTCGTAGCTGATGCCCGc                                                                                 <  6:121884/149‑1 (MQ=255)                                                                    aGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGc                                                                                 >  8:141675/1‑149 (MQ=255)                                                                    aGGCAAACCGGGTCAGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGc                                                                                 <  7:63485/149‑1 (MQ=255)                                                                      gCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCa                                                                               >  3:83263/1‑149 (MQ=255)                                                                      gCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCa                                                                               >  6:20878/1‑149 (MQ=255)                                                                             gggTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAg                                                                        <  7:141675/149‑1 (MQ=255)                                                                             gggTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAg                                                                        <  5:20878/149‑1 (MQ=255)                                                                               gTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTCGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAgg                                                                       <  1:60596/149‑1 (MQ=255)                                                                                    acaacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGATGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTg                                                                 <  4:244302/149‑1 (MQ=255)                                                                                      aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGc                                                               >  5:123397/1‑149 (MQ=255)                                                                                        caGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                              >  2:59307/1‑148 (MQ=255)                                                                                         aGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCa                                                            <  4:83263/149‑1 (MQ=255)                                                                                          gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCTGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                              <  3:3505/146‑1 (MQ=255)                                                                                          gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCTGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                              >  4:3505/1‑146 (MQ=255)                                                                                          gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGAt                                                                                       >  1:118494/1‑121 (MQ=255)                                                                                          gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGAt                                                                                       <  2:118494/121‑1 (MQ=255)                                                                                             cGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGata                                                        >  1:8547/1‑149 (MQ=255)                                                                                                   gCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGc                                                                                 <  4:302477/118‑1 (MQ=255)                                                                                                   gCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGc                                                                                 >  3:302477/1‑118 (MQ=255)                                                                                                     tctTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAg                                                >  5:28248/1‑149 (MQ=255)                                                                                                       ttCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACtt                                                                            >  5:171623/1‑119 (MQ=255)                                                                                                       ttCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACtt                                                                            <  6:171623/119‑1 (MQ=255)                                                                                                          ggCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTc                                            <  1:59307/148‑1 (MQ=255)                                                                                                          ggCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTc                                            <  5:293328/148‑1 (MQ=255)                                                                                                            cATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGc                                         <  2:8547/149‑1 (MQ=255)                                                                                                             aTAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAg                                                                                            <  2:264856/97‑1 (MQ=255)                                                                                                             aTAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAg                                                                                            >  1:264856/1‑97 (MQ=255)                                                                                                              tAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTc                                       <  6:123397/149‑1 (MQ=255)                                                                                                              tAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTc                                       >  6:268139/1‑149 (MQ=255)                                                                                                                 aTGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATcgc                                                    >  7:153442/1‑133 (MQ=255)                                                                                                                 aTGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATcgc                                                    <  8:153442/133‑1 (MQ=255)                                                                                                                  tGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                              >  5:201037/1‑122 (MQ=255)                                                                                                                  tGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                              <  6:201037/122‑1 (MQ=255)                                                                                                                  ggAGGACCGGCA‑GCGCGCCACTCTGGGTGACCAGCCCACCCAAGGGGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATCTCGCGCAGTTTCTGCTCCACt                                   <  7:262845/148‑1 (MQ=255)                                                                                                                   gAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTg                                  >  8:139300/1‑149 (MQ=255)                                                                                                                        ccGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                             >  2:238266/1‑149 (MQ=255)                                                                                                                         cGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                             >  2:268945/1‑148 (MQ=255)                                                                                                                                   gcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCAt                   <  6:28248/149‑1 (MQ=255)                                                                                                                                        tCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCa                                                                         >  4:198089/1‑90 (MQ=255)                                                                                                                                        tCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCa                                                                         <  3:198089/90‑1 (MQ=255)                                                                                                                                          aacccgACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAAt            >  1:153450/6‑149 (MQ=255)                                                                                                                                           acccgACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                             >  5:235953/5‑131 (MQ=255)                                                                                                                                           acccgACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                             <  6:235953/127‑1 (MQ=255)                                                                                                                                               gACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCt       >  1:144135/1‑149 (MQ=255)                                                                                                                                                aCCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCTGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTc      <  7:165073/149‑1 (MQ=255)                                                                                                                                                    gCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATc  <  1:33220/149‑1 (MQ=255)                                                                                                                                                    gCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACGTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGAt   >  2:239305/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CGGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATC  >  NC_002947/1622760‑1623054

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