Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,857,535	T→A	5.1%	A23A (GCT→GCA)	PP_3407 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,857,535	0	T	A	5.1%	240.1 / 4.0	79	A23A (GCT→GCA)	PP_3407	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base A (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

CAGAAACTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAGGGCCGGACAATGCCGGTGCTGC  >  NC_002947/3857387‑3857683                                                                                                                                                     |                                                                                                                                                      cAGAAACTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACctgctg                                                                                                                                                      >  1:931289/1‑150 (MQ=255)    aaaCTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCCGGCCTGAGCTTGTCGTACCTGCTGTT‑c                                                                                                                                                  <  1:913791/150‑1 (MQ=255)    aaaCTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑c                                                                                                                                                  <  2:1508689/150‑1 (MQ=255)      aCTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGa                                                                                                                                                <  1:777292/150‑1 (MQ=255)       cTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGAc                                                                                                                                               <  1:993350/150‑1 (MQ=255)       cTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCCACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGAc                                                                                                                                               <  1:1365257/150‑1 (MQ=255)        tGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACa                                                                                                                                              <  1:1703032/150‑1 (MQ=255)        tGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCATGGGCCTGTCGTACCTGCTGTT‑CGACa                                                                                                                                              <  1:1202718/150‑1 (MQ=255)           gggCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGcc                                                                                                                                           <  1:622327/150‑1 (MQ=255)           gggCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGcc                                                                                                                                           <  1:1670888/150‑1 (MQ=255)            tccggccctgctgaagggcctgggacgcaccggtttcccgtttcaccccctgaggggcacctgaacgcgaggcataaggagagcatggggcgggcttggcgtcaatgtcatacagcgtgcGAGCCTGTCGTACCTGCTGTT‑CGACAGcc                                                                                                                                           <  1:468408‑M1/29‑1 (MQ=255)            tccggccctgctgaagggcctgggacgcaccggtttcccgtttcaccccctgaggggcacctgaacgcgaggcataaggagagcatggggcgggcttggcgtcaatgtcatacagcgtgcGAGCCTGTCGTACCTGCTGTT‑CGACAGcc                                                                                                                                           >  2:468408‑M1/121‑149 (MQ=255)            ggCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCt                                                                                                                                          >  1:1625807/1‑150 (MQ=255)              ccGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCg                                                                                                                                        >  2:12431/1‑150 (MQ=255)               cGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGc                                                                                                                                       >  1:913758/1‑150 (MQ=255)                   gcgcAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGc                                                                                                                                   >  2:1430579/1‑150 (MQ=255)                    cgcAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCt                                                                                                                                  <  2:116118/150‑1 (MQ=255)                    cgcAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCGGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCt                                                                                                                                  <  1:1314603/150‑1 (MQ=255)                      cAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGt                                                                                                                                <  2:652405/150‑1 (MQ=255)                        ggATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAGGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTAt                                                                                                                              >  1:1275624/1‑150 (MQ=255)                            gAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACa                                                                                                                          <  1:1430579/150‑1 (MQ=255)                                 ggTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGCAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCGGTATGACATTGAc                                                                                                                     >  2:1032020/1‑150 (MQ=255)                                 ggTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGAc                                                                                                                     >  1:44860/1‑150 (MQ=255)                                  gTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACg                                                                                                                    <  2:609389/150‑1 (MQ=255)                                     cgcTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCa                                                                                                                 >  1:1776890/1‑150 (MQ=255)                                        tACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCaa                                                                                                                >  1:1350298/1‑148 (MQ=255)                                        tACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCaa                                                                                                                <  2:1350298/148‑1 (MQ=255)                                                       gACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCtctc                                                                                               >  1:1207999/1‑150 (MQ=255)                                                       gACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCtctc                                                                                               >  1:1264705/1‑150 (MQ=255)                                                       gACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTCCTGTT‑CGACAGCCTCGCACGCTGTATCACATTGACTGAAAGCCCGCCCCATGCtctc                                                                                               >  2:344448/1‑150 (MQ=255)                                                        aCAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑c                                                                                                                                                  <  1:1193644/98‑1 (MQ=255)                                                        aCAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑c                                                                                                                                                  >  2:1193644/1‑98 (MQ=255)                                                            gCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCACCCCATGCTCTCCTTAt                                                                                          <  2:77869/150‑1 (MQ=255)                                                             cTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATg                                                                                         >  1:1007754/1‑150 (MQ=255)                                                                 cggcATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTc                                                                                     <  2:1007754/150‑1 (MQ=255)                                                                        cGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCa                                                                              <  1:866741/150‑1 (MQ=255)                                                                         gCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTGTGGCATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAg                                                                             <  2:739195/150‑1 (MQ=255)                                                                         gCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAg                                                                             <  1:1417812/150‑1 (MQ=255)                                                                          cccGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAgg                                                                            >  1:32676/1‑150 (MQ=255)                                                                          cccGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAgg                                                                            >  1:1695577/1‑150 (MQ=255)                                                                           ccGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGCAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGt                                                                           <  2:1176102/150‑1 (MQ=255)                                                                                ggTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGcccc                                                                      >  2:580250/1‑150 (MQ=255)                                                                                 gTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCt                                                                     >  2:131596/1‑150 (MQ=255)                                                                                 gTGATCGCCCGGTTTCTGCACTTCTGCCTGCTAGAAGAAGACTTCATTCCTGAGCCTGTCGTACCTGCTGGC‑CGACAGCCTCGAACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCt                                                                     >  1:1077844/1‑150 (MQ=255)                                                                                    aTCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGACGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAg                                                                  <  2:1335398/150‑1 (MQ=255)                                                                                     tagcaccgTTGCTGCGCTTATGCCTGATGGAAGAAGACTTCAGGTCTGCGCATTTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTTAgg                                                                 <  2:447814/143‑1 (MQ=255)                                                                                     tCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAgg                                                                 <  1:1797103/150‑1 (MQ=255)                                                                                     tCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAgg                                                                 >  1:61197/1‑150 (MQ=255)                                                                                     ataaggagagcatggggcgggcttggcgtcaatgtcatacagcgtgcGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTc                                                                                                                                         <  1:1410699‑M1/31‑1 (MQ=255)                                                                                     ataaggagagcatggggcgggcttggcgtcaatgtcatacagcgtgcGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTc                                                                                                                                         >  2:1410699‑M1/48‑78 (MQ=255)                                                                                          cGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTg                                                            >  1:464082/1‑150 (MQ=255)                                                                                          cGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTg                                                            >  2:1803393/1‑150 (MQ=255)                                                                                                cTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACgg                                                      <  1:1032020/150‑1 (MQ=255)                                                                                                cTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACgg                                                      <  1:580250/150‑1 (MQ=255)                                                                                                 tGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACggg                                                     <  2:913758/150‑1 (MQ=255)                                                                                                 tGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACggg                                                     <  2:1776890/150‑1 (MQ=255)                                                                                                 tGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACggg                                                     >  2:577731/1‑150 (MQ=255)                                                                                                 tGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACggg                                                     <  1:1527958/150‑1 (MQ=255)                                                                                                      ttCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAAcc                                                >  1:1420989/1‑150 (MQ=255)                                                                                                         tgcctgcTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                             <  2:61197/150‑1 (MQ=255)                                                                                                         tgcctgcTGGAAGAAGACTTCAGGCCTGAGCCTGTAGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGGAACGGGAAACCGGt                                             >  1:238593/1‑150 (MQ=255)                                                                                                            ctgctgGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGt                                                                           >  1:359693/1‑117 (MQ=255)                                                                                                            ctgctgGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGt                                                                           <  2:359693/117‑1 (MQ=255)                                                                                                             ggcgtcaatgtcatacagcgtgcGAGGCTGTCGAACA‑GCAGGTACGACAGGCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                             <  1:746488‑M1/123‑1 (MQ=255)                                                                                                             ggcgtcaatgtcatacagcgtgcGAGGCTGTCGAACA‑GCAGGTACGACAGGCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                             >  2:746488‑M1/24‑146 (MQ=255)                                                                                                                tgGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTccc                                      >  1:1449970/1‑150 (MQ=255)                                                                                                                 ggAAGAAGACTTCAGGGCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCa                                     <  2:95046/150‑1 (MQ=255)                                                                                                                 ggAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCAAACCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCa                                     >  1:1060707/1‑150 (MQ=255)                                                                                                                  gaagaaGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAg                                    >  2:646765/1‑150 (MQ=255)                                                                                                                   aagaagACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAgg                                   <  2:1449970/150‑1 (MQ=255)                                                                                                                     gaagaCTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                             >  2:64146/1‑138 (MQ=255)                                                                                                                     gaagaCTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                             <  1:64146/138‑1 (MQ=255)                                                                                                                     gaagaCTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGcc                                 >  2:1101957/1‑150 (MQ=255)                                                                                                                          tacagcgtgcGAGGCTGTCGAACA‑GCAGGTACGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                             <  1:172632‑M1/123‑1 (MQ=255)                                                                                                                          tacagcgtgcGAGGCTGTCGAACA‑GCAGGTACGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGt                                             >  2:172632‑M1/11‑133 (MQ=255)                                                                                                                          cTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGTTTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTCca                            <  1:1040276/150‑1 (MQ=255)                                                                                                                               ggCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAgg                       <  2:44860/150‑1 (MQ=255)                                                                                                                                  cTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAGGGcc                    >  1:900062/1‑150 (MQ=255)                                                                                                                                       ccTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAGGGCCGGACa               >  1:862197/1‑150 (MQ=255)                                                                                                                                                ccTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAGGGCCGGACAATGCCGGtg      <  1:12431/150‑1 (MQ=255)                                                                                                                                                    ctgTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTCCGTCCCAGGCCCTTCAGCAGGGCCGGCCAATGCCGGtgctgc  >  2:365000/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CAGAAACTGGGGGCCGAGGCGCAGGATGAGTGGGTGCGCTACCGGCTGGATGGCGACAAGCTGGCGGCATTCGCCCGCGGGTGATCGCCCGGTTTCTGCACTTCTGCCTGCTGGAAGAAGACTTCAGGCCTGAGCCTGTCGTACCTGCTGTT‑CGACAGCCTCGCACGCTGTATGACATTGACGCCAAGCCCGCCCCATGCTCTCCTTATGCCTCGCGTTCAGGTGCCCCTCAGGGGGTGAAACGGGAAACCGGTGCGTCCCAGGCCCTTCAGCAGGGCCGGACAATGCCGGTGCTGC  >  NC_002947/3857387‑3857683

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑