Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,949,101:1	+A	6.9%	coding (433/642 nt)	PP_3485 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,949,101	1	.	A	6.9%	350.5 / 17.5	87	coding (433/642 nt)	PP_3485	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (39/42);  new base A (2/4);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

GCCTCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGC  >  NC_002947/3948956‑3949244                                                                                                                                                    |                                                                                                                                                  gcctCCCGCTCCCCGATGAGCGTTGCCCGACACAAGGCCAGGCAGTGTTCGCGTATGACCAATGGATTATCACTTTCGTTCTGCGATCCTAGCAGTCTGTCCTGGTAGGCATCAAAAAATGACTCCACTA‑CTCCCTCGCACGGCGC‑Gt                                                                                                                                                 >  4:295052/1‑148 (MQ=255)   ctCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGa                                                                                                                                              <  7:100495/149‑1 (MQ=255)   ctCCCGATGCACGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCCTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCAGTCCTGGTTGGCATCAAAAAATCACTCCACAT‑CGCCCTCGCTCGGCGC‑GTCg                                                                                                                                               >  6:356631/1‑148 (MQ=255)      ccGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACt                                                                                                                                            <  2:317029/148‑1 (MQ=255)      ccGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTg                                                                                                                                           <  3:85495/149‑1 (MQ=255)          tcgccGATGAGCGTTGCCCGAGCCAAGGCCTGGCCGTGTTCGCGTATGACCAATTGTTTATCAGTTTGGTTCTGCTGTCCTTGCAGTCTGTCCTTGTTGGCATCAAAAAATTCCTCGACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCt                                                                                                                                       <  5:356631/146‑1 (MQ=255)          tGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCt                                                                                                                                       >  1:9559/1‑149 (MQ=255)            cccGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGcgc‑g                                                                                                                                                  >  3:55389/1‑136 (MQ=255)            cccGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGcgc‑g                                                                                                                                                  <  4:55389/136‑1 (MQ=255)               gATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAg                                                                                                                                    >  2:53611/1‑147 (MQ=255)                 tGAGCGTTGCCCGAGGCAAGGCCAGGCCGTGGTAGCGAATGACCAATGCATTATCACTTTGGTTCTGCGCTCCTTGCTGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGcgagg                                                                                                                               <  7:211468/149‑5 (MQ=37)                 tGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑Agg                                                                                                                               <  4:38650/149‑1 (MQ=255)                 tGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGcgagg                                                                                                                               <  8:346036/149‑5 (MQ=38)                  gatcGTTGGCCGAGCCAAGTCCAGGCCGTGTTCGCTTATGACCAATGGTTTATCACTTTGGTTCTGCTCTCCTTGCAGTCTGTCCTGGTTTGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCTCTGAGC‑Agg                                                                                                                               <  8:177066/145‑1 (MQ=255)                   aGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGcgaggg                                                                                                                              <  1:239362/148‑6 (MQ=38)                   aGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGcgaggg                                                                                                                              <  5:69851/148‑6 (MQ=38)                      gTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTa                                                                                                                          >  4:309667/1‑149 (MQ=255)                      gTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTa                                                                                                                          <  5:238897/149‑1 (MQ=255)                       ttGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTac                                                                                                                         <  8:274042/149‑1 (MQ=255)                         gCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTaca                                                                                                                        <  3:309667/148‑1 (MQ=255)                           ccGAGCCAATGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCTCTGAGC‑AGGGCGTACACtt                                                                                                                     <  3:327923/149‑1 (MQ=255)                            cGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTg                                                                                                                    <  2:149197/149‑1 (MQ=255)                               gCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTg                                                                                                                 <  2:129124/149‑1 (MQ=255)                                ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAGAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTgg                                                                                                                >  8:243970/1‑149 (MQ=255)                                ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTgg                                                                                                                >  3:256370/1‑149 (MQ=255)                                ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTgg                                                                                                                >  7:235764/1‑149 (MQ=255)                                ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTgg                                                                                                                >  5:391968/1‑149 (MQ=255)                                 cAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTggg                                                                                                               <  2:314704/149‑1 (MQ=255)                                   aggccaggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGgcc                                                                                                             <  6:359469/149‑1 (MQ=255)                                     gccaggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGgccgc                                                                                                           >  2:404153/1‑149 (MQ=255)                                       caggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCa                                                                                                         >  3:107302/1‑149 (MQ=255)                                         ggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCAtc                                                                                                       <  8:235764/149‑1 (MQ=255)                                          gccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCAtct                                                                                                      >  5:180027/1‑149 (MQ=255)                                           ccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCAtct                                                                                                      >  7:41235/1‑148 (MQ=255)                                                  cgcgTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACtt                                                                                                                     >  6:323791/1‑126 (MQ=255)                                                  cgcgTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACtt                                                                                                                     <  5:323791/126‑1 (MQ=255)                                                    cgTATGACCAATTGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                            <  2:167531/149‑1 (MQ=255)                                                     gTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgtg                                                                                           <  4:271718/149‑1 (MQ=255)                                                        tGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTAt                                                                                        <  2:64393/149‑1 (MQ=255)                                                          aCCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTAtt                                                                                       <  6:214512/148‑1 (MQ=255)                                                          aCCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTAtt                                                                                       <  4:256370/148‑1 (MQ=255)                                                            cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                            <  4:30655/141‑1 (MQ=255)                                                            cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                            >  3:30655/1‑141 (MQ=255)                                                            cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑Ag                                                                                   >  2:313355/1‑149 (MQ=255)                                                            cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCACTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTA‑Ag                                                                                   >  5:365527/1‑149 (MQ=255)                                                                    tATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTgc                                                                           >  2:319788/1‑149 (MQ=255)                                                                      tCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCAt                                                                                                   <  7:20825/124‑1 (MQ=255)                                                                      tCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCAt                                                                                                   >  8:20825/1‑124 (MQ=255)                                                                       cACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGc                                                                                                                   >  7:54081/1‑107 (MQ=255)                                                                       cACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGc                                                                                                                   <  8:54081/107‑1 (MQ=255)                                                                        aCTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGcc                                                                        <  8:41235/148‑1 (MQ=255)                                                                         cTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTACGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGccc                                                                       >  1:319616/1‑149 (MQ=255)                                                                         cTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCt                                                                      <  2:9559/149‑1 (MQ=255)                                                                              gTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCg                                                                  <  1:53611/148‑1 (MQ=255)                                                                                   gcgcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGccc                                                            >  3:75893/1‑149 (MQ=255)                                                                                     gcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGc                                                          >  8:144800/1‑149 (MQ=255)                                                                                         cTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCctc                                                       <  3:295052/148‑1 (MQ=255)                                                                                          ttGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTg                                                     >  5:179818/1‑149 (MQ=255)                                                                                            gccgTCGGTCCTGGTGGGCATCAAAAAATGCCTCCACTG‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGTGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGc                                                    <  8:158721/145‑1 (MQ=255)                                                                                            gCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑Agt                                                                                  <  8:322563/118‑1 (MQ=255)                                                                                            gCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑Agt                                                                                  >  7:322563/1‑118 (MQ=255)                                                                                            gCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGc                                                    <  7:144800/148‑1 (MQ=255)                                                                                             cAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGcc                                                   <  1:404153/148‑1 (MQ=255)                                                                                                    tCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGcc                                            <  4:75893/148‑1 (MQ=255)                                                                                                            ggCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGA‑AGCGCGTCCAATTGCTGGGCCGCCATCTCATCCCGGGTGTATTC‑AGTGTACTGCTCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAAt                                    >  1:268122/1‑148 (MQ=255)                                                                                                             gCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATaa                                  <  2:319616/149‑1 (MQ=255)                                                                                                             gCATCAAAAAATGCATGCACTT‑GGCCCTCGCTCGGGGG‑GGGGACTGCGCTTAGC‑AGGGCTTACACTTGCTGGGCAGCCATCTCATGCCTGGTGTAATC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATa                                   <  5:349401/148‑1 (MQ=255)                                                                                                                tCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCTTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGc                               >  4:217011/1‑149 (MQ=255)                                                                                                                tCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGc                               >  8:174697/1‑149 (MQ=255)                                                                                                                  aaaaaaTGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGcc                              <  6:179818/148‑1 (MQ=255)                                                                                                                    aaaaTGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCggg                           <  6:391968/149‑1 (MQ=255)                                                                                                                    aaaaTGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCggg                           >  2:379822/1‑149 (MQ=255)                                                                                                                       aTGCCTCCAGTT‑CGCCCTCGCTCGGCGG‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTTTTCTAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCt                         <  1:319788/149‑1 (MQ=255)                                                                                                                        tGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTAc                                                                              <  6:155023/94‑1 (MQ=255)                                                                                                                        tGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTAc                                                                              >  5:155023/1‑94 (MQ=255)                                                                                                                          ccTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGcc                      >  3:103096/1‑148 (MQ=255)                                                                                                                            tCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑a                                                                                    >  4:147284/1‑84 (MQ=255)                                                                                                                            tCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑a                                                                                    <  3:147284/84‑1 (MQ=255)                                                                                                                            tCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCaaa                   >  7:52960/1‑149 (MQ=255)                                                                                                                            tCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCaaa                   >  1:373353/1‑149 (MQ=255)                                                                                                                             ccACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAAt                  >  2:216529/1‑149 (MQ=255)                                                                                                                              cACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATg                 >  4:196683/1‑149 (MQ=255)                                                                                                                              cACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCACTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGTCTGCCAAATg                 >  6:207848/1‑149 (MQ=255)                                                                                                                                    cGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGc             <  4:107302/148‑1 (MQ=255)                                                                                                                                     gCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCt            <  2:373353/148‑1 (MQ=255)                                                                                                                                       ccTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCATTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTc          >  2:73221/1‑148 (MQ=255)                                                                                                                                        cctgctcaGCGCAGTCGAC‑GCGCCGAGCGAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTc          >  1:56772/4‑148 (MQ=38)                                                                                                                                          tgctcaGCGCAGTCGAC‑GCGCCGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGAtc      >  6:388222/2‑149 (MQ=255)                                                                                                                                               ggCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGc  >  5:388402/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCCTCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGC  >  NC_002947/3948956‑3949244

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