Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,543,961	Δ1 bp	7.4%	coding (1517/2607 nt)	PP_3130 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,543,961	0	C	.	7.4%	247.9 / 15.2	81	coding (1517/2607 nt)	PP_3130	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/36);  new base . (3/3);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

CGTCGCTGGTCAGCGGGTTCTTCTCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAAGGCCTTGGCCGCCGTCATGTAGC  >  NC_002947/3543812‑3544108                                                                                                                                                      |                                                                                                                                                        cGTCGCTGGTCAGCGGGTTCTTCTCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGc                                                                                                                                                         >  1:796948/1‑150 (MQ=255)       tGGTCAGCGGGTTCTTCTCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATTGGGACTTTCACGCCAACGCTgcgcc                                                                                                                                                   >  1:869619/1‑149 (MQ=255)         gTCAGCGGGTTCTTCTCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCTGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCacg                                                                                                                                                 >  1:845471/1‑148 (MQ=255)             gCGGGTTCTTCTCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCCGTGAGCTAAGACCCTTCCCTCGCCAAGGGCGTTCAGTACCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTggg                                                                                                                                           <  1:36216/150‑1 (MQ=255)              cGGGTTCTTCTCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑Agtgggc                                                                                                                                           >  2:1826364‑M2/1‑144 (MQ=255)              cGGGTTCTTCTCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑Agtggg                                                                                                                                            >  1:1174459‑M2/1‑144 (MQ=255)              cGGGTTCTTCTCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑Agtggg                                                                                                                                            <  2:1174459‑M2/149‑6 (MQ=255)               gggTTCTTCTCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑Agtgggcg                                                                                                                                          <  2:382384‑M2/150‑8 (MQ=255)                ggTTCTTCTCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGa                                                                                                                                        <  2:796948/150‑1 (MQ=255)                  ttcttcTCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATc                                                                                                                                      >  1:490908/1‑150 (MQ=255)                    cttctCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCAt                                                                                                                                    >  2:98488/1‑150 (MQ=255)                     ttctCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTAGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCCCGCCCACGCTGCGCA‑G‑GTGGGCGCTCAtt                                                                                                                                   >  1:1904024/1‑149 (MQ=255)                      tctcCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACAGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGt                                                                                                                                  >  2:59413/1‑150 (MQ=255)                        tcCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCg                                                                                                                                >  2:1904336/1‑150 (MQ=255)                         ccAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCgg                                                                                                                               >  2:1868/1‑150 (MQ=255)                         caaggcctttgccaaacagtttgacctacccgccgacctgacctggcgtgcctgggagcatgggccgtcgaagctgttgctcaacgagctggaacaacgcttcaatgccgacatgatCGCCCACGCTGCGCA‑A‑GTggg                                                                                                                                           >  2:501299‑M1/118‑138 (MQ=255)                         caaggcctttgccaaacagtttgacctacccgccgacctgacctggcgtgcctgggagcatgggccgtcgaagctgttgctcaacgagctggaacaacgcttcaatgccgacatgatCGCCCACGCTGCGCA‑A‑GTggg                                                                                                                                           <  1:501299‑M1/21‑1 (MQ=255)                                 ccTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAg                                                                                                                       <  2:1689299/150‑1 (MQ=255)                                       tCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTgc                                                                                                                                                      <  2:924797/115‑1 (MQ=255)                                       tCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCgttgt                                                                                                                 <  1:1184346/150‑1 (MQ=255)                                       tCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGACAGCAATGGGGACTTTCACGCCCACGCTgc                                                                                                                                                      >  1:924797/1‑115 (MQ=255)                                          ccGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTcc                                                                                                              >  2:1422307/1‑150 (MQ=255)                                          ccGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTcc                                                                                                              >  1:671675/1‑150 (MQ=255)                                                 tgtgGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGTGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACT‑TGCGCA‑Gcgtgggcgtgaaagtcc                                                                                                                                >  1:612629‑M2/1‑108 (MQ=255)                                                 tgtgGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACT‑TGCGCAGCgtgggcgtgaaagtcc                                                                                                                                 <  2:612629‑M2/125‑17 (MQ=255)                                                  gtgGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGtt                                                                                                    <  1:382384/150‑1 (MQ=255)                                                    ggggGTCGATGACATCGCCAACGCTCAATGCCGGAAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCAGTGAAGCGTTGTTCCAGC‑‑TCGttta                                                                                                  <  2:845471/150‑3 (MQ=255)                                                      gCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGc                                                                                                <  2:35490/150‑1 (MQ=255)                                                        gTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAg                                                                                                            <  2:671970/138‑1 (MQ=255)                                                        gTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAg                                                                                                            >  1:671970/1‑138 (MQ=255)                                                           gATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAg                                                                                           >  2:249571/1‑150 (MQ=255)                                                               aCATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCAGTACCCCCCGCGGGGCCGGTTGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTc                                                                                       <  2:1506414/150‑1 (MQ=255)                                                                cATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCg                                                                                      <  1:1868/150‑1 (MQ=255)                                                                 aTCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGa                                                                                     >  1:284624/1‑150 (MQ=255)                                                                  tCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGAc                                                                                    >  1:1489446/1‑150 (MQ=255)                                                                     ccAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACT‑TGCGCA‑Gcgtgggcgtgaaagtccccatt                                                                                                                           >  2:1605973‑M2/1‑88 (MQ=255)                                                                     ccAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACT‑TGCGCA‑Gcgtgggcgtgaaagtccccatt                                                                                                                           <  1:1605973‑M2/110‑23 (MQ=255)                                                                     ccAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGc                                                                                 >  1:1385935/1‑150 (MQ=255)                                                                      cAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACgg                                                                                  <  2:616215/148‑1 (MQ=255)                                                                      cAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACgg                                                                                  >  1:616215/1‑148 (MQ=255)                                                                      cAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGcc                                                                                <  2:490908/150‑1 (MQ=255)                                                                           cTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGc                                                                           >  2:201021/1‑150 (MQ=255)                                                                            tCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATc                                                                                                                                      <  1:321301/92‑1 (MQ=255)                                                                            tCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATc                                                                                                                                      >  2:321301/1‑92 (MQ=255)                                                                                     gCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCa                                                                                                                             <  1:72573/92‑1 (MQ=255)                                                                                     gCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCa                                                                                                                             >  2:72573/1‑92 (MQ=255)                                                                                     gCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCAc                                                                 >  1:1805158/1‑150 (MQ=255)                                                                                          gcgcTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCcag                                                            <  1:1904336/150‑1 (MQ=255)                                                                                            gcTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCcaggt                                                          >  1:1341745/1‑150 (MQ=255)                                                                                             cTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCcaggtc                                                         >  2:870983/1‑150 (MQ=255)                                                                                              tAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCcaggtca                                                        <  2:442886/150‑1 (MQ=255)                                                                                               aggCCCAGTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCGCGGCCACGCTGCTCA‑A‑GTGGGCGGTCATGTCGTGAGTGAAGCGTGGTTGCAGC‑‑TCGTGGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGGcaggtcag                                                       <  2:119664/148‑1 (MQ=255)                                                                                                  cccATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCaggtcaggtc                                                    >  2:1786365/1‑150 (MQ=255)                                                                                                     aTTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTcggc                                                 >  1:1346669/1‑150 (MQ=255)                                                                                                     aTTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTcggc                                                 <  2:671675/150‑1 (MQ=255)                                                                                                       tACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTcggcgg                                               <  2:1341745/150‑1 (MQ=255)                                                                                                       tACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTcggcgg                                               <  2:1805158/150‑1 (MQ=255)                                                                                                         tcaacgagctggaacaacgcttcaatgccgacatgatCGCCCACGCTGCGCA‑A‑GTGGGCGAt                                                                                                                                       <  1:1086790‑M1/25‑1 (MQ=255)                                                                                                         tcaacgagctggaacaacgcttcaatgccgacatgatCGCCCACGCTGCGCA‑A‑GTGGGCGAt                                                                                                                                       >  2:1086790‑M1/38‑62 (MQ=255)                                                                                                         ccGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTcc                                                                        >  1:641394/1‑123 (MQ=255)                                                                                                         ccGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTcc                                                                        <  2:641394/123‑1 (MQ=255)                                                                                                          cGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAAGAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTa                                            >  1:1876909/1‑150 (MQ=255)                                                                                                          cGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTa                                            <  2:284624/150‑1 (MQ=255)                                                                                                                      tggtggTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGtt                                <  1:59413/150‑1 (MQ=255)                                                                                                                        gtggtAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTg                              <  1:97340/150‑1 (MQ=255)                                                                                                                          ggtAGCAATCGGGACTTTTACGTCCACGCTTCGCC‑A‑GTCGGCGAGCATGTCCGCAATAAAGCGTTCTTCCAGG‑‑CCCTTGAGCAACAGCTTCGACCGCCCATCCTCCCAGACACGCCAGATCAGATCGGCAGGCAGCTCAAAATGTTTCGc                            <  2:869619/150‑1 (MQ=255)                                                                                                                            cttcaatgccgacatgatCGCCCACT‑TGCGCA‑GCGTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGtt                                >  2:1821312‑M1/19‑144 (MQ=255)                                                                                                                            cttcaatgccgacatgatCGCCCACT‑TGCGCA‑GCGTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGtt                                <  1:1821312‑M1/126‑1 (MQ=255)                                                                                                                             aGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCaaa                         <  2:1904024/150‑1 (MQ=255)                                                                                                                                  tGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAAGGCCt                    >  1:692662/1‑150 (MQ=255)                                                                                                                                    gggACTTTCACGCCCACGCTGCGCC‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAAGGCCTTg                  >  1:1022116/1‑150 (MQ=255)                                                                                                                                    gggACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAAgg                       <  2:1216703/145‑1 (MQ=255)                                                                                                                                    gggACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAAgg                       >  1:1216703/1‑145 (MQ=255)                                                                                                                                       aCTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAAGGCCTTGgcc               <  2:1264408/150‑1 (MQ=255)                                                                                                                                         tttCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAAGGCCTTGgccgc             >  1:91733/1‑150 (MQ=255)                                                                                                                                         tttCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCACC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAAGGCCTTGgccgc             >  2:150496/1‑150 (MQ=255)                                                                                                                                         tttCACGCCCACACTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTCAGCAACAGCTTCAACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGCAGGTGAAACTGTTTGGCAAAGGCATTGgccgc             >  1:950622/1‑150 (MQ=255)                                                                                                                                            cACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAAGGCCTTGGCCGCCGt          <  2:1346669/150‑1 (MQ=255)                                                                                                                                             aCGCCCACGCTGCGCA‑C‑GTGGGCGATCATGTCGGCATTGAAGCTTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAAGGCCTTGGCCGCCGTc         >  2:21734/1‑150 (MQ=255)                                                                                                                                              cGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTG‑‑TTGTTTTTCAAGCAGACGGTGGGCAACAGCTTCGCAGATCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAAGGCCTTGGCCGCCGTCa        <  1:435543/150‑1 (MQ=255)                                                                                                                                                 CCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAATGCCTTGGCCGCCGTCATGt     <  2:77757‑M1/150‑1 (MQ=255)                                                                                                                                                    CGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAAGGCCTTGGCCGCCGTCATGTAGc  >  2:1489446‑M1/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                        CGTCGCTGGTCAGCGGGTTCTTCTCCAGTTGGCCTGTGTCGCCGTAGCTGTGGGCGTCGATGACATCGCCAACGCTCAATGCCGGCAGTGCGCTAAGCCCATTACCGCCCCAGGTGCTGGTGGTAGCAATGGGGACTTTCACGCCCACGCTGCGCA‑A‑GTGGGCGATCATGTCGGCATTGAAGCGTTGTTCCAGC‑‑TCGTTGAGCAACAGCTTCGACGGCCCATGCTCCCAGGCACGCCAGGTCAGGTCGGCGGGTAGGTCAAACTGTTTGGCAAAGGCCTTGGCCGCCGTCATGTAGC  >  NC_002947/3543812‑3544108

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