Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,341,588	A→C	5.2%	K232T (AAG→ACG)	panB →	3‑methyl‑2‑oxobutanoate hydroxymethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,341,588	0	A	C	5.2%	225.5 / 2.0	72	K232T (AAG→ACG)	panB	3‑methyl‑2‑oxobutanoate hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (36/32);  new base C (2/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CCAGCGAACTGGCGGCGGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCC  >  NC_002947/5341444‑5341736                                                                                                                                                  |                                                                                                                                                        ccAGCGAACTGGCGGCGGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAAc                                                                                                                                                     >  5:140151/1‑149 (MQ=255)    gCGAACTGGCGGCGGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTc                                                                                                                                                  <  3:38429/149‑1 (MQ=255)    gCGAACTGGCGGCGGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTc                                                                                                                                                  <  7:208743/149‑1 (MQ=255)     cGAACTGGCGGCGGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCTTCACGAACTTcg                                                                                                                                                 <  3:216374/149‑2 (MQ=255)              ggcggAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCAGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGc                                                                                                                                         >  8:170554/1‑148 (MQ=255)                     aTCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGa                                                                                                                                 >  6:133844/1‑149 (MQ=255)                       cACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGata                                                                                                                               >  8:230158/1‑149 (MQ=255)                        aCCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGatat                                                                                                                              >  1:146835/1‑149 (MQ=255)                         ccAATTCCGTTGGCATCCCGGTAATCGGCATTGTCGCAGGCAACGCCACAGAAGTCAAGGTACTAGTTCTAAACGACATTAACGGACTGT‑CGCTTAGCTGCAGCGAACAGAAGTTCGTGCAGAACGACATGCTCGGCCAGCCGGagaga                                                                                                                             >  4:174047/1‑145 (MQ=255)                          cAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATc                                                                                                                             >  6:92805/1‑148 (MQ=255)                             tGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTATCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATa                                                                                                                         <  1:72534/149‑1 (MQ=255)                             tGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATa                                                                                                                         <  2:189104/149‑1 (MQ=255)                               ccGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAg                                                                                                                        <  2:163868/148‑1 (MQ=255)                                     gCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCtc                                                                                                                 >  5:311356/1‑149 (MQ=255)                                           cGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCt                                                                                                          >  7:47768/1‑149 (MQ=255)                                              tAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATc                                                                                                                             >  5:35596/1‑128 (MQ=255)                                              tAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATc                                                                                                                             <  6:35596/128‑1 (MQ=255)                                                   ggCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCg                                                                                                   <  8:47768/148‑1 (MQ=255)                                                   atcggtGGCGCTGCCGGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGAGATCCATAGCGCCCTCGTC‑GCTTACGTc                                                                                                    >  2:14587/6‑147 (MQ=255)                                                       gtGGCGCTGCCGGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCtcgtc‑g                                                                                                            >  7:192792/2‑135 (MQ=255)                                                       gtGGCGCTGCCGGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCtcgtc‑g                                                                                                            <  8:192792/134‑1 (MQ=255)                                                          gcgcCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGTAAGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCACCCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGt                                                                                            >  1:267841/1‑149 (MQ=255)                                                            gcCGGCAGCGCCACCGATGGCCAGGTACTGTTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGTTACCGAAGTTCGTGAATAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTcaa                                                                                         <  5:311362/149‑1 (MQ=255)                                                              cGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTcaagc                                                                                       <  6:140151/149‑1 (MQ=255)                                                                     gcCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAg                                                                                >  1:109859/1‑149 (MQ=255)                                                                       cACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGtt                                                                              >  5:20933/1‑149 (MQ=255)                                                                          cGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCAAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTcc                                                                           >  5:274654/1‑149 (MQ=255)                                                                           gATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTccc                                                                          <  5:92805/149‑1 (MQ=255)                                                                           gATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTcc                                                                           >  7:220776/1‑148 (MQ=255)                                                                             tGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAg                                                                        <  6:311356/149‑1 (MQ=255)                                                                               gCCAGGTAGTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGCCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCCAGGCGGTCAAGCAAGTCAGTTTccagg                                                                        >  2:81512/1‑144 (MQ=255)                                                                               gCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGc                                                                      >  2:264401/1‑149 (MQ=255)                                                                                 cAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAg                                                                    >  5:349707/1‑149 (MQ=255)                                                                                  aGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑c                                                                  >  7:276406/1‑149 (MQ=255)                                                                                   ggtactggtactGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑c                                                                  >  4:124848/1‑148 (MQ=255)                                                                                          gtactgCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTCGGCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑Aacac                                                           <  1:254314/148‑1 (MQ=255)                                                                                          gtactgCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACgg                                                         <  2:267841/149‑1 (MQ=255)                                                                                            actgCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGt                                                       >  6:57580/1‑149 (MQ=255)                                                                                             ctgCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGtt                                                      >  1:236980/1‑149 (MQ=255)                                                                                             ctgCACGACATGCTCGGCCTGT‑CGCTGAGCAGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGtt                                                      <  7:170554/149‑1 (MQ=255)                                                                                               gCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGtt                                                      >  3:120445/1‑147 (MQ=255)                                                                                                cACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTcaagcaa                                                                                     <  5:24313/117‑1 (MQ=255)                                                                                                cACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTcaagcaa                                                                                     >  6:24313/1‑117 (MQ=255)                                                                                                 aCGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGt                                                  <  7:230158/149‑1 (MQ=255)                                                                                                   gACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGc                                                <  6:20933/149‑1 (MQ=255)                                                                                                   gACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACAGGATATCCATAGCGCCCTCGTA‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTg                                                 >  2:258503/1‑148 (MQ=255)                                                                                                    aCATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCa                                               <  6:349707/149‑1 (MQ=255)                                                                                                     cATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAgg                                                                       >  2:209918/1‑126 (MQ=255)                                                                                                     cATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAgg                                                                       <  1:209918/126‑1 (MQ=255)                                                                                                     cATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCAt                                              >  3:336177/1‑149 (MQ=255)                                                                                                     cATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCAt                                              <  4:26715/149‑1 (MQ=255)                                                                                                      aTGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATg                                             <  3:124848/149‑1 (MQ=255)                                                                                                        gCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCtcgtc                                                                                                              <  5:96396/85‑1 (MQ=255)                                                                                                        gCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCtcgtc                                                                                                              >  6:96396/1‑85 (MQ=255)                                                                                                             gCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCTGTCAAGCAAGTCAGTTTCCCAGGCAGACG‑AACACGGGTTCAGTGCATGAATACa                                       <  1:14587/149‑1 (MQ=255)                                                                                                               cTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAGTc                                    <  8:220776/149‑1 (MQ=255)                                                                                                                tGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAGTc                                    >  3:136755/1‑148 (MQ=255)                                                                                                                  t‑CGCTGAGCGGCCGTGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAGTCAAg                                 <  5:133844/149‑1 (MQ=255)                                                                                                                  t‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAGTCaa                                  <  2:109859/148‑1 (MQ=255)                                                                                                                    cGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTcaagcaa                                                                                     <  6:264821/98‑1 (MQ=255)                                                                                                                    cGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTcaagcaa                                                                                     >  5:264821/1‑98 (MQ=255)                                                                                                                     gCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGt                                                  >  5:172324/1‑130 (MQ=255)                                                                                                                     gCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGt                                                  <  6:172324/130‑1 (MQ=255)                                                                                                                               cGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGaa                     <  5:57580/149‑1 (MQ=255)                                                                                                                                gggTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAGTCAAGACCGTGCGCGAAc                    >  4:205071/1‑149 (MQ=255)                                                                                                                                  tgcatgAAGTTCTTCACGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAGTc                                    >  5:264152/6‑131 (MQ=255)                                                                                                                                  tgcatgAAGTTCTTCACGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAGTc                                    <  6:264152/126‑1 (MQ=255)                                                                                                                                     atgAAGTTCTTCACGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTgcgc               >  6:214762/3‑149 (MQ=255)                                                                                                                                       gAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAg                                      <  6:292653/124‑1 (MQ=255)                                                                                                                                       gAAGTTCGTGAAGAACTTCAAGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAg                                      >  5:292653/1‑124 (MQ=255)                                                                                                                                        aaGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CGCAACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTgcgcgc             <  1:258503/149‑1 (MQ=255)                                                                                                                                        aaGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGcgcgcg            >  3:71066/1‑149 (MQ=255)                                                                                                                                                  agaaCTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGgccc  <  2:5031/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                        CCAGCGAACTGGCGGCGGAAATCACCAATGCCGTTGGCATCCCGGTAATCGGCATTGGCGCCGGCAGCGCCACCGATGGCCAGGTACTGGTACTGCACGACATGCTCGGCCTGT‑CGCTGAGCGGCCGGGTACCGAAGTTCGTGAAGAACTTCATGCAGGGCCAACCGGATATCCATAGCGCCCTCGTC‑GCTTACGTCGAGGCGGTCAAGCAAGTCAGTTTCCCAGGCAG‑CG‑AACACGGGTTCAGTGCATGAATACAGTCAAGACCGTCCGCGAACTGCGCGCGGCAGTGGCCC  >  NC_002947/5341444‑5341736

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