Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	807,900	A→T	5.7%	L45Q (CTG→CAG)	PP_0694 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	807,900	0	A	T	5.7%	270.0 / 4.4	88	L45Q (CTG→CAG)	PP_0694	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (44/39);  new base T (3/2);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGG  >  NC_002947/807753‑808041                                                                                                                                                     |                                                                                                                                               cccTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAg                                                                                                                                               >  6:170569/1‑149 (MQ=255)    tGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCCACAGGCGCTGCGCGGTCTTGATGACCTGCGCGGCGGTCATGTCCCGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGATCAGG‑Gt                                                                                                                                           >  7:165315/1‑149 (MQ=255)         aGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCgg                                                                                                                                       >  6:4389/1‑148 (MQ=255)          gCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTa                                                                                                                                     >  2:310399/1‑149 (MQ=255)            aaGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGTACAGG‑GTGCGGTACt                                                                                                                                   <  7:349652/149‑1 (MQ=255)            aaGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGCTCATGTCCGGCACGATGTAGTAGAAACCACGGGTATCGCACACAAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACt                                                                                                                                   >  6:366928/1‑149 (MQ=255)              ggTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTgg                                                                                                                                 >  3:319049/1‑149 (MQ=255)                 gCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTcc                                                                                                                              >  1:151453/1‑149 (MQ=255)                  cTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                             <  4:61754/149‑1 (MQ=255)                   tCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACAACAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTg                                                                                                                            >  5:183270/1‑149 (MQ=255)                     gACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                             >  8:252142/1‑146 (MQ=255)                     gACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                             <  7:252142/146‑1 (MQ=255)                      aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTgg                                                                                                                           >  3:32476/1‑147 (MQ=255)                      aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATAccc                                                                                                                                  >  4:223967/1‑138 (MQ=25)                      aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATAccc                                                                                                                                  <  3:223967/140‑3 (MQ=25)                         aGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGtt                                                                                                                      <  4:32476/149‑1 (MQ=255)                            aTGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGa                                                                                                                   <  2:116847/149‑1 (MQ=255)                             tGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGa                                                                                                                   >  2:336574/1‑148 (MQ=255)                             tGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACAAGTTCG‑AAGCCCAGGCCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGa                                                                                                                   >  5:194476/1‑148 (MQ=255)                                     aaCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTggg                                                                                                           <  4:319049/148‑1 (MQ=255)                                     aaCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGc                                                                                                          <  5:170569/149‑1 (MQ=255)                                     aaCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGc                                                                                                          >  1:355153/1‑149 (MQ=255)                                      aCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGc                                                                                                          <  1:171383/148‑1 (MQ=255)                                         gCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTAc                                                                                                      <  6:265899/149‑1 (MQ=255)                                          cGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTa                                                                                                       >  3:142164/1‑147 (MQ=255)                                                gcgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGTCAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCa                                                                                                <  1:310399/149‑1 (MQ=255)                                                 cgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAgg                                                                                              <  2:151453/149‑1 (MQ=255)                                                 cgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAgg                                                                                              <  8:76081/149‑1 (MQ=255)                                                    tGCGCGGTCTCGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCGGGGTc                                                                                           >  3:253965/1‑149 (MQ=255)                                                      cgcgGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTcgc                                                                                         <  2:84123/149‑1 (MQ=255)                                                       ggggTGTTGTTGACCTGCGCGGGGGTCATGTCCGGCACGAAGTATTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTcgcg                                                                                        <  6:194476/147‑1 (MQ=255)                                                         ggTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCAGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGcc                                                                                      >  2:143636/1‑149 (MQ=255)                                                             ttgttgACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                                  >  7:10279/1‑149 (MQ=255)                                                             gacatgACCGCCGCGCAGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                                  >  4:93474/5‑149 (MQ=255)                                                             gacatgACCGCCGCGCAGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                                  >  2:221444/5‑149 (MQ=255)                                                                   aCCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTg                                                                            <  1:336574/149‑1 (MQ=255)                                                                   aCCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTAG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGCAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTg                                                                            >  6:289396/1‑149 (MQ=255)                                                                     cTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTgg                                                                           >  3:35609/1‑148 (MQ=255)                                                                     cTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTgg                                                                           >  8:67872/1‑148 (MQ=255)                                                                        cgcgGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAg                                                                        <  5:66872/148‑1 (MQ=255)                                                                          cggcggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCt                                                                     >  8:270216/1‑149 (MQ=255)                                                                            gcggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                                   <  7:67872/149‑1 (MQ=255)                                                                            gcggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                                   <  2:138920/149‑1 (MQ=255)                                                                             cggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCGGGGTCGCGCCGGCTGACGTGGAAGCCTTTg                                                                  <  4:253965/149‑1 (MQ=255)                                                                             cggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GGGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTg                                                                  <  4:51058/149‑1 (MQ=255)                                                                                tCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGaa                                                                >  2:60983/1‑148 (MQ=255)                                                                                   tGTCCGGAACTAAGTAGTAGAAACCCCGGGTATCGCACACCCGTTCG‑AAGCCCCGGGCCTTGAACAGG‑GGGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGc                                                            <  8:28121/149‑1 (MQ=255)                                                                                    gTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGc                                                            >  4:97630/1‑148 (MQ=255)                                                                                     tCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                                   <  5:87628/140‑1 (MQ=255)                                                                                     tCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                                   >  6:87628/1‑140 (MQ=255)                                                                                       cGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAgg                                                                                                                                              >  5:26893/1‑64 (MQ=255)                                                                                       cGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTcgc                                                        <  1:60983/149‑1 (MQ=255)                                                                                        ggCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTTACGTGGAAGCCTTTGAACAGCTcgcg                                                       <  6:183270/149‑1 (MQ=255)                                                                                          cACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                                      >  7:284163/1‑148 (MQ=255)                                                                                          cACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                                      >  2:286983/1‑148 (MQ=255)                                                                                          cACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGa                                                     <  5:4389/149‑1 (MQ=255)                                                                                             gAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                                      >  2:141098/1‑145 (MQ=255)                                                                                                  agtagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCgg                                              <  4:35609/148‑1 (MQ=255)                                                                                                   gtagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcg                                            >  7:220357/1‑149 (MQ=255)                                                                                                   gtagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgc                                             <  7:270216/148‑1 (MQ=255)                                                                                                    tagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcgc                                           >  1:149557/1‑149 (MQ=255)                                                                                                      gAAACCACGGGTATCGCACACCGGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTGGCGGAAGATCGgcgc                                           >  7:59424/1‑147 (MQ=255)                                                                                                      gAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcgc                                           >  6:33613/1‑147 (MQ=255)                                                                                                         aCCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCt                                      >  2:342658/1‑149 (MQ=255)                                                                                                          ccACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCTGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTg                                     <  1:286983/149‑1 (MQ=255)                                                                                                           cACGGGCATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGGCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTgg                                    <  8:158044/149‑1 (MQ=255)                                                                                                              cagTACCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGa                                  >  1:56034/3‑148 (MQ=25)                                                                                                               agTACCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGa                                  >  3:95580/2‑147 (MQ=25)                                                                                                                 tATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGt                              >  5:21947/1‑149 (MQ=255)                                                                                                                   tCGCACACCATAGCG‑TACCCCAGGGCCTTGCACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACACCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGtt                             <  8:165315/148‑1 (MQ=255)                                                                                                                   tCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGTCAGtt                             <  8:284163/148‑1 (MQ=255)                                                                                                                     gCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTc                            >  2:82327/1‑147 (MQ=255)                                                                                                                     gCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTc                            <  1:82327/147‑1 (MQ=255)                                                                                                                     gCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCgg                          <  1:141098/149‑1 (MQ=255)                                                                                                                      cacaCCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGa                         <  3:97630/149‑1 (MQ=255)                                                                                                                         ccctgTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGaaa                       <  4:95580/145‑1 (MQ=37)                                                                                                                         acCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGTAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAg                      <  8:220357/149‑1 (MQ=255)                                                                                                                         acCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGa                         >  3:37701/1‑146 (MQ=255)                                                                                                                          ccAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGa                     <  2:355153/149‑1 (MQ=255)                                                                                                                           cAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTCAACAGCTCGCGGAAGATCGGCGCGAGCTGGTACAGTTCGGAAAGa                     >  4:49877/1‑148 (MQ=255)                                                                                                                              ttCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCaa                 <  5:33613/149‑1 (MQ=255)                                                                                                                                 g‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCAGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGac              >  3:128662/1‑148 (MQ=255)                                                                                                                                   aaGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATg             >  8:50674/1‑149 (MQ=255)                                                                                                                                       ccAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAg                <  7:190377/142‑1 (MQ=255)                                                                                                                                       ccAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAg                >  8:190377/1‑142 (MQ=255)                                                                                                                                              cTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGTGAGCTGGGACAGTTCGGAAAGATCAAGATGCAtgggctggg  >  2:28929/1‑149 (MQ=255)                                                                                                                                              cTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCAtgggctggg  <  4:37701/149‑1 (MQ=255)                                                                                                                                              cTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCAtgggctggg  >  4:101239/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                               CCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGG  >  NC_002947/807753‑808041

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