Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,128,236	A→T	7.7%	S1203C (AGC→TGC)	PP_1887 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,128,236	0	A	T	7.7%	230.3 / 7.1	78	S1203C (AGC→TGC)	PP_1887	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/37);  new base T (4/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

GACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGATGAGCAGGG  >  NC_002947/2128095‑2128384                                                                                                                                              |                                                                                                                                                     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:385746/1‑149 (MQ=255)          ggTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGcc                                                                                                                                      <  5:257578/149‑1 (MQ=255)           gTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGc                                                                                                                                       >  5:52014/1‑147 (MQ=255)           gTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGc                                                                                                                                       <  6:52014/147‑1 (MQ=255)              aaaCCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGAt                                                                                                                                  >  8:123026/1‑149 (MQ=255)              aaaCCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGAt                                                                                                                                  >  2:354365/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                                                                                                      >  6:332538/1‑136 (MQ=255)                                                       aacaTCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATc                                                                                                      <  5:332538/136‑1 (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                                                                                           >  7:254311/1‑140 (MQ=255)                                                              tctcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTggag                                                                                           <  8:254311/140‑1 (MQ=255)                                                              tctcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAgttgtt                                                                                   >  1:344240/1‑148 (MQ=255)                                                                tcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTc                                                                                  <  6:103888/147‑1 (MQ=255)                                                                tcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTc                                                                                  <  8:44454/147‑1 (MQ=255)                                                                tcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTc                                                                                  >  6:81604/1‑147 (MQ=255)                                                                tcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTc                                                                                  <  7:223836/147‑1 (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                                                               >  2:10682/1‑135 (MQ=255)                                                                                               ggCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATc                                                               <  1:10682/135‑1 (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                 >  2:319369/1‑147 (MQ=255)                                                                                                                                 tCCGGCTCATTTGTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGAt                 >  6:316752/1‑147 (MQ=255)                                                                                                                                 tCCGGCTCATTTGTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCa               >  3:265627/1‑149 (MQ=255)                                                                                                                                 tCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCa               <  7:286430/149‑1 (MQ=255)                                                                                                                                 tCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCa               <  3:318101/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)                                                                                                                                        caAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGtgatga        <  8:300397/149‑1 (MQ=255)                                                                                                                                              aGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGATGAGCAggg  <  1:336836/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                     GACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGATGAGCAGGG  >  NC_002947/2128095‑2128384

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