Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	4,023,563	T→C	L4P (CTC→CCC)	rnpA →	ribonuclease P

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	4,023,563	0	T	C	97.4%	122.0 / ‑4.3	39	L4P (CTC→CCC)	rnpA	ribonuclease P
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (22/16);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

CAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCT  >  CP006584/4023281‑4023854                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                      cAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGa                                                                                                                                                                                                                                                                                        >  1:174037/1‑298 (MQ=255)                   atcttttgttaCCCATGGCTTCGGTTTCCCTCTTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGGATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACGTAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑cc                                                                                                                                                                                                                                                                    <  2:29440/288‑1 (MQ=255)                                 aTGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACAt                                                                                                                                                                                                                                                      <  1:295769/298‑1 (MQ=255)                                           ggtccagacttttATTTTTAGGTGGCAATCAAACAAAGCCATTGAATTTCTTAAAATTTACGTAGAACTCCACATTTCACGTCGTTTTCATCCGTGTGTACTTAACCGCAGCCGTCCTCACTGCTTCCGGGCTCGTATGGCTACTAAAAATGTTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCTGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCca                                                                                                                                                                                                                                            <  2:40634/289‑1 (MQ=255)                                              ccATCGTTTCATTTTTCGGCGGATATCCAATACAGCCATTGAATCTATTCAAGTTTAGGTCGAAATCGCCATGAACCGCCCTTTTCAACCGTCTGTACTCAAGCGCAACCATTCTCACGGCTTCTGTGCCAGTAGGACTACTAAACATGGTCTCAAGGTTACGGCACGGCCTCGTGCCTAAGGACGGG‑CGCGCCTCAACGCTTCTAAGTAATAAACCTAACCCCCGAGTTGTTCGGCCCCCCTTTCCCAGCGAATTCGCCTTGTTCCACT‑‑CCAGTCCATTCACTTTcctgttacacca                                                                                                                                                                                                                                         >  2:56455/1‑287 (MQ=255)                                              ccATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCcagca                                                                                                                                                                                                                                         <  1:141676/298‑1 (MQ=255)                                              ccATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGCCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCcagca                                                                                                                                                                                                                                         <  1:295945/298‑1 (MQ=255)                                                          ttttCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGc                                                                                                                                                                                                                             <  1:118316/298‑1 (MQ=255)                                                                    atatCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCAcgccg                                                                                                                                                                                                                   >  1:133350/1‑298 (MQ=255)                                                                         cAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTCAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACTCCCCAGAC‑ATTCACATTCGTCTTCACTCCGCCACAACGGGCTGGCACCCCACAAAt                                                                                                                                                                                                              >  2:217157/1‑298 (MQ=255)                                                                              aaGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTA‑CACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCAcccccacaattccca                                                                                                                                                                                                         >  2:127358/1‑287 (MQ=255)                                                                                        aTTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGgccg                                                                                                                                                                                               <  1:6054/298‑1 (MQ=255)                                                                                             ttCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGa                                                                                                                                                                                          >  1:101807/1‑298 (MQ=255)                                                                                                                      gcaaCGCACTTTCCAGCCGTCTGTACGAAACCGCATCCGTTCTACCGTCTTCCTGTCTCGTATGGTAAATAAAAATGGTCTTCAGGTTCTGGCACTGCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCgg                                                                                                                                                                 <  2:132431/296‑1 (MQ=255)                                                                                                                               ttttCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGt                                                                                                                                                        <  1:127358/298‑1 (MQ=255)                                                                                                                               ttttCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGt                                                                                                                                                        >  1:308977/1‑298 (MQ=255)                                                                                                                                    aaCCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCa                                                                                                                                                   >  1:293860/1‑298 (MQ=255)                                                                                                                                    aaCCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTAGTCAGGTTCTCGCACGTCGTCGTGCTAAAGGCCGCG‑CCCGTCTGACCGTTTCTAAGTAATAAATCTAACCCCTGAGTTGTTAAGCCCGCACTTCCCAGGGAGTTACACTTGTT‑AACT‑CCCAGTCAATTCACATCCGTCTTCTCGCAGCCACACCGGGCTGGCAAGCCCAAAATCACCTTTCCACGCTCCCCGAATTCTCTGGAGAGTCCCCTGATCGGTccacacatcccca                                                                                                                                                   >  2:223724/1‑286 (MQ=255)                                                                                                                                                      gcAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAACGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑CACT‑CCCAGTCAATTCACATTCGTCTTCCAGCCACCACAACGGGCCGGCCCGCCGCCACTTACCATTCCAGGCCGCATGAACTCCCTGGGGCAACCCAGTATCGTGCTTACAGTCCCCAATACAAAAGTTccagccc                                                                                                                                 >  2:110443/1‑292 (MQ=255)                                                                                                                                                            gTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTc                                                                                                                                                                                      <  2:45481/239‑1 (MQ=255)                                                                                                                                                            gTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTc                                                                                                                                                                                      >  1:45481/1‑239 (MQ=255)                                                                                                                                                                        tCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGa                                                                                                               <  1:294676/298‑1 (MQ=255)                                                                                                                                                                               gttctgttcctctaAAATTGGCGTTCAGGTTAGTGCACGTCGTCGGGC‑CAAGTGCGCGTCTTGTTTGACCGTTTCTCAGTAATAAAGTGATCCCCTGGGTGTTTAAGCCCGCATTTCCCAGGTGGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGTCTGGCACGCCGCAAATTACCATTCCCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACg                                                                                                        <  2:317023/287‑1 (MQ=255)                                                                                                                                                                                     ggCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCAACAACGGGCTGGCACGCCGAAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCTGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGtagatc                                                                                                  >  1:317023/1‑293 (MQ=255)                                                                                                                                                                                               aaTGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGc                                                                                        >  1:270637/1‑298 (MQ=255)                                                                                                                                                                                                  ggtcgtgtgtgtgtGGCACTTCGTCGGTGTAACGGACGCT‑CTCGTCTGACCGTTTCTAAGGAATAAAGCTAGCCCCGGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCAACACGGGCTGGCACGCCGCTAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTc                                                                                     <  2:18721/285‑1 (MQ=255)                                                                                                                                                                                                         aGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTgc                                                                              >  1:50537/1‑298 (MQ=255)                                                                                                                                                                                                                      tcgtcgtGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCATCAGCCACAACGGGCTGGCACGCCGCAAATTCCCATTCTCGGCCGCCTGAATTCGCTGGGCCACCCCCGTACCGGTCTTACAGTCGCCACGAAAAACGTTCGACGCGCCCACGACCCCAATCGGACTAAACTTCCGACGCGTTACAGCTTCACGCCGCGCCAACCTGaccc                                                                 >  2:336839/1‑295 (MQ=255)                                                                                                                                                                                                                                           cTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCgtgg                                             >  1:231072/1‑298 (MQ=255)                                                                                                                                                                                                                                            tCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCgtggt                                            >  1:133395/1‑298 (MQ=255)                                                                                                                                                                                                                                                 tGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCgtggt                                            >  2:109356/1‑293 (MQ=255)                                                                                                                                                                                                                                                 tGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCgtggt                                            <  1:109356/293‑1 (MQ=255)                                                                                                                                                                                                                                                            aaGTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCGCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAggg                            >  1:335625/1‑298 (MQ=255)                                                                                                                                                                                                                                                             attaataaAGTTACCCCCTGAGTGGTTCACCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCCCCTTCGTCTTCCCGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAgggg                           <  2:50537/296‑1 (MQ=255)                                                                                                                                                                                                                                                              gTAATAAAGCTAACCCCTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGt                          <  1:250350/298‑1 (MQ=255)                                                                                                                                                                                                                                                                aataaAGCTAACCCCTGAGTGGTTAAGCCCGCCTTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGCCTGGCACGACGCAAATTACCATTCTCCGCCGCCTGAATCCCCCGGGGCACCCCCGTACCGGTGTTACATGCCCCAAGAAAAAAGTTCGACGCGCCCATGAACGCACTAGGATTAAACGTCAGACGCGTAAAAGCTACCGCCTGCGCCAACACGAACACCCGGCTATCGACTCTCTGGTGGTGGCAAAAAAAGGGGtag                        >  2:251725/1‑296 (MQ=255)                                                                                                                                                                                                                                                                             ccTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATaa           <  1:217157/298‑1 (MQ=255)                                                                                                                                                                                                                                                                              cTGAGTGGTTAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAAc          >  1:224127/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                      ttAAGCCCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTCTCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGctct  >  1:174524/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                      CAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCG‑CTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTT‑AACT‑CCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCT  >  CP006584/4023281‑4023854

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑