Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,900,065	A→G	10.7%	intergenic (+26/‑472)	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,065	0	A	G	10.7%	195.3 / 12.7	75	intergenic (+26/‑472)	PP_4308/PP_4309	AsnC family transcriptional regulator/NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base A (37/30);  new base G (5/3);  total (42/33)
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.72e-01

GCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGAT  >  NC_002947/4899919‑4900212                                                                                                                                                    |                                                                                                                                                    gCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACAAATACCGCGATCATCAAGAAGACCACCGTCAAGCGGCGGTTGACGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTACGCAgg                                                                                                                                                   >  7:150233/1‑149 (MQ=255)  cATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGc                                                                                                                                                  <  3:263037/149‑1 (MQ=255)  cATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGc                                                                                                                                                  <  6:119436/149‑1 (MQ=255)   aTGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGCGGGca                                                                                                                                                 >  5:298368/1‑148 (MQ=255)    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                                                                                                                       >  1:368575/1‑147 (MQ=255)                               ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaa                                                                                                                       >  6:209205/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‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGc                                                                                                                                                  >  1:375416/1‑113 (MQ=255)                                      gACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGc                                                                                                                                                  <  2:375416/113‑1 (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                                                                                                        >  5:335814/1‑149 (MQ=255)                                            aTGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAttttt                                                                                                        >  8:41557/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                                                                      >  8:266255/1‑149 (MQ=255)                                                                              tcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAttttt                                                                      >  2:341986/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                                                                                                 <  3:338055/61‑1 (MQ=255)                                                                                                                                            ctctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTacac                                                                                                 >  4:338055/1‑61 (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)                                                                                                                                                    |                                                                                                                                                    GCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGAT  >  NC_002947/4899919‑4900212

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