Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,565,158	G→A	6.1%	S57S (TCC→TCT)	pitB ←	phosphate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,565,158	0	G	A	6.1%	188.2 / 2.4	65	S57S (TCC→TCT)	pitB	phosphate transporter
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base A (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

AGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCTCACCAACACCCAGGCATC  >  NC_002947/1565010‑1565304                                                                                                                                                     |                                                                                                                                                    aGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCg                                                                                                                                                     <  3:61807/149‑1 (MQ=255) aGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCg                                                                                                                                                     >  4:61807/1‑149 (MQ=255)          ccGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAac                                                                                                                                            >  7:85056/1‑149 (MQ=255)               gccagcAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCg                                                                                                                                                     >  5:251158/1‑135 (MQ=255)               gccagcAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCg                                                                                                                                                     <  6:251158/135‑1 (MQ=255)                 cagcAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑Ga                                                                                                                                    <  5:351079/149‑1 (MQ=255)                 cagcAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑Ga                                                                                                                                    <  5:337856/149‑1 (MQ=255)                    cAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGACCGCCAAGGAAGTTGACCACGCCGGAGAAGACCACGGC‑GAGGt                                                                                                                                 >  8:105364/1‑149 (MQ=255)                      aCGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTgc                                                                                                                               >  7:139888/1‑149 (MQ=255)                           aagaCCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGACGGAGAAGAACACGCC‑GAGGTGcggcgg                                                                                                                          >  7:76070/1‑149 (MQ=255)                              aCCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCAt                                                                                                                       <  3:296301/149‑1 (MQ=255)                                    gCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCttt                                                                                                                 <  8:139888/149‑1 (MQ=255)                                      cAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTgg                                                                                                               <  8:193682/149‑1 (MQ=255)                                       aGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGAGGCATGGCTTTgg                                                                                                               >  6:310287/1‑148 (MQ=255)                                            gTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGa                                                                                                         <  5:310287/149‑1 (MQ=255)                                                 cGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACg                                                                                                    >  5:70322/1‑149 (MQ=255)                                                   gtgtTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCAGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACgg                                                                                                   >  4:344908/1‑148 (MQ=255)                                                   gtgtTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCAGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACgg                                                                                                   >  4:339445/1‑148 (MQ=255)                                                      ttCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACAGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGc                                                                                               >  7:81028/1‑149 (MQ=255)                                                          caTTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTAcc                                                                                           >  7:32408/1‑149 (MQ=255)                                                           aTTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCGCTACACCGCCGAGCAGCACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACcg                                                                                          >  1:77689/1‑149 (MQ=255)                                                           aTTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACcg                                                                                          >  2:357080/1‑149 (MQ=255)                                                               aTCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGtt                                                                                      <  2:329399/149‑1 (MQ=255)                                                               aTCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGCCTTTGGTATACATGACGGTAGCTACCGCGtg                                                                                      >  4:271962/1‑148 (MQ=255)                                                                tccgcagcTCTCCCCGCAACCGGTGCCCGATGGCATAGGCCACCCCAACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAACACCCCGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTg                                                                                     <  8:121490/146‑1 (MQ=255)                                                                  agcagcTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTTCACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTgg                                                                                    <  3:114229/148‑1 (MQ=255)                                                                  agcagcTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGc                                                                                   <  2:245391/149‑1 (MQ=255)                                                                   gcagcTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCg                                                                                  >  6:171333/1‑149 (MQ=255)                                                                     agcTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGgt                                                                                <  8:81028/149‑1 (MQ=255)                                                                        tctACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTa                                                                                                 <  4:205504/129‑1 (MQ=255)                                                                        tctACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTa                                                                                                 >  3:205504/1‑129 (MQ=255)                                                                          tACCGGAAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCAt                                                                           <  5:56887/149‑1 (MQ=255)                                                                            ccGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCAt                                                                           <  6:15716/147‑1 (MQ=255)                                                                                 aaCAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGcc                                                                    <  7:105364/149‑1 (MQ=255)                                                                                 aaCAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGc                                                                     >  2:287217/1‑148 (MQ=255)                                                                                      gTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATtga                                                               >  6:164086/1‑149 (MQ=255)                                                                                      gTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATtga                                                               <  2:98/149‑1 (MQ=255)                                                                                       tGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATtgat                                                              <  3:271962/149‑1 (MQ=255)                                                                                         cACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATtgatg                                                             <  6:70322/148‑1 (MQ=255)                                                                                          aCGATGGCATAGGCCCCCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGaa                                                           >  1:235400/1‑149 (MQ=255)                                                                                                    aGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGCACGGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGc                                                  >  1:107496/1‑149 (MQ=255)                                                                                                          ccccTACACCGCCGAGCAGAACGCCAAGTAAGTTTAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACg                                           <  1:354701/149‑1 (MQ=255)                                                                                                           cccTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACg                                           >  4:334201/1‑148 (MQ=255)                                                                                                             cTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCAACCGCGTTCGCGCTGTCATGAAAGCCATtgat                                                              >  5:367325/1‑127 (MQ=255)                                                                                                               acacCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGc                                                  <  2:320157/137‑1 (MQ=255)                                                                                                               acacCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGc                                                  >  1:320157/1‑137 (MQ=255)                                                                                                                cacCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGTGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGt                                      >  7:57599/1‑148 (MQ=255)                                                                                                                           aGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGagcag                          >  1:333459/1‑149 (MQ=255)                                                                                                                               ggCCAAGGAAGTTGAACTCGCCGGAGATTTACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAgg                      <  2:77689/148‑1 (MQ=255)                                                                                                                               cGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAgg                      <  3:334201/149‑1 (MQ=255)                                                                                                                                  cAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCt                    >  4:208980/1‑148 (MQ=255)                                                                                                                                  cAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCt                    >  1:46318/1‑148 (MQ=255)                                                                                                                                  cAAGGAAGTTGAACACGCCAGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCTc                   <  3:344908/149‑1 (MQ=255)                                                                                                                                  cAAGGAAGTTGAACACGCCAGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCTc                   <  3:339445/149‑1 (MQ=255)                                                                                                                                   aaGGAAGTTGAACACGCCGGAGAAGAACACGTC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGa                              <  3:63732/137‑1 (MQ=255)                                                                                                                                   aaGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGACATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGa                              >  4:63732/1‑137 (MQ=255)                                                                                                                                    aGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCTCAc                 >  3:54976/1‑149 (MQ=255)                                                                                                                                     ggAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCTCAcc                <  1:287217/149‑1 (MQ=255)                                                                                                                                         gTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCTCACCAaca            <  8:32408/149‑1 (MQ=255)                                                                                                                                            gAACACGCCGGAGAAGAACACGGC‑GAGGTGGGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGATGGTCAAGGCGAGCAGCAGGCTCACCaccccc          >  4:199452/1‑143 (MQ=255)                                                                                                                                             aaCACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCTCACCAACAccc         <  5:164086/148‑1 (MQ=255)                                                                                                                                              acacGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCATCAGGCTCACCCaccacaag       >  6:5440/1‑143 (MQ=255)                                                                                                                                              acacGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCTCACCAACACCCa        <  1:357080/148‑1 (MQ=255)                                                                                                                                              acacGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCTCACCAACACCCAg       >  6:10569/1‑149 (MQ=255)                                                                                                                                              acacGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCTCACCAACACCCAg       <  5:171333/149‑1 (MQ=255)                                                                                                                                                 cGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTTGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCTCACCAACACCCAGGCa    <  2:46318/149‑1 (MQ=255)                                                                                                                                                    cGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGt                                              >  3:187351/1‑104 (MQ=255)                                                                                                                                                    cGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGt                                              <  4:187351/104‑1 (MQ=255)                                                                                                                                                    cGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCTCACCAACACCCAGGCATc  >  4:328985/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTACCGGCAACAGGTGCACGATGGCATAGGCCACCCCTACACCGCCGAGCAGAACGCCAAGGAAGTTGAACACGCCGGAGAAGAACACGGC‑GAGGTGCGGCGGCATGGCTTTGGTATAGATGACGGTAGCTACCGCGTTGGCGGTGTCATGAAAGCCATTGATGAACTCGAAGGCGAGTACGAAGGTCAAGGCGAGCAGCAGGCTCACCAACACCCAGGCATC  >  NC_002947/1565010‑1565304

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