Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,900,051	G→A	7.5%	intergenic (+12/‑486)	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,051	0	G	A	7.5%	191.0 / 3.1	66	intergenic (+12/‑486)	PP_4308/PP_4309	AsnC family transcriptional regulator/NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base A (2/3);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

TTCGTCGCCCGACTGCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCA  >  NC_002947/4899905‑4900190                                                                                                                                                    |                                                                                                                                            ttCGTCGCCCGACTGCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCt                                                                                                                                           >  6:334807/1‑149 (MQ=255)    gtcgCCCGACTGCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATT‑ATTGGCTCCTGTGCTGGc                                                                                                                                        <  1:105390/148‑1 (MQ=255)        cccGACTGCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCctct                                                                                                                                    >  6:287994/1‑149 (MQ=255)        cccGACTGCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCctct                                                                                                                                    >  2:25608/1‑149 (MQ=255)               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)                  tgtgCGCTCGATGGAGCAGCTGGACACCCTTCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGcggtcac                                                                                                                          <  1:367598/149‑7 (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                                                                                                 >  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‑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                                                                                                   <  4:122003/134‑1 (MQ=255)                                                        tCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAgg                                                                                                   >  3:122003/1‑134 (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                                                                                                               <  7:33549/100‑1 (MQ=38)                                                                              cAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCGGGCACGCCTGCGAaga                                                                                                               >  8:33549/1‑100 (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                                                <  6:199226/149‑1 (MQ=255)                                                                                            tcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAttttt                                                >  8:266255/1‑149 (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)                                                                                                                                                    |                                                                                                                                            TTCGTCGCCCGACTGCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGC‑CAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCA  >  NC_002947/4899905‑4900190

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