Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	990,352	T→A	7.6%	V78D (GTC→GAC)	PP_0856 →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	990,352	0	T	A	7.6%	194.6 / 5.4	66	V78D (GTC→GAC)	PP_0856	lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base A (2/3);  total (33/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.29e-01

GGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGA  >  NC_002947/990206‑990499                                                                                                                                                   |                                                                                                                                                    ggCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCt                                                                                                                                                   >  1:235583/1‑149 (MQ=255)   ccGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTc                                                                                                                                                 <  6:225448/149‑1 (MQ=255)           aCCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGACTTCGccgtt                                                                                                                                           <  7:293624/147‑4 (MQ=255)            ccAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGAc                                                                                                                                        <  2:235583/149‑1 (MQ=255)             ccaGTTCACCGAGGCAGTGGGGCTGAAGACGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGccc                                                                                                                                             <  6:332519/141‑1 (MQ=255)             cAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGccc                                                                                                                                             >  5:332519/1‑143 (MQ=255)             cAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGAcc                                                                                                                                       >  3:148685/1‑149 (MQ=255)               aGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCg                                                                                                                                      >  2:243104/1‑148 (MQ=255)                   cACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGACTTCGccgt                                                                                                                                            <  3:234525/138‑3 (MQ=255)                   cACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGACTTCGccgt                                                                                                                                            >  4:234525/1‑136 (MQ=255)                      cGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCg                                                                                                                              >  6:338364/1‑149 (MQ=255)                                gtgCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACgtggt                                                                                                                    <  7:258259/149‑1 (MQ=255)                                   cTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGtgtg                                                                                                                 >  6:285375/1‑149 (MQ=255)                                      aagaagCAGTGGAGCCGTTCGATCGGTGACGTCCATGGTGAAAGCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATc                                                                                                                                  <  6:24792/129‑1 (MQ=255)                                      aagaagCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCTGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTATTCGCCCTCGACCGCATc                                                                                                                                  >  5:24792/1‑129 (MQ=255)                                       agaagCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGgaag                                                                                                             <  2:293354/149‑1 (MQ=255)                                          agCAGTGGTGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGaagaag                                                                                                          <  5:290722/149‑1 (MQ=255)                                           gCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGaagaag                                                                                                          >  5:314021/1‑148 (MQ=255)                                            cAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGa                                                                                                        >  6:154303/1‑149 (MQ=255)                                             aGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACgtggt                                                                                                                    >  2:149039/1‑136 (MQ=255)                                             aGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACgtggt                                                                                                                    <  1:149039/136‑1 (MQ=255)                                             aGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGa                                                                                                        >  3:115045/1‑148 (MQ=255)                                                 gAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTg                                                                                                    <  6:314021/148‑1 (MQ=255)                                                  aGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTg                                                                                                    >  7:25272/1‑147 (MQ=255)                                                       ttCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGc                                                                                             <  4:148685/149‑1 (MQ=255)                                                         cGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCa                                                                                            <  5:285375/148‑1 (MQ=255)                                                               gTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCa                                                                                            <  3:66471/142‑1 (MQ=255)                                                               gTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCa                                                                                            >  4:66471/1‑142 (MQ=255)                                                                tGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCTACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCg                                                                                    <  4:115045/149‑1 (MQ=255)                                                                tGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCg                                                                                    <  5:154303/149‑1 (MQ=255)                                                                    ggCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTg                                                                                 >  4:66428/1‑148 (MQ=255)                                                                         gggTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTg                                                                           <  8:25272/149‑1 (MQ=255)                                                                              aaaCCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCAGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCg                                                                        >  6:123075/1‑147 (MQ=255)                                                                                aCCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAgg                                                                    >  1:159325/1‑149 (MQ=255)                                                                                  cTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGc                                                                   >  3:1482/1‑148 (MQ=255)                                                                                  cTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGAACTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGc                                                                   >  6:30153/1‑148 (MQ=255)                                                                                      aacaCCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGtt                                                                            >  4:48033/1‑135 (MQ=255)                                                                                      aacaCCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGtt                                                                            <  3:48033/135‑1 (MQ=255)                                                                                      aacaCCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATgg                                                              >  4:25946/1‑149 (MQ=255)                                                                                               gTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCAt                                                     <  2:159325/149‑1 (MQ=255)                                                                                                tACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTcc                                                                                     <  4:317155/116‑1 (MQ=255)                                                                                                tACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTcc                                                                                     >  3:317155/1‑116 (MQ=255)                                                                                                tACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATg                                                    >  4:67233/1‑149 (MQ=255)                                                                                                 aCCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATg                                                    >  8:102317/1‑148 (MQ=255)                                                                                                    ggCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCg                                                <  3:67233/149‑1 (MQ=255)                                                                                                         tCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTAcc                                           <  3:25946/149‑1 (MQ=255)                                                                                                              aaCGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTc                                        >  1:165873/1‑147 (MQ=255)                                                                                                                cGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAgg                                    >  1:253496/1‑149 (MQ=255)                                                                                                                cGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAgg                                    <  7:102317/149‑1 (MQ=255)                                                                                                                cGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAgg                                    <  2:165873/149‑1 (MQ=255)                                                                                                                cGACCGTATCTACGCTTCAGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAgg                                    <  5:123075/149‑1 (MQ=255)                                                                                                                     gTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTAGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGaa                               <  1:243104/149‑1 (MQ=255)                                                                                                                          tACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCAt                          >  1:41865/1‑149 (MQ=255)                                                                                                                          tACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGGGGATGCAGGACCTGGAGCTGCAGGTTTCCGGGGCAGTTGGCCTAGGCTATGGCCTGGTCATGATCGGTACCCTCATGGGTGAAGTCAt                          >  6:263413/1‑149 (MQ=255)                                                                                                                              cTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCCCCGGTGACGTGGTGTGGTAGAACGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGcc                      <  6:349912/149‑1 (MQ=255)                                                                                                                                   gACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGa                 <  2:41865/149‑1 (MQ=255)                                                                                                                                     ggtcgagggCGAAGACTTCGCCGTTGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACt               >  5:88219/8‑149 (MQ=255)                                                                                                                                      gTCATCGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTc              <  5:338364/149‑1 (MQ=255)                                                                                                                                       tCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTcc             <  2:253496/149‑1 (MQ=255)                                                                                                                                         aaCGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTCCAg           >  1:112520/1‑149 (MQ=255)                                                                                                                                           cGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGCAGGCCCAGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCATGGACTCCAGCa         >  6:88596/1‑149 (MQ=255)                                                                                                                                               gAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCgg     >  4:322606/1‑149 (MQ=255)                                                                                                                                                aaGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGt    <  3:322606/149‑1 (MQ=255)                                                                                                                                                aaGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGt    >  4:329521/1‑149 (MQ=255)                                                                                                                                                aaGACTGCGCCGTTGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGt    <  6:88219/149‑1 (MQ=255)                                                                                                                                                  gTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTga  <  3:329521/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GGCCGAGCTGACCAAGTTCACCGAGGAAGTGGTGCTGAAGAAGCAGTGGAGCCGTTCGATCGGTGACGGCCAGGGTGAAACCTACAACACCTTGGTACCGGCCATCGAGAACGACCGTATCTACGCTTCCGACGTCAACGGCGAAGTCTTCGCCCTCGACCGCATCACCGGTGACGTGGTGTGGAAGAAGGACCTGGAGCTGCAGGTTTCCGGTGCAGTTGGCGTAGGCTATGGCCTGGTCATGCTCGGTACCCTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGA  >  NC_002947/990206‑990499

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