Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,949,108	Δ1 bp	6.7%	coding (426/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,108	0	T	.	6.7%	278.2 / 13.9	89	coding (426/642 nt)	PP_3485	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (41/42);  new base . (2/4);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

CCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTGCGC  >  NC_002947/3948961‑3949254                                                                                                                                                      |                                                                                                                                                     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)          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                                                                                                                                    <  5:238897/149‑1 (MQ=255)                 gTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTa                                                                                                                                    >  4:309667/1‑149 (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                                                                                                                          >  5:391968/1‑149 (MQ=255)                           ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTgg                                                                                                                          >  7:235764/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                                                                                                                               <  5:323791/126‑1 (MQ=255)                                             cgcgTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACtt                                                                                                                               >  6:323791/1‑126 (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                                                                                                      >  3:30655/1‑141 (MQ=255)                                                       cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                                      <  4:30655/141‑1 (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                                                                                                                             <  8:54081/107‑1 (MQ=255)                                                                  cACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGc                                                                                                                             >  7:54081/1‑107 (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                                                                                            >  7:322563/1‑118 (MQ=255)                                                                                       gCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑Agt                                                                                            <  8:322563/118‑1 (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                                     >  2:379822/1‑149 (MQ=255)                                                                                                               aaaaTGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCggg                                     <  6:391968/149‑1 (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                                                                                              <  3:147284/84‑1 (MQ=255)                                                                                                                       tCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑a                                                                                              >  4:147284/1‑84 (MQ=255)                                                                                                                       tCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCaaa                             >  1:373353/1‑149 (MQ=255)                                                                                                                       tCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCaaa                             >  7:52960/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)                                                                                                                                                gTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCt      >  7:229245/1‑149 (MQ=255)                                                                                                                                                 tCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTg     >  7:386978/1‑149 (MQ=255)                                                                                                                                                 tCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGAGCTGCCAAATGCCGCTTCGATCTCGCCAACCTg     >  1:381245/1‑149 (MQ=255)                                                                                                                                                    aCTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTgcgc  <  4:103096/149‑1 (MQ=255)                                                                                                                                                     cTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCt      >  2:26244/1‑144 (MQ=255)                                                                                                                                                     cTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCt      <  1:26244/144‑1 (MQ=255)                                                                                                                                                      tGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGcc           >  6:55709/1‑138 (MQ=255)                                                                                                                                                      tGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGcc           <  5:55709/138‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTGCGC  >  NC_002947/3948961‑3949254

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