Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,395,520	T→C	6.3%	D85G (GAC→GGC)	PP_3871 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,395,520	0	T	C	6.3%	182.7 / 2.5	65	D85G (GAC→GGC)	PP_3871	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base C (2/2);  total (32/33)
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.09e-01

CAGCCGAGAAGCTGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCA  >  NC_002947/4395374‑4395667                                                                                                                                                   |                                                                                                                                                     cagccGAGAAGCTGAAGCGGTAGCCCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTcg                                                                                                                                                    >  7:120087/1‑149 (MQ=255) cagccGAGAAGCTGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGCcg                                                                                                                                                    >  1:322070/1‑149 (MQ=255)          aGCTGAAGCTGTAGACCACCCTCGCAGGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCctgc                                                                                                                                           <  4:210960/149‑1 (MQ=255)             tGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCCTCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCTGCGGCCTGCTgg                                                                                                                                        >  2:269533/1‑149 (MQ=255)              gAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGCCGCCGACctgctgc                                                                                                                                        >  2:44411/1‑147 (MQ=255)              gAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGCCGCCGACctgctgc                                                                                                                                        <  1:44411/148‑2 (MQ=255)               aaGCGTTAGACCAGCCGCTCAATGTTGATCAGTACCACCCCGACAACGTCGTAGACGATGTGGTGGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCCTGCAGCAGGGCGGCGGCCTGCTgg                                                                                                                                        <  6:374879/147‑1 (MQ=255)               aaGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGcc                                                                                                                                      >  4:216435/1‑149 (MQ=255)               aaGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGCCGCCGACctgctgca                                                                                                                                       <  2:322070/148‑3 (MQ=255)               aaGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGAAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTgg                                                                                                                                        <  5:88771/147‑1 (MQ=255)                   ggTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCcgcg                                                                                                                                   >  8:66205/1‑148 (MQ=255)                         ccaccCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTcc                                                                                                                            <  1:113558/149‑1 (MQ=255)                          caccCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCt                                                                                                                           <  3:216435/149‑1 (MQ=255)                            cccGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTg                                                                                                                         <  7:14092/149‑1 (MQ=255)                                  cACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTcggcg                                                                                                                                                 >  7:316540/1‑119 (MQ=255)                                  cACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTcggcg                                                                                                                                                 <  8:316540/119‑1 (MQ=255)                                          aTCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGaa                                                                                                           >  7:356566/1‑149 (MQ=255)                                             aGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAacac                                                                                                        <  5:328431/149‑1 (MQ=255)                                                 ccAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAAcaccac                                                                                                     <  7:238950/148‑1 (MQ=255)                                                 ccAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACaccacc                                                                                                    <  1:387216/149‑1 (MQ=255)                                                     cccGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCtt                                                                                                                          <  3:10461/123‑1 (MQ=255)                                                     cccGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCtt                                                                                                                          >  4:10461/1‑123 (MQ=255)                                                      ccGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCAGAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCTTCCTTGGTTTCCATGGCGAACACCACC‑TCaaa                                                                                              <  4:158214/149‑1 (MQ=255)                                                       cGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAg                                                                                             >  1:82438/1‑149 (MQ=255)                                                        gCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCGTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCGCGGCCTTGGTGTCCATGGCGACCTCCACC‑TCAAAGc                                                                                            >  8:276824/1‑149 (MQ=255)                                                        gCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACAACC‑TCAAAGc                                                                                            >  8:203375/1‑149 (MQ=255)                                                           aCCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGACCACCACC‑TCAAAGCGAt                                                                                         >  5:241441/1‑149 (MQ=255)                                                                 gtAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGa                                                                                   <  3:61787/149‑1 (MQ=255)                                                                   aGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATgg                                                                                                               >  1:13/1‑120 (MQ=255)                                                                   aGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATgg                                                                                                               <  2:13/120‑1 (MQ=255)                                                                      cGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGgtgt                                                                                                                     >  4:319338/1‑111 (MQ=255)                                                                      cGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGgtgt                                                                                                                     <  3:319338/111‑1 (MQ=255)                                                                       gATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCt                                                                             <  7:66205/149‑1 (MQ=255)                                                                          ggggTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTggg                                                                          >  4:128224/1‑149 (MQ=255)                                                                             gTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTaa                                                                       >  2:351708/1‑149 (MQ=255)                                                                                 tACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCAGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGc                                                                   >  5:83259/1‑149 (MQ=255)                                                                                   cTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGt                                                                 >  1:280150/1‑149 (MQ=255)                                                                                    tcatcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGt                                                                 >  8:92911/1‑148 (MQ=255)                                                                                      atcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCt                                                              >  2:187428/1‑149 (MQ=255)                                                                                       tcaCCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTg                                                             >  4:136595/1‑149 (MQ=255)                                                                                        caCCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTgc                                                            >  2:99680/1‑149 (MQ=255)                                                                                         aCCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTgcg                                                           >  1:403372/1‑149 (MQ=255)                                                                                         aCCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTgc                                                            <  2:82438/148‑1 (MQ=255)                                                                                          ccAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTgcgc                                                          <  8:95819/149‑1 (MQ=255)                                                                                           cAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCaaa                                                                                              >  6:293273/1‑112 (MQ=255)                                                                                           cAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCaaa                                                                                              <  5:293273/112‑1 (MQ=255)                                                                                            aGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCCTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCt                                                        >  8:87746/1‑149 (MQ=255)                                                                                                cgccATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGtt                                                    >  2:107528/1‑149 (MQ=255)                                                                                                    aTCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACCATCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTtctc                                                 <  1:99680/149‑1 (MQ=255)                                                                                                     tCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGc                                               <  1:351708/149‑1 (MQ=255)                                                                                                         aCCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGc                                           >  7:249612/1‑149 (MQ=255)                                                                                                            aaGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCttt                                        <  3:128224/149‑1 (MQ=255)                                                                                                             aGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTg                                       >  7:88423/1‑149 (MQ=255)                                                                                                                    ccACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAg                                >  4:31073/1‑149 (MQ=255)                                                                                                                        aaCTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCACCGGCTTTGTCATCAGCGGt                            <  8:356566/149‑1 (MQ=255)                                                                                                                           tCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTtctc                                                 <  2:207602/125‑1 (MQ=255)                                                                                                                           tCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTtctc                                                 >  1:207602/1‑125 (MQ=255)                                                                                                                            ccTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGAt                         >  1:282511/1‑148 (MQ=255)                                                                                                                             cTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGAt                         <  6:83259/147‑1 (MQ=255)                                                                                                                             cTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATca                       <  1:107528/149‑1 (MQ=255)                                                                                                                                 cGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATa                   <  3:136595/149‑1 (MQ=255)                                                                                                                                          tGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAggg          <  8:120087/149‑1 (MQ=255)                                                                                                                                            cagcagGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGt        <  8:88423/149‑1 (MQ=255)                                                                                                                                                 ggTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGttcttc   <  2:280150/149‑1 (MQ=255)                                                                                                                                                  gTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCa  <  2:282511/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CAGCCGAGAAGCTGAAGCGGTAGACCACCCGCGCACGGTTGATCAGTACCAACCCGCCACCGTCGTAGACGATGGGGTCGTACTCATCACCAGGCCGCCATCCGACCAAGGCCCGCCACAACTCCTTGCGAAAGTCGTGCAGCAGGTCGGCGGCCTGCTGGCCGCGCTCGTCCTTGGTGTCCATGGCGAACACCACC‑TCAAAGCGATCCTCGATGTCCTGGGTAACTGCGTTCTGCGCCTTGTTCTCGCCGGCTTTGTCATCAGCGGTGATCACATAGGCCGAGGGGTTCTTCA  >  NC_002947/4395374‑4395667

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