Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,390,101	T→A	7.7%	intergenic (+190/‑532)	PP_4738 → / → PP_mr56	hypothetical protein/C4 ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,390,101	0	T	A	7.7%	172.0 / 4.4	65	intergenic (+190/‑532)	PP_4738/PP_mr56	hypothetical protein/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base T (29/31);  new base A (2/3);  total (31/34)
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.62e-01

ATGGCGTATACTCCGTCTCACGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGTAAACCCGAACTGATTGTCATACAGAA  >  NC_002947/5389957‑5390239                                                                                                                                                 |                                                                                                                                           aTGGCGTATACTCCGTCTCACGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGtt                                                                                                                                        <  7:98618/149‑1 (MQ=255) aTGGCGTATACTCCGTCTCACGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGtt                                                                                                                                        <  3:342898/149‑1 (MQ=255)     cGTATACTCCGTCTCACGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTATAAAAGTTGcaa                                                                                                                                    <  3:129669/149‑3 (MQ=255)       tataCTCCGTCTCACGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAAt                                                                                                                                  <  8:53456/149‑1 (MQ=255)       tataCTCCGTCTCACGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAAt                                                                                                                                  <  3:173850/149‑1 (MQ=255)        ataCTCCGTCTCACGTACAGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCAATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAAt                                                                                                                                  <  8:77322/148‑1 (MQ=255)         taCTCCGTCTCACGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGa                                                                                                                                <  4:162586/149‑1 (MQ=255)                tctcACGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGAAACTTTTATAGTTGCCAATGACGAAAAc                                                                                                                         >  7:95264/1‑149 (MQ=255)                  tcACGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTa                                                                                                                       >  8:241668/1‑149 (MQ=255)                    aCGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACg                                                                                                                     >  7:109308/1‑149 (MQ=255)                    aCGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACg                                                                                                                     <  7:241668/149‑1 (MQ=255)                    aCGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACg                                                                                                                     >  5:67470/1‑149 (MQ=255)                    aCGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTAAAGTTGCCAATGACGAAAACTACg                                                                                                                     >  6:356670/1‑149 (MQ=255)                         ccGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCGGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGcgcg                                                                                                                <  3:262277/149‑1 (MQ=255)                           gTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTATAAAAGTTGc                                                                                                                                      <  7:329442/125‑1 (MQ=255)                           gTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTATAAAAGTTGc                                                                                                                                      >  8:329442/1‑125 (MQ=255)                             aaGGTGCGGTAGCGTAAAGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTATGTGCATGCCGAGTTGGTATCAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGAATAAAACTACGGAGCTGAg                                                                                                             >  6:120370/1‑148 (MQ=255)                                           tAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCAACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCg                                                                                              >  8:101320/1‑149 (MQ=255)                                               gCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGa                                                                                                              >  3:37147/1‑129 (MQ=255)                                               gCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGa                                                                                                              <  4:37147/129‑1 (MQ=255)                                                 tAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGa                                                                                                              <  2:32481/127‑1 (MQ=255)                                                 tAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAg                                                                                          >  6:332259/1‑147 (MQ=255)                                                            cGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGAc                                                                                                                               <  8:196743/99‑1 (MQ=255)                                                            cGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGAc                                                                                                                               >  7:196743/1‑99 (MQ=255)                                                               ttAGGATTCTACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCg                                                                          <  2:196671/149‑1 (MQ=255)                                                                   gATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGt                                                                                                            >  8:18529/1‑111 (MQ=255)                                                                   gATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGt                                                                                                            <  7:18529/111‑1 (MQ=255)                                                                    aTTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACATAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGc                                                                                               <  3:281709/123‑1 (MQ=255)                                                                    aTTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCt                                                                     >  4:182734/1‑149 (MQ=255)                                                                    aTTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGACAATGACGAAAACTACGGCTCGGAGTACAAgaaggctga                                                                                               >  4:281709/1‑122 (MQ=255)                                                                       cGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCt                                                                  >  5:246986/1‑149 (MQ=255)                                                                        gACGCCGGTAGCGAAACTCTAGGTACATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTg                                                                 >  1:38818/1‑149 (MQ=255)                                                                         aCGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTgg                                                                <  5:332259/149‑1 (MQ=255)                                                                           gCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTa                                                              <  8:109308/149‑1 (MQ=255)                                                                                tAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTAAAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTc                                                         <  5:356670/149‑1 (MQ=255)                                                                                    gAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGCCGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCt                                                     >  5:257773/1‑149 (MQ=255)                                                                                    gAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTc                                                                   <  5:313908/135‑1 (MQ=255)                                                                                    gAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTc                                                                   >  6:313908/1‑135 (MQ=255)                                                                                    gAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGc                                                      >  6:214405/1‑148 (MQ=255)                                                                                    gAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCt                                                     >  5:226257/1‑149 (MQ=255)                                                                                    gAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTATGGTAGCTTCGGCt                                                     >  1:236002/1‑149 (MQ=255)                                                                                     aaaCTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCATCTTGACCCGTTGCTTCTGGTAGCTTCGGCtg                                                    >  8:238894/1‑148 (MQ=255)                                                                                                  aTGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTgg                                                                <  8:62219/124‑1 (MQ=255)                                                                                                  aTGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTgg                                                                >  7:62219/1‑124 (MQ=255)                                                                                                  aTGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAACAATCACTAg                                       >  1:338510/1‑149 (MQ=255)                                                                                                         tttGGTAACAGAACTCGTAAATCCACTGTAGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAATCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGc                                <  5:120370/148‑1 (MQ=255)                                                                                                         gTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGc                                <  6:246986/149‑1 (MQ=255)                                                                                                          ttGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGc                                <  7:101320/148‑1 (MQ=255)                                                                                                            ggTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTTGCTGCGTAAGCAGTCCTCTTGACCCGTTCCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTg                             >  3:200571/1‑149 (MQ=255)                                                                                                             gTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGt                            <  2:38818/149‑1 (MQ=255)                                                                                                             gTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTCCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGt                            <  6:226257/149‑1 (MQ=255)                                                                                                              tAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGt                            >  6:99211/1‑148 (MQ=255)                                                                                                              tAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGt                            >  6:280128/1‑148 (MQ=255)                                                                                                                      cTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCGGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGTAAACCCGaa                   <  8:95264/149‑1 (MQ=255)                                                                                                                         gTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGt                            <  8:97875/137‑1 (MQ=255)                                                                                                                         gTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGt                            >  7:97875/1‑137 (MQ=255)                                                                                                                           aaaTCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGTAAACCCGAACTGAt              >  1:190310/1‑149 (MQ=255)                                                                                                                                 aCTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGTAAACCCGAACTGAtt             <  3:302212/144‑1 (MQ=255)                                                                                                                                 aCTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGTAAACCCGAACTGAtt             >  4:302212/1‑144 (MQ=255)                                                                                                                                  cTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGTAAACCCGAACTGATTGTCATa       <  6:257773/149‑1 (MQ=255)                                                                                                                                     ttGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGTAAACCCGAACTGATTGTCATAc      <  5:280128/147‑1 (MQ=255)                                                                                                                                     ttGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGTAAACCCGAACTGATTGTCATAc      <  5:99211/147‑1 (MQ=255)                                                                                                                                        cAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGTAAACCCGaa                   <  8:105316/131‑1 (MQ=255)                                                                                                                                        cAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGTAAACCCGaa                   >  7:105316/1‑131 (MQ=255)                                                                                                                                        cAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGTAAACCCGAACTGATTGTCATACAGaa  <  3:182734/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                           ATGGCGTATACTCCGTCTCACGTACCGTAAGGTGCGGTAGCGTAACGCTAGCGTGGGGCCGATTAGGATTCGACGCCGGTAGCGAAACTCTAGGTGCATGCCGAGTTGGTAACAGAACTCGTAAATCCACTGTTGCAACTTTTATAGTTGCCAATGACGAAAACTACGGCGCGGAGTACAAGCTGGCTGCGTAAGCAGTCCTCTTGACCCGTTGCTTCTGGTAGCTTCGGCTCCAGCAATCACTAGGGGATGCCTGTAAACCCGAACTGATTGTCATACAGAA  >  NC_002947/5389957‑5390239

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