Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,622,921	C→G	7.3%	A146P (GCC→CCC)	PP_1422 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,622,921	0	C	G	7.3%	209.3 / 5.7	68	A146P (GCC→CCC)	PP_1422	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base G (3/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01

GCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAGGCCGATGGC  >  NC_002947/1622773‑1623068                                                                                                                                                      |                                                                                                                                                     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                                                                                                                                                     >  5:69910/1‑137 (MQ=255)              tGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGcc                                                                                                                                                     <  6:69910/137‑1 (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                                                                                                                                         <  1:297642/124‑1 (MQ=255)                                       cTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTg                                                                                                                                         >  2:297642/1‑124 (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                                                                                             >  6:20878/1‑149 (MQ=255)                                                         gCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCa                                                                                             >  3:83263/1‑149 (MQ=255)                                                                gggTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAg                                                                                      <  5:20878/149‑1 (MQ=255)                                                                gggTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAg                                                                                      <  7:141675/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                                                                            >  4:3505/1‑146 (MQ=255)                                                                             gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCTGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                            <  3:3505/146‑1 (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                                                                                                          >  1:264856/1‑97 (MQ=255)                                                                                                aTAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAg                                                                                                          <  2:264856/97‑1 (MQ=255)                                                                                                 tAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTc                                                     >  6:268139/1‑149 (MQ=255)                                                                                                 tAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTc                                                     <  6:123397/149‑1 (MQ=255)                                                                                                    aTGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATcgc                                                                  <  8:153442/133‑1 (MQ=255)                                                                                                    aTGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATcgc                                                                  >  7:153442/1‑133 (MQ=255)                                                                                                     tGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                            <  6:201037/122‑1 (MQ=255)                                                                                                     tGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                            >  5:201037/1‑122 (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                                           <  6:235953/127‑1 (MQ=255)                                                                                                                              acccgACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                                           >  5:235953/5‑131 (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)                                                                                                                                           acccaAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATc                <  4:275240/145‑1 (MQ=255)                                                                                                                                           acccaAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATc                >  3:275240/1‑145 (MQ=255)                                                                                                                                             ccaAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAgg          >  1:303403/1‑149 (MQ=255)                                                                                                                                                     gCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAGGCCGATGGc  >  1:102971/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCA‑GCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGT‑GAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAGGCCGATGGC  >  NC_002947/1622773‑1623068

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