Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,955,965	T→A	8.2%	A201A (GCT→GCA)	PP_1753 →	DMT superfamily permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,955,965	0	T	A	8.2%	178.7 / 3.5	60	A201A (GCT→GCA)	PP_1753	DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base T (30/25);  new base A (2/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATTGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGAGGCGTCGCG  >  NC_002947/1955824‑1956108                                                                                                                                              |                                                                                                                                                  ttCATTGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAgcagcgcagcgcagc                                                                                                                                            >  7:287712/1‑149 (MQ=255)     ttGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCAGCGCAGCGCAGCTGag                                                                                                                                        <  8:287712/149‑2 (MQ=255)     ttGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCAGCGCAGCGCAGCTGag                                                                                                                                        <  1:302101/149‑2 (MQ=255)     ttGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGAGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGac                                                                                                                                        <  4:248795/149‑1 (MQ=255)         gTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCt                                                                                                                                   >  8:150582/1‑149 (MQ=255)             cTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCgctgct                                                                                                                                           <  5:308379/138‑1 (MQ=255)             cTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCgctgct                                                                                                                                           >  6:308379/1‑138 (MQ=255)               acacCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTg                                                                                                                             >  5:188781/1‑149 (MQ=255)                caccacCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGt                                                                                                                              <  7:163440/147‑1 (MQ=255)                          gtgtTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGa                                                                                                                 <  1:341847/149‑1 (MQ=255)                               gACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGcc                                                                                                             <  5:265454/148‑1 (MQ=255)                                aCCTTCTTGCTGGCGTCCAGCATCTGGGTGTGCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTTACT‑CCTTGGGTGCA‑CCCGCTGGAGGCCtt                                                                                                           <  1:364629/149‑1 (MQ=255)                                 ccATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCTTTg                                                                                                          <  4:7870/149‑1 (MQ=255)                                 ccATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTg                                                                                                          <  3:148603/149‑1 (MQ=255)                                  cATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGCCCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTg                                                                                                          >  2:266121/1‑148 (MQ=255)                                  cATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGc                                                                                                         <  7:310650/149‑1 (MQ=255)                                       tGCTGGCGACCAGCATCTGGGTGTTCTTCGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACt                                                                                                     >  1:16337/1‑148 (MQ=255)                                        gCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCaa                                                                                                       <  3:125237/145‑1 (MQ=255)                                          tGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGAcc                                                                                                 >  1:164866/1‑149 (MQ=255)                                             cGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCg                                                                                              >  1:277049/1‑149 (MQ=255)                                              gACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCgg                                                                                             >  5:143335/1‑149 (MQ=255)                                               aCCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGTGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGt                                                                                            <  2:164866/149‑1 (MQ=255)                                                ccAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACt                                                                                                     >  7:178272/1‑139 (MQ=255)                                                ccAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACt                                                                                                     <  8:178272/139‑1 (MQ=255)                                                   gCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCaa                                                                                        <  8:101379/149‑1 (MQ=255)                                                     aTCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAgg                                                                                      >  5:203624/1‑149 (MQ=255)                                                             gtTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGgctggctg                                                                               <  8:206410/148‑1 (MQ=255)                                                             gtTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGgctggctg                                                                               >  7:206410/1‑148 (MQ=255)                                                              ttCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGgctgct                                                                             >  4:315875/1‑149 (MQ=255)                                                                cTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGgctgctgc                                                                           >  2:291796/1‑149 (MQ=255)                                                                cTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCCGTGCAAGGCTGgctgctgc                                                                           >  5:101481/1‑149 (MQ=255)                                                                 tACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGgctgctgct                                                                          >  3:85824/1‑149 (MQ=255)                                                                   cGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGctgctgctg                                                                         <  7:150582/148‑1 (MQ=255)                                                                    gCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCACCCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAg                                                                                       <  1:239690/134‑1 (MQ=255)                                                                    gCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAg                                                                                       >  2:239690/1‑133 (MQ=255)                                                                    gCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGc                                                                       <  2:277049/149‑1 (MQ=255)                                                                        tGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCtgt                                                                   >  3:185483/1‑149 (MQ=255)                                                                              aGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCGTGTTGTCt                                                              >  2:289316/1‑148 (MQ=255)                                                                               gAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGa                                                            <  6:143335/149‑1 (MQ=255)                                                                                     gTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACAccc                                                      >  8:319910/1‑149 (MQ=255)                                                                                        gCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACACCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCtgg                                                   >  8:201708/1‑149 (MQ=255)                                                                                            aCCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGgctggctg                                                                               <  1:46084/117‑1 (MQ=255)                                                                                            aCCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGgctggctg                                                                               >  2:46084/1‑117 (MQ=255)                                                                                            aCCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCtggtgg                                                <  3:315875/148‑1 (MQ=255)                                                                                             ccGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGc                                                                       >  4:15327/1‑124 (MQ=255)                                                                                             ccGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGc                                                                       <  3:15327/124‑1 (MQ=255)                                                                                                    gCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGgcg                                       >  4:169320/1‑149 (MQ=255)                                                                                                     cACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGgcgc                                      >  3:227880/1‑149 (MQ=255)                                                                                                               agGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGCCCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGaa                            >  4:80004/1‑149 (MQ=255)                                                                                                                ggTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGCCGCGTTCGCCGAAg                           <  7:319910/149‑1 (MQ=255)                                                                                                                    atgatgGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCt                       <  4:85824/149‑1 (MQ=255)                                                                                                                      gatgGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTgg                     >  6:344052/1‑149 (MQ=255)                                                                                                                                 gtgtCAGCAGCGCAGCGCAGCTGACAGCCTTGGGTGCA‑CCCGCTGGAGGCCTTGCACCTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTCACGGCGCGTTCGCCGAAGCGCTGGCGCACTgg             >  1:195616/4‑147 (MQ=255)                                                                                                                                 aCCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCt                                              <  6:55013/113‑1 (MQ=255)                                                                                                                                 aCCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCt                                              >  5:55013/1‑113 (MQ=255)                                                                                                                                  tgtCAGCAGCGCAGCGCAGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGAgg         >  6:50293/3‑149 (MQ=255)                                                                                                                                  ccacAGCTGCCCAGCGCTGCTGACA‑CCTTCGGTCCA‑CCCGCTGCAGGCCATGCAACTCACCCCGGAGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTg                                                             <  4:222781/94‑1 (MQ=255)                                                                                                                                  ccTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTg                                                             >  3:222781/1‑97 (MQ=255)                                                                                                                                  ccTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTgg             <  5:344052/145‑1 (MQ=255)                                                                                                                                  ccTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGAgg         >  2:187346/1‑149 (MQ=255)                                                                                                                                     cAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGAGGCGt      >  5:62614/1‑149 (MQ=255)                                                                                                                                         tgcgctgcgctgcTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGAGGCGTcgcg  <  3:282430/149‑1 (MQ=255)                                                                                                                                           cgctgcgctgcTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTACTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGAGGCGTcgcg  <  4:153482/147‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                  TTCATTGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCACTCGCAGCAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACA‑CCTTGGGTGCA‑CCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGAGGCGTCGCG  >  NC_002947/1955824‑1956108

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