Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,375,433	C→T	6.6%	A35T (GCC→ACC)	PP_5637 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,375,433	0	C	T	6.6%	177.5 / 2.3	60	A35T (GCC→ACC)	PP_5637	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base T (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

GCCCTGCACCACCACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAA  >  NC_002947/4375285‑4375581                                                                                                                                                       |                                                                                                                                                         gCCCTGCACCACCACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGc                                                                                                                                                          >  6:2587/1‑149 (MQ=255) gCCCTGCACCACCACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGc                                                                                                                                                          >  3:176060/1‑149 (MQ=255)     tGCACCACCACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTgg                                                                                                                                                      >  2:340647/1‑149 (MQ=255)       cacctccacGACTGCGCATTTAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTgg                                                                                                                                                      <  7:139752/147‑1 (MQ=255)       caccaccacGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTgg                                                                                                                                                      <  2:127415/147‑1 (MQ=255)         ccaccacGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCAGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTcg                                                                                                                                                  >  8:94763/1‑149 (MQ=255)            ccacGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGa                                                                                                                                               >  2:359246/1‑149 (MQ=255)             cacGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGAc                                                                                                                                              <  8:217233/149‑1 (MQ=255)              acGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACt                                                                                                                                             >  7:89342/1‑149 (MQ=255)              acGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACt                                                                                                                                             <  7:195805/149‑1 (MQ=255)                gACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCg                                                                                                                                           >  7:21158/1‑149 (MQ=255)                 aCTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑a                                                                                                                                      <  4:176060/149‑1 (MQ=255)                  cTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑Ac                                                                                                                                     >  7:131186/1‑149 (MQ=255)                     cgcATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑Acc                                                                                                                                    <  8:223129/147‑1 (MQ=255)                        atatAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCagag                                                                                                                               <  5:161113/149‑1 (MQ=255)                            aGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTc                                                                                                                           >  1:186410/1‑149 (MQ=255)                             gCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCg                                                                                                                          <  8:21158/149‑1 (MQ=255)                                 gCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCAAGGCCTGGGTCGCGACTc                                                                                                                                            >  5:30836/1‑131 (MQ=255)                                   cTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGa                                                                                                                    <  5:2587/149‑1 (MQ=255)                                    ttCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTATGACCTTCGCATCGTGCAGCCGGGGGATAATCTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGa                                                                                                                    >  1:276023/1‑149 (MQ=255)                                    ttCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGTCTGGGTCGAGTCGCG‑‑‑‑ACCCaggcctgg                                                                                                                           <  7:299302/141‑7 (MQ=255)                                    ttCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGTCTGGGTCGAGTCGCG‑‑‑‑ACCCaggcctgg                                                                                                                           >  8:299302/1‑135 (MQ=255)                                       gACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGTCTGGGTCGAGTCGCG‑‑‑‑ACCCaggcctgg                                                                                                                           >  7:315873/1‑132 (MQ=255)                                       gACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGTCTGGGTCGAGTCGCG‑‑‑‑ACCCaggcctgg                                                                                                                           <  8:315873/138‑7 (MQ=255)                                       gACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑CGCCAGAGGCTCGCTGGGATGGt                                                                                                                >  4:109456/1‑149 (MQ=255)                                        aCCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTc                                                                                                               <  8:89342/149‑1 (MQ=255)                                        aCCTCTGCTACCCCGTTGATGTTGATGCCGAGCGGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATACTGTCTGG‑CCCAGGCGTAGGTTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTc                                                                                                               <  4:107921/149‑1 (MQ=255)                                               ctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAAt                                                                                                         <  6:170923/148‑1 (MQ=255)                                                  ccccTTTGATTGTGATGCGTTGCGGGGCCTTGCGGT‑TTACCTTCTCATCGTGCAGCCGGGGGATAATGTGGGG‑CCCAGGCGTTGGGTTAAACACCCAGGCCTTGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAAtgtg                                                                                                      <  8:131186/148‑1 (MQ=255)                                                    ccGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGAc                                                                                                   >  7:302769/1‑149 (MQ=255)                                                          atgttgatgCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgc                                                                                             >  4:135923/1‑149 (MQ=255)                                                             ttgatgCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgcgcg                                                                                          <  1:148182/149‑1 (MQ=255)                                                                atgCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAg                                                                                       <  4:43746/149‑1 (MQ=255)                                                                  gCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgcgc                                                                                           >  1:278776/1‑143 (MQ=255)                                                                  gCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgcgc                                                                                           <  2:278776/143‑1 (MQ=255)                                                                  gCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAggg                                                                                     >  7:212732/1‑149 (MQ=255)                                                                   ccGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTTGGGCCCCCGGCGTAGGGTTAAACACCCAGGCATGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAggg                                                                                     >  1:96417/1‑149 (MQ=255)                                                                   ccGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGTGGTGGAGAGAGGCTCGCTGTGATGGTCGACAATGTGGACGTCcgcgcg                                                                                          >  1:308123/1‑147 (MQ=255)                                                                         tGGGCCTTGCCGT‑TGACATTGGCATGGTGCTGCGGGGGTATAATGTCGGG‑CCCAGGCTTTGGTTAAAACACCTAGGCGTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGAc                                                                              <  2:186410/149‑1 (MQ=255)                                                                          gggCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACTCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACt                                                                             >  8:47563/1‑149 (MQ=255)                                                                            gCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACttt                                                                           >  1:28975/1‑149 (MQ=255)                                                                             ccATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTc                                                                          <  1:340647/149‑1 (MQ=255)                                                                                  ccGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACTCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTa                                                                      <  7:47563/148‑1 (MQ=255)                                                                                   cGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAgg                                                                    <  8:302769/149‑1 (MQ=255)                                                                                         aCCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACttt                                                                           <  5:90962/137‑1 (MQ=255)                                                                                         aCCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACttt                                                                           >  6:90962/1‑137 (MQ=255)                                                                                              cGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAg                                                           <  8:212732/148‑1 (MQ=255)                                                                                               gCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAgg                                                          <  2:96417/148‑1 (MQ=255)                                                                                                 aTCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGc                                                        >  6:6707/1‑148 (MQ=255)                                                                                                       cAGCCGGGGGATAATGTCGGG‑CCCAGGCGTTGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCa                                                 >  1:30712/1‑149 (MQ=255)                                                                                                                     tGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCgcgg                                   >  2:34037/1‑149 (MQ=255)                                                                                                                        cGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGt                                <  7:94763/149‑1 (MQ=255)                                                                                                                           g‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGa                             <  6:320910/149‑1 (MQ=255)                                                                                                                              ccAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTTGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAAc                           <  5:178601/149‑1 (MQ=255)                                                                                                                                aGGCGTAGTGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACga                         <  2:308123/149‑1 (MQ=255)                                                                                                                                    gTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaag                     <  5:6707/149‑1 (MQ=255)                                                                                                                                       gggTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaagg                    >  4:24445/1‑147 (MQ=255)                                                                                                                                        tgTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaaggaag                 <  2:276023/148‑1 (MQ=255)                                                                                                                                        ggTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaaggaag                 >  6:13023/1‑149 (MQ=255)                                                                                                                                        ggTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACCTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaaggaag                 >  4:231533/1‑149 (MQ=255)                                                                                                                                         gTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACgag                        <  6:201308/141‑1 (MQ=255)                                                                                                                                         gTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACgag                        >  5:201308/1‑141 (MQ=255)                                                                                                                                                       ccTGGGTCGCGACTCG‑‑‑‑ACCGAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCaa  <  3:109456/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                         GCCCTGCACCACCACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGT‑TGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGG‑CCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCG‑‑‑‑ACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAA  >  NC_002947/4375285‑4375581

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑