Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,128,229	A→T	8.7%	A1200A (GCA→GCT)	PP_1887 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,128,229	0	A	T	8.7%	223.3 / 7.8	78	A1200A (GCA→GCT)	PP_1887	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (36/35);  new base T (4/3);  total (40/38)
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.98e-01

GCGCCCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAG  >  NC_002947/2128086‑2128372                                                                                                                                                |                                                                                                                                                gcgcCCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAt                                                                                                                                            >  8:246479/1‑149 (MQ=255)  cgcCCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGctcatttg                                                                                                                                           >  5:40445/1‑145 (MQ=255)    cccGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATg                                                                                                                                           >  5:30529/1‑147 (MQ=255)    cccGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATg                                                                                                                                           <  6:30529/147‑1 (MQ=255)          gACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCTCATTTGTGCCGGa                                                                                                                                    <  6:40445/148‑1 (MQ=255)          gACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGAt                                                                                                                                   <  7:246479/149‑1 (MQ=255)             cccTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCTCATTTGTGCCGGATcgc                                                                                                                                <  1:75949/149‑3 (MQ=255)             cccTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCAt                                                                                                                                >  7:32857/1‑149 (MQ=255)              ccTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCAtt                                                                                                                               >  8:125922/1‑149 (MQ=255)               cTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTc                                                                                                                              <  8:255421/149‑1 (MQ=255)                tCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTcc                                                                                                                             >  8:223836/1‑149 (MQ=255)                  cGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTcc                                                                                                                             <  8:268133/147‑1 (MQ=255)                   ggTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGcc                                                                                                                          <  5:257578/149‑1 (MQ=255)                   ggTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGcc                                                                                                                          >  5:385746/1‑149 (MQ=255)                    gTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGc                                                                                                                           >  5:52014/1‑147 (MQ=255)                    gTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGc                                                                                                                           <  6:52014/147‑1 (MQ=255)                       aaaCCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGAt                                                                                                                      >  2:354365/1‑149 (MQ=255)                       aaaCCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGAt                                                                                                                      >  8:123026/1‑149 (MQ=255)                       aaaCCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGAt                                                                                                                      >  2:300675/1‑149 (MQ=255)                           cATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCTATCCGGCACAAATGAGCCGGATCATTCCGCCTGATgcgc                                                                                                                  <  2:93568/149‑1 (MQ=255)                           cATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATgcgc                                                                                                                  >  3:237394/1‑149 (MQ=255)                            aTTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACTTAAAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCTGATCATTCCGCCTGATgcgc                                                                                                                  >  8:195429/1‑148 (MQ=255)                            aTTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGa                                                                                                                                    >  8:343575/1‑130 (MQ=255)                            aTTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGa                                                                                                                                    <  7:343575/130‑1 (MQ=255)                                ccAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTg                                                                                                             >  4:211708/1‑149 (MQ=255)                                 cAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTgg                                                                                                            <  5:75565/149‑1 (MQ=255)                                  aGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGga                                                                                                           <  1:354365/149‑1 (MQ=255)                                  aGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGCATCATTCAGCCTGATGCGCTACTGga                                                                                                           >  6:33327/1‑149 (MQ=255)                                       gATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCt                                                                                                      <  3:211708/149‑1 (MQ=255)                                         tatTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTAt                                                                                                    <  5:47432/149‑1 (MQ=255)                                          atTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTAt                                                                                                    >  6:116300/1‑148 (MQ=255)                                             cGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGAc                                                                                                 <  4:237394/148‑1 (MQ=255)                                                cTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTc                                                                                              >  4:281944/1‑148 (MQ=255)                                                 ttCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGa                                                                                            >  5:103888/1‑149 (MQ=255)                                                              acaacaTCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTggag                                                                               <  5:116300/149‑1 (MQ=255)                                                                aacaTCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGga                                                                                                           >  4:15273/1‑119 (MQ=255)                                                                aacaTCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGga                                                                                                           <  3:15273/119‑1 (MQ=255)                                                                aacaTCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATc                                                                                          <  5:332538/136‑1 (MQ=255)                                                                aacaTCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATc                                                                                          >  6:332538/1‑136 (MQ=255)                                                                aacaTCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTggagga                                                                             <  2:242701/149‑1 (MQ=255)                                                                   aTCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAgtt                                                                          <  7:123026/149‑1 (MQ=255)                                                                      ctctCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTTCTGGAGACCTATGACGTCGATCACGACCTGGAGGAgttgtt                                                                       >  5:321396/1‑149 (MQ=255)                                                                      ctctCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAgttgt                                                                        <  1:300675/148‑1 (MQ=255)                                                                       tctcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTggag                                                                               <  8:254311/140‑1 (MQ=255)                                                                       tctcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTggag                                                                               >  7:254311/1‑140 (MQ=255)                                                                       tctcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAgttgtt                                                                       >  1:344240/1‑148 (MQ=255)                                                                         tcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTc                                                                      <  6:103888/147‑1 (MQ=255)                                                                         tcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTc                                                                      <  8:44454/147‑1 (MQ=255)                                                                         tcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTc                                                                      <  7:223836/147‑1 (MQ=255)                                                                         tcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTc                                                                      >  6:81604/1‑147 (MQ=255)                                                                           ggTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAg                                                                   <  5:81604/148‑1 (MQ=255)                                                                             tGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATc                                                                <  5:328115/149‑1 (MQ=255)                                                                                  aCCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGTGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCa                                                            >  6:352488/1‑148 (MQ=255)                                                                                  aCCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCa                                                            >  2:336836/1‑148 (MQ=255)                                                                                       cAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATa                                                        <  5:33327/147‑1 (MQ=255)                                                                                            cagACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCg                                                 >  1:54665/3‑149 (MQ=255)                                                                                              gACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCa                                               <  8:32857/149‑1 (MQ=255)                                                                                                   gTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTgg                                          <  6:385746/149‑1 (MQ=255)                                                                                                        ggCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATc                                                   <  1:10682/135‑1 (MQ=255)                                                                                                        ggCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATc                                                   >  2:10682/1‑135 (MQ=255)                                                                                                           aGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGAt                                  <  7:125922/149‑1 (MQ=255)                                                                                                             cGGACCTTGTACGAATTCGAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGAttt                                >  8:286430/1‑149 (MQ=255)                                                                                                              ggACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTg                               >  8:273397/1‑149 (MQ=255)                                                                                                               gACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGc                              <  2:344240/149‑1 (MQ=255)                                                                                                                aCCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCa                             <  7:195429/149‑1 (MQ=255)                                                                                                                      tACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTacga                       >  4:318101/1‑149 (MQ=255)                                                                                                                        cGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTacgacg                     >  7:300397/1‑149 (MQ=255)                                                                                                                               aaGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTAAGACGGCAATgg              >  8:70921/1‑149 (MQ=255)                                                                                                                                 gAATAGCGATCCTGCTCTATTGAGCCGGATCATTGCGCCTGATTCGCTACTGTAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTAt                                      <  3:34792/123‑1 (MQ=255)                                                                                                                                     aGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCt        >  8:20424/1‑149 (MQ=255)                                                                                                                                        gATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGAt     >  1:334331/1‑149 (MQ=255)                                                                                                                                         aTCCGGCTCATTTGTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGAt     >  4:318473/1‑148 (MQ=255)                                                                                                                                         aTCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGAt     <  5:352488/148‑1 (MQ=255)                                                                                                                                          tCCGGCTCATTTGTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGAt     >  6:316752/1‑147 (MQ=255)                                                                                                                                          tCCGGCTCATTTGTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGAt     >  2:319369/1‑147 (MQ=255)                                                                                                                                          tCCGGCTCATTTGTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCa   >  3:265627/1‑149 (MQ=255)                                                                                                                                          tCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCa   <  3:318101/149‑1 (MQ=255)                                                                                                                                          tCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCa   <  7:286430/149‑1 (MQ=255)                                                                                                                                           ccGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAg  >  7:17929/1‑149 (MQ=255)                                                                                                                                                acaAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATgg              >  7:90215/1‑132 (MQ=255)                                                                                                                                                acaAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATgg              <  8:90215/132‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                GCGCCCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAG  >  NC_002947/2128086‑2128372

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