Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,900,070	T→C	10.1%	intergenic (+31/‑467)	PP_4308 → / → PP_4309	AsnC family transcriptional regulator/NCS1 family transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,900,070	0	T	C	10.1%	177.4 / 10.0	69	intergenic (+31/‑467)	PP_4308/PP_4309	AsnC family transcriptional regulator/NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base T (35/27);  new base C (4/3);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.77e-01

TGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAAC  >  NC_002947/4899922‑4900216                                                                                                                                                      |                                                                                                                                                   tgtgCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGt                                                                                                                                                    <  1:338100/149‑1 (MQ=255)  gtgCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTg                                                                                                                                                   >  7:288005/1‑149 (MQ=255)              gAGCTGCTGGACACCCTGCTCGACCTCATGAATGTGCTGCCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGCGGGCACGCCTGCGAagag                                                                                                                                        <  7:176895/148‑1 (MQ=38)              gAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAagag                                                                                                                                        <  1:144756/148‑1 (MQ=255)                ggagcTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTTGCAGGCGTGCCCGCGAAGAGGc                                                                                                                                      <  5:364611/146‑1 (MQ=255)                 cagcTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCa                                                                                                                                    <  5:334807/149‑1 (MQ=255)                  agcTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAg                                                                                                                                   >  7:218050/1‑149 (MQ=255)                       gACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTAc                                                                                                                                <  5:287994/147‑1 (MQ=255)                          acCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGa                                                                                                                            >  4:263910/1‑148 (MQ=255)                            ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaa                                                                                                                           >  6:209205/1‑147 (MQ=255)                            ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaa                                                                                                                           >  1:368575/1‑147 (MQ=255)                            ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGCGGGCACGCCTGCGAAGAGGCCAGcaca                                                                                                                               >  4:64749/1‑143 (MQ=38)                            ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGCGGGCACGCCTGCGAAGAGGCCAGcaca                                                                                                                               <  3:64749/143‑1 (MQ=38)                            ccTGCTCGACCTCATGAATGTGCTGGCCGACACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaa                                                                                                                           >  1:32750/1‑147 (MQ=255)                             cTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCTGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaa                                                                                                                           >  2:122651/1‑146 (MQ=255)                                   gACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGa                                                                                                                    <  6:30228/147‑1 (MQ=255)                                   gACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAg                                                                                                                   >  3:113209/1‑148 (MQ=255)                                       tCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAgg                                                                                                                             >  3:122003/1‑134 (MQ=255)                                       tCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAgg                                                                                                                             <  4:122003/134‑1 (MQ=255)                                         aTGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAttttt                                                                                                            >  8:41557/1‑149 (MQ=255)                                         aTGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAttttt                                                                                                            >  5:335814/1‑149 (MQ=255)                                          tGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAtttttt                                                                                                           <  2:32750/149‑1 (MQ=255)                                                   tGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCAAAATGGATTGATTGCCTCATGTGATGCCCTCTACGCAGGCGTGCCCGCGAAGAGGCCAGTAAAGGAAAAGATGAGTATTTTTTTTAcaca                                                                                                    >  1:303797/1‑148 (MQ=255)                                                     gCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACaaaaa                                                                                                >  5:199226/1‑149 (MQ=255)                                                      ccGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTTATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAAc                                                                                               <  5:209205/149‑1 (MQ=255)                                                          gaCCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTaa                                                                                           <  4:113209/149‑1 (MQ=255)                                                          gaCCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTaa                                                                                           <  3:263910/149‑1 (MQ=255)                                                             cAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCGGGCACGCCTGCGAaga                                                                                                                                         >  8:33549/1‑100 (MQ=38)                                                             cAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCGGGCACGCCTGCGAaga                                                                                                                                         <  7:33549/100‑1 (MQ=38)                                                             cAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCa                                                                                        >  8:299807/1‑149 (MQ=255)                                                               aTACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATa                                                                                      <  6:335814/149‑1 (MQ=255)                                                                  ccGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGt                                                                                   <  2:368575/149‑1 (MQ=255)                                                                     cgATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGTTgg                                                                                 >  4:111098/1‑147 (MQ=255)                                                                      gATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACAATAGGTGCCa                                                                               >  5:130655/1‑149 (MQ=255)                                                                          atcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAttt                                                                            <  1:25608/148‑1 (MQ=255)                                                                           tcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAttttt                                                                          >  2:341986/1‑149 (MQ=255)                                                                           tcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAttttt                                                                          >  8:266255/1‑149 (MQ=255)                                                                           tcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAttttt                                                                          <  6:199226/149‑1 (MQ=255)                                                                           gcaAGAAGACCCCGGTCAAGAGGCGGTTGCCGC‑CAATTGATTGATTGGCTGCTGTGCTGTCCTGTTCGCAGGCGTGCCCTCGATGAGGCGAGTTCTGGAAATGCTGAGTATTTTTTGTACACTAATAGGTTACCATAGGTGGCAttttt                                                                          <  8:150233/148‑1 (MQ=255)                                                                             aagaagACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgt                                                                        <  8:288005/149‑1 (MQ=255)                                                                              agaagaCCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgtg                                                                       >  5:386588/1‑149 (MQ=255)                                                                                     cccGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACtt                                                                >  5:226187/1‑149 (MQ=255)                                                                                     cccGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACtt                                                                >  4:99646/1‑149 (MQ=255)                                                                                               cggcggTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTgg                                                      >  6:128818/1‑149 (MQ=255)                                                                                                  cggTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGgcgc                                                   <  6:226187/149‑1 (MQ=255)                                                                                                   ggTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGcc                                                  <  2:303797/149‑1 (MQ=255)                                                                                                    gTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGccc                                                 <  6:130655/149‑1 (MQ=255)                                                                                                              aaTGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGtt                                        <  7:266255/149‑1 (MQ=255)                                                                                                                 ggATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTac                                      >  4:78624/1‑148 (MQ=255)                                                                                                                  gattgattGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTaca                                     <  5:128818/148‑1 (MQ=255)                                                                                                                   attgattgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTacac                                    >  8:11302/1‑148 (MQ=255)                                                                                                                    ttgattgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCACCAATTCGTTacacac                                  >  6:239526/1‑149 (MQ=255)                                                                                                                        ttgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACttt                               >  7:340902/1‑148 (MQ=255)                                                                                                                        ttgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTc                              <  6:386588/149‑1 (MQ=255)                                                                                                                          ggCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCa                            >  1:128456/1‑149 (MQ=255)                                                                                                                                      ggCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTAtttt                 >  5:199735/1‑148 (MQ=255)                                                                                                                                       gCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTATTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTGAGTGGCGCCCCAATTCGTTACACACTTTCGAAGCTCTAtataga               >  6:190678/1‑144 (MQ=255)                                                                                                                                         ctctTCGCGGGCACGCCTGCGAAGAGGCCAGCACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAc              >  8:131958/1‑148 (MQ=39)                                                                                                                                         ctctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTacac                                                                                                     >  6:97414/1‑61 (MQ=255)                                                                                                                                         ctctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTacac                                                                                                     <  5:97414/61‑1 (MQ=255)                                                                                                                                         ctctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTacac                                                                                                     >  4:338055/1‑61 (MQ=255)                                                                                                                                         ctctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTacac                                                                                                     <  3:338055/61‑1 (MQ=255)                                                                                                                                         ctctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAc              >  5:94090/1‑148 (MQ=255)                                                                                                                                          tctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCGATTCGTTACACACTTTCCACCCTCTATTTTGAcc             >  2:170932/1‑148 (MQ=255)                                                                                                                                          tctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAcc             >  4:194008/1‑148 (MQ=255)                                                                                                                                           ctTCGCGGGCACGCCTGCGAAGAGGCCAGCACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTGTCCACCCTCTATTTTGACCa            >  2:269729/1‑148 (MQ=39)                                                                                                                                           ctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTaca                                     <  3:174407/123‑1 (MQ=255)                                                                                                                                           ctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTaca                                     >  4:174407/1‑123 (MQ=255)                                                                                                                                                cAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGAt      <  7:299807/149‑1 (MQ=255)                                                                                                                                                    cGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGACATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAAc  <  8:340902/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAAC  >  NC_002947/4899922‑4900216

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