Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	796,845:2	+C	6.6%	intergenic (+44/‑146)	PP_0681 → / → PP_0682	hypothetical protein/inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	796,845	2	.	C	6.6%	240.4 / 11.4	61	intergenic (+44/‑146)	PP_0681/PP_0682	hypothetical protein/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (29/28);  new base C (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

GATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCAT  >  NC_002947/796697‑796992                                                                                                                                                         |                                                                                                                                                    gATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGCGCATCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑t                                                                                                                                                    >  1:1121884/1‑150 (MQ=255)         gCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGTGCCGCTTTGCGGCACATCGCGACACAAGGCC‑‑GC‑‑TCCCACAgg                                                                                                                                            <  1:297046/150‑1 (MQ=255)         gCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAgg                                                                                                                                            >  1:1193601/1‑150 (MQ=255)         gCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAgg                                                                                                                                            >  1:1193631/1‑150 (MQ=255)            gATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTAc                                                                                                                                         >  2:3168/1‑150 (MQ=255)             aTCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACa                                                                                                                                        <  2:1193631/150‑1 (MQ=255)             aTCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACa                                                                                                                                        <  2:1193601/150‑1 (MQ=255)             aTCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGCGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACa                                                                                                                                        <  2:113122/150‑1 (MQ=255)                 ggTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCAt                                                                                                                                    >  1:1268416/1‑150 (MQ=255)                     tCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGc                                                                                                                                >  1:458387/1‑150 (MQ=255)                         cATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGc                                                                                                                            <  1:3168/150‑1 (MQ=255)                          aTTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCt                                                                                                                           <  2:957697/150‑1 (MQ=255)                           tttGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCtg                                                                                                                          >  1:561551/1‑150 (MQ=255)                             tacgttccatgaagttcaggactcatcaggtgactaaaaaagcgacggcgagtatagccagccacctgttacagcaaaagcccgggggctgctccgcagcctgATCGCGACACAAGGCC‑‑ACACTCCCACAGGTACAGCATCGGCTTGCtg                                                                                                                          >  2:749493‑M1/104‑150 (MQ=255)                              gCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGAAACTTGGGGGATTACCAGTTTCATCCTGAGACCTGAAAACAATGTAGGTCTATTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAACATAGGCTTGCTGTgg                                                                                                                       >  1:611400/1‑150 (MQ=25)                                  gcgcAGTTGCGGCGGATCAATGGTGAGCGGTGTAGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAgt                                                                                                                   >  2:1124657/1‑150 (MQ=255)                                    gcAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAgtgt                                                                                                                 >  2:726453/1‑150 (MQ=255)                                    gcAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAgtgt                                                                                                                 >  2:599334/1‑150 (MQ=255)                                     cAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGtgtg                                                                                                                <  1:128568/150‑1 (MQ=255)                                      tgaagttcaggactcatcaggtgactanaaaagcgacggcnagtatanccagccacctgttacagcaaaagcccgggggctgctccgcagcctgATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTgg                                                                                                               <  1:783584‑M1/56‑1 (MQ=255)                                        ttGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGcc                                                                                                             <  1:486764/150‑1 (MQ=255)                                            ggcggATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGc                                                                                                                            >  1:628282/1‑131 (MQ=255)                                            ggcggATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGc                                                                                                                            <  2:628282/131‑1 (MQ=255)                                            ggcggATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTtgt                                                                                                         <  1:1088093/150‑1 (MQ=255)                                                gactcatcaggtgactaaaaaagcgacggcgagtatagccagccacctgttacagcaaaagcccgggggctgctccgcagcctgATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTcg                                                                                                     >  2:783584‑M1/85‑150 (MQ=255)                                                  tCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTcgcg                                                                                                   >  1:677063/1‑150 (MQ=255)                                                     atcaggtgactaaaaaagcgacggcgagtatagccagccacctgttacagcaaaagcccgggggctgctccgcagcctgATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCtt                                                                                                           >  2:794089‑M1/80‑139 (MQ=255)                                                     atcaggtgactaaaaaagcgacggcgagtatagccagccacctgttacagcaaaagcccgggggctgctccgcagcctgATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCtt                                                                                                           <  1:794089‑M1/60‑1 (MQ=255)                                                     atcaggtgactaaaaaagcgacggcgagtatagccagccacctgttacagcaaaagcccgggggctgctccgcagcctgATCGCGACACAAGGCC‑‑ACACTCCCACAGGTACAGCATCGGCTTACTGTGGGa                                                                                                                     <  1:390378‑M1/52‑1 (MQ=255)                                                     atcaggtgactaaaaaagcgacggcgagtatagccagccacctgttacagcaaaagcccgggggctgctccgcagcctgATCGCGACACAAGGCC‑‑ACACTCCCACAGGTACAGCATCGGCTTACTGTGGGa                                                                                                                     >  2:390378‑M1/80‑131 (MQ=255)                                                     aTGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATc                                                                                                <  2:677063/150‑1 (MQ=255)                                                      tGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCa                                                                                               <  2:1268416/150‑1 (MQ=255)                                                              ggTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCgg                                                                                       >  1:718337/1‑150 (MQ=255)                                                                      aaCCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGAt                                                                                                 <  2:1022232/132‑1 (MQ=255)                                                                      aaCCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGAt                                                                                                 >  1:1022232/1‑132 (MQ=255)                                                                               gATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTTCGGAGCAGCCCCCGGGCttt                                                                      >  1:1025869/1‑150 (MQ=255)                                                                                aTTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCtttt                                                                     >  2:874688/1‑150 (MQ=255)                                                                                   aCCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCt                                                                  <  2:458387/150‑1 (MQ=255)                                                                                    ccAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTg                                                                 >  1:572143/1‑150 (MQ=255)                                                                                    ccAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTg                                                                 >  1:607413/1‑150 (MQ=255)                                                                                     cAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGt                                                                <  2:78116/150‑1 (MQ=255)                                                                                          tCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAg                                                           <  2:561551/150‑1 (MQ=255)                                                                                                 ttacagcaaaagcccgggggctgctccgcagcctgATCGCGACACAAGGCCACAC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAgt                                                                                                                   >  1:84683‑M1/36‑89 (MQ=255)                                                                                                 ttacagcaaaagcccgggggctgctccgcagcctgATCGCGACACAAGGCC‑‑ACACTCCCACAGGTACAGCATCGGCTTGCTGTGGGAgt                                                                                                                   <  2:84683‑M1/54‑1 (MQ=255)                                                                                                  agaCCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAACATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGtggctgg                                                   <  2:611400/150‑1 (MQ=255)                                                                                                             aTGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGAC‑‑CC‑‑TCCCACAGGTATAGCATCGGCTTGCTGTGGGAGTGTTGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGcc                                        <  1:726453/150‑1 (MQ=255)                                                                                                                   gggCCGCTGTGCGGCCCATCGCGACACAAGGCC‑‑GT‑‑TCCCACAGGAACAGCATTTGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCt                                  <  1:1048671/150‑1 (MQ=39)                                                                                                                     gCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCttt                                >  1:812945/1‑150 (MQ=255)                                                                                                                          tttGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAg                           >  1:870534/1‑150 (MQ=255)                                                                                                                           ttGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTAGCAGTCTATACTCGCCGTCGCTTTTTTAGt                          <  1:874688/150‑1 (MQ=255)                                                                                                                              cGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCGGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCAc                       <  1:1124657/150‑1 (MQ=255)                                                                                                                                gCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCt                     <  2:572143/150‑1 (MQ=255)                                                                                                                                  ccATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATAc                                             <  1:266383/124‑1 (MQ=255)                                                                                                                                  ccATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATAc                                             >  2:266383/1‑124 (MQ=255)                                                                                                                                    ATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCtgatg                 <  2:1256580‑M1/150‑1 (MQ=255)                                                                                                                                           CACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTc                                     >  2:657064‑M1/1‑123 (MQ=255)                                                                                                                                           CACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTc                                     <  1:657064‑M1/123‑1 (MQ=255)                                                                                                                                              AAGGCC‑‑GT‑‑TCCCACTGTTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCTCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAAc       >  1:177995‑M1/1‑150 (MQ=255)                                                                                                                                                 GCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCta                                                <  2:772054‑M1/106‑1 (MQ=255)                                                                                                                                                 GCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCta                                                >  1:772054‑M1/1‑106 (MQ=255)                                                                                                                                                  CC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCa   >  2:870534‑M1/1‑150 (MQ=255)                                                                                                                                                   C‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCAt  >  2:718337‑M1/1‑150 (MQ=255)                                                                                                                                                         |                                                                                                                                                    GATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCC‑‑GC‑‑TCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCAT  >  NC_002947/796697‑796992

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