Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,013,339	G→A	12.5%	A38T (GCG→ACG)	PP_5678 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,013,339	0	G	A	12.5%	150.9 / 14.0	64	A38T (GCG→ACG)	PP_5678	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base A (4/4);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-01

CGAGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTACTCTCA  >  NC_002947/5013190‑5013485                                                                                                                                                      |                                                                                                                                                   cGAGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAgcg                                                                                                                                                    >  1:1164249/1‑150 (MQ=255)     cAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCacggc                                                                                                                                                <  2:440552/150‑1 (MQ=255)      aGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCgcggcg                                                                                                                                               >  1:819985/1‑150 (MQ=255)       gTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGt                                                                                                                                              <  2:672510/150‑1 (MQ=255)        tGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTg                                                                                                                                             <  2:91705/150‑1 (MQ=255)          caacaaGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCACGGCGTgct                                                                                                                                           <  2:1382665/150‑3 (MQ=255)           aacaaGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTg                                                                                                                                             >  2:451115/1‑147 (MQ=255)           aacaaGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTg                                                                                                                                             <  1:451115/147‑1 (MQ=255)                 ccTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAgg                                                                                                                                    <  1:989120/150‑1 (MQ=255)                 ccTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAgg                                                                                                                                    <  1:792638/150‑1 (MQ=255)                  cTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGc                                                                                                                                   >  2:1052034/1‑150 (MQ=255)                      gACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCgc                                                                                                                               >  1:810283/1‑150 (MQ=255)                         ttAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCACGGCGTgctgt                                                                                                                                         >  2:121569/1‑133 (MQ=255)                         ttAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCACGGCGTgctgt                                                                                                                                         <  1:121569/137‑5 (MQ=255)                            ggggAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAgg                                                                                                                                    <  2:523656/139‑1 (MQ=255)                            ggggAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAgg                                                                                                                                    >  1:523656/1‑139 (MQ=255)                            ggggAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCaaa                                                                                                                         >  1:1029926/1‑150 (MQ=255)                                    ccAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGaa                                                                                                                 >  2:539210/1‑150 (MQ=255)                                      aaTGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCAAATACCAGCTCGCTGTTGGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGc                                                                                                               <  1:247797/150‑1 (MQ=255)                                         gatgatCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTa                                                                                                            <  1:1290443/150‑1 (MQ=255)                                          atgatCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTAt                                                                                                           <  1:1090122/150‑1 (MQ=255)                                                aaaGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAGCCCCGAAGCCTATGAACTc                                                                                                     >  2:800291/1‑150 (MQ=255)                                                aaaGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAACGCGACGTGGAAATCAGACTCGCAACAAAACCCCAAAGCCTATAAACTc                                                                                                     >  2:435574/1‑150 (MQ=255)                                                 aaGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCa                                                                                                                                      <  1:279617/116‑1 (MQ=255)                                                 aaGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCa                                                                                                                                      >  2:279617/1‑116 (MQ=255)                                                 aaGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTcc                                                                                                    >  2:1056409/1‑150 (MQ=255)                                                           ccTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCga                                                                                          >  1:427070/1‑150 (MQ=255)                                                           ccTCGCAGCCGTTGCAGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGAAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCgc                                                                                          >  2:124668/1‑149 (MQ=255)                                                             tcacAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCacggc                                                                                                                                                <  1:558388/91‑1 (MQ=255)                                                             tcacAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCacggc                                                                                                                                                >  2:558388/4‑94 (MQ=255)                                                             tCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCAt                                                                                                  <  1:641462/140‑1 (MQ=255)                                                             tCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCAt                                                                                                  >  2:641462/1‑140 (MQ=255)                                                               gCAGCCGTTGCAGTCAGCGAAGCTAGCGACCGCTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAgc                                                                                      >  2:642056/1‑150 (MQ=255)                                                                        gCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCTAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCttt                                                                             <  2:819985/150‑1 (MQ=255)                                                                            tCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCa                                                                         <  1:1056409/150‑1 (MQ=255)                                                                                          gACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGc                                                                                                               <  2:1349330/98‑1 (MQ=255)                                                                                          gACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGc                                                                                                               >  1:1349330/1‑98 (MQ=255)                                                                                          gACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATa                                                           >  1:301335/1‑150 (MQ=255)                                                                                          gACCGTTGCGGCCAGACTTGTCCTTATGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATa                                                           >  2:511762/1‑150 (MQ=255)                                                                                            ccGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATaaa                                                         <  2:1364081/150‑1 (MQ=255)                                                                                            ccGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATaaa                                                         <  2:427070/150‑1 (MQ=255)                                                                                               ttGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTa                                                                                                            >  2:576847/1‑96 (MQ=255)                                                                                               ttGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTa                                                                                                            <  1:576847/96‑1 (MQ=255)                                                                                                      cAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGc                                               >  1:257860/1‑150 (MQ=255)                                                                                                        gACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCTAGCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTg                                             <  1:1169909/150‑1 (MQ=255)                                                                                                        gACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTg                                             <  2:333227/150‑1 (MQ=255)                                                                                                           ttGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTgg                                            >  1:207190/1‑148 (MQ=255)                                                                                                           ttGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTgg                                            <  2:207190/148‑1 (MQ=255)                                                                                                                 ttCTGATGACGGCCATCCACAGCACGCCGTGCTGCGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAAtt                                    <  1:435574/150‑1 (MQ=255)                                                                                                                 ttCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAAtt                                    >  2:1105920/1‑150 (MQ=255)                                                                                                                 ttCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAAtt                                    >  2:832059/1‑150 (MQ=255)                                                                                                                    tgatgaCGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGAt                                 >  2:1023436/1‑150 (MQ=255)                                                                                                                                tgatttccACGCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGa                     >  2:906291/8‑150 (MQ=255)                                                                                                                                 gatttccACGCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGaa                    <  1:906291/144‑1 (MQ=255)                                                                                                                                 ccACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGaa                    <  1:1023436/150‑1 (MQ=255)                                                                                                                                 ccACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGaa                    <  2:301335/150‑1 (MQ=255)                                                                                                                                   tttccACGCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCttt                                                                             <  1:1285772/87‑1 (MQ=38)                                                                                                                                   tttccACGCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCttt                                                                             >  2:1285772/5‑91 (MQ=38)                                                                                                                                      gCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGa               <  1:539210/150‑1 (MQ=255)                                                                                                                                      gCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGa               <  1:1052034/150‑1 (MQ=255)                                                                                                                                      gCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTAAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGa               <  1:1105920/150‑1 (MQ=255)                                                                                                                                      ccACGCCGTGCGGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGACTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGa               <  1:124668/149‑1 (MQ=255)                                                                                                                                        aCGCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGAAAAGCTTTGAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACt             >  2:903214/1‑150 (MQ=255)                                                                                                                                           ccGTGCTGAGCGCGGCGTGGAAATCAGGNTCGCAGCANAACCCCGAANCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAACGTGAACGCTGGGCCCAATTGATGTTCTGTGAGGAAGTTGGCTCCt          >  1:522583/1‑150 (MQ=255)                                                                                                                                            cGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAGCCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCtt         <  1:800291/150‑1 (MQ=255)                                                                                                                                            cGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATTTTCTCTGAGGAAGCTGACTCCtt         <  1:389991/150‑1 (MQ=255)                                                                                                                                            cGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCtt         >  2:389991/1‑150 (MQ=255)                                                                                                                                                 tGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACCCGTTActct    >  1:241775/1‑150 (MQ=255)                                                                                                                                                   aGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTACTCTCa  <  2:1029926/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CGAGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTACTCTCA  >  NC_002947/5013190‑5013485

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