Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,199,963:1	+C	11.0%	coding (829/1278 nt)	PP_1944 →	aminomethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,199,963	1	.	C	11.0%	256.8 / 25.1	73	coding (829/1278 nt)	PP_1944	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base . (33/32);  new base C (4/4);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

ACTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTG  >  NC_002947/2199819‑2200109                                                                                                                                                    |                                                                                                                                                   aCTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑Tc                                                                                                                                                  >  3:181390/1‑147 (MQ=255) aCTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGCGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGATATCGCAGCCCAG‑Tc                                                                                                                                                  >  6:250759/1‑147 (MQ=255)  cTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGa                                                                                                                                                <  6:261488/148‑1 (MQ=255)     aCTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGTATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgct                                                                                                                                             <  1:55068/148‑1 (MQ=255)       tACATGCCGGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctggg                                                                                                                                          <  3:44084/149‑1 (MQ=255)          aTGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGCTCGa‑ctgggctg                                                                                                                                       <  1:93286/149‑1 (MQ=255)             cccGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCtgggctgggc                                                                                                                                    >  4:257067/1‑149 (MQ=255)              ccGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCtgggctgggc                                                                                                                                    >  8:96263/1‑148 (MQ=255)               cgcgATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTCCTCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCCCAGCCCAGCTCGa‑ctgggctgcga                                                                                                                                    <  1:91944/149‑4 (MQ=255)               cgcgATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGcc                                                                                                                                   >  3:300494/1‑148 (MQ=255)                         gACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCa                                                                                                                         <  3:257067/148‑1 (MQ=255)                            ccAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGtt                                                                                                                     >  3:153936/1‑149 (MQ=255)                             cAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTc                                                                                                                    >  8:207000/1‑149 (MQ=255)                              agagATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCCCAGCCCAGCTCGa‑ctgggctg                                                                                                                                       >  1:148793/1‑129 (MQ=255)                              agagATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCCCAGCCCAGCTCGa‑ctgggctg                                                                                                                                       <  2:148793/129‑1 (MQ=255)                                agaTGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGa                                                                                                                  >  4:112279/1‑148 (MQ=255)                                  aTGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATcg                                                                                                               >  1:263128/1‑148 (MQ=255)                                  aTGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCa                                                                                                               >  7:126747/1‑149 (MQ=255)                                      cGGTTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGTAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGAc                                                                                                           <  5:67210/149‑1 (MQ=255)                                                cAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATc                                                                                                     <  4:139170/145‑1 (MQ=255)                                                cAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATc                                                                                                     >  3:139170/1‑145 (MQ=255)                                                cAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTc                                                                                                 >  6:34255/1‑149 (MQ=255)                                                  aGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGt                                                                                               <  2:263128/149‑1 (MQ=255)                                                       ttCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCcgc                                                                                          >  4:122446/1‑149 (MQ=255)                                                          cgcgcgGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGa                                                                                                            <  1:270820/128‑1 (MQ=255)                                                          cgcgcgGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGa                                                                                                            >  2:270820/1‑128 (MQ=255)                                                          cgcgcgGGTATGCCAAGTTTTGCCTAGACTTTCAATGTT‑‑GCGCGGAGCTTAGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGATCGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGtt                                                                                       >  7:135243/1‑149 (MQ=255)                                                           gcgcgGGTATGCCAAGTTTTGCCGCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTg                                                                                      <  5:260812/149‑1 (MQ=255)                                                            cgcgGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTgg                                                                                     <  8:126747/149‑1 (MQ=255)                                                                 gTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGaaaaa                                                                                <  3:112279/149‑1 (MQ=255)                                                                     gCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctgg                                                                                                                                           <  4:121275/86‑1 (MQ=255)                                                                     gCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctgg                                                                                                                                           >  3:121275/1‑86 (MQ=255)                                                                     gCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAagacg                                                                            >  4:104515/1‑147 (MQ=255)                                                                     gCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAgg                                                                            >  7:144823/1‑149 (MQ=255)                                                                        aaGGTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTg                                                                         <  4:181390/149‑1 (MQ=255)                                                                         aGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTAGCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGc                                                                        <  5:250759/149‑1 (MQ=255)                                                                                ccTCGACTTTCAATGTT‑‑GCGGGGGGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGGCGCGTTGGAAAAAGAGGTTGTAGAGCCt                                                                 >  7:194820/1‑149 (MQ=255)                                                                                  tCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTc                                                                >  2:14241/1‑148 (MQ=255)                                                                                  tCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTc                                                                >  3:71355/1‑148 (MQ=255)                                                                                      cTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACgg                                                           <  7:96263/149‑1 (MQ=255)                                                                                         tCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGt                                                         <  4:300494/148‑1 (MQ=255)                                                                                          cAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACttt                                                                                                        <  6:223338/100‑1 (MQ=255)                                                                                            aTGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATg                                                     >  5:29244/1‑149 (MQ=255)                                                                                             tggt‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATgc                                                    <  7:207000/146‑1 (MQ=255)                                                                                               tt‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATgcgc                                                  <  8:144823/149‑1 (MQ=255)                                                                                                   gCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCa                                                 >  3:198379/1‑148 (MQ=255)                                                                                                     ggggAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCt                                              >  3:78642/1‑149 (MQ=255)                                                                                                        gAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGAGTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTgg                                           <  3:104515/149‑1 (MQ=255)                                                                                                           cTTTCAGCCGAATCAGGGAACCGATTCGCATCCCAGCCCAG‑TCGACCTGGGCGGGGCCCAAACCATCAAGTTCGAACACGCCTATATCCGGCGTGCCGCGTTCCAAAACGCGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGttt                                        <  8:135243/149‑1 (MQ=255)                                                                                                            tttGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                       >  4:263412/1‑149 (MQ=255)                                                                                                              tGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                       <  4:71355/147‑1 (MQ=255)                                                                                                                 gTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTg                                                        <  4:237286/127‑1 (MQ=255)                                                                                                                 gTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTg                                                        >  3:237286/1‑127 (MQ=255)                                                                                                                  tCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTgg                                                                                     <  7:36549/97‑1 (MQ=255)                                                                                                                  tCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTgg                                                                                     >  8:36549/1‑97 (MQ=255)                                                                                                                  tCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATgg                                 <  3:122446/149‑1 (MQ=255)                                                                                                                    gAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtg                               <  6:29244/149‑1 (MQ=255)                                                                                                                           gTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGt                        >  4:244600/1‑149 (MQ=255)                                                                                                                           gTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGt                        >  4:240338/1‑149 (MQ=255)                                                                                                                            tAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGtt                       >  3:214315/1‑149 (MQ=255)                                                                                                                                cGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAg                   <  1:14241/149‑1 (MQ=255)                                                                                                                                   ttGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGAt                 <  4:214315/148‑1 (MQ=255)                                                                                                                                      gTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtgatga                           >  1:60303/1‑135 (MQ=255)                                                                                                                                      gTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtgatga                           <  2:60303/135‑1 (MQ=255)                                                                                                                                      gTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGtt              <  4:78642/148‑1 (MQ=255)                                                                                                                                       tATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGTAGAGCCTCGACGGGTGATGCGCACCTTGTTTTGGAATGGTGATGACGTTGCAGATGTTCa            <  8:194820/149‑1 (MQ=255)                                                                                                                                       tATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                       >  4:278893/1‑122 (MQ=255)                                                                                                                                       tATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                       <  3:278893/122‑1 (MQ=255)                                                                                                                                         gcccAGCCCAGCTCGA‑CTGGGCTGCGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGAGTTGGAAAAAGAGGTTGCAGAGACACGACGGGTGATGCGCACCTTGGTTTGGAATCGTGATGACGTTGCAGATGGTCaag          >  6:163569/4‑147 (MQ=255)                                                                                                                                          cccAGCCCAGCTCGA‑CTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtg                               <  7:178828/125‑1 (MQ=255)                                                                                                                                          cccAGCCCAGCTCGA‑CTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtg                               >  8:178828/3‑127 (MQ=255)                                                                                                                                          cccAGCCCAGCTCGA‑CTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGc         >  3:34378/3‑149 (MQ=255)                                                                                                                                                 cAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATCCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTg  >  3:37654/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   ACTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTT‑‑GCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTG  >  NC_002947/2199819‑2200109

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