Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,342,103	T→G	8.6%	intergenic (+18/‑46)	yraO → / → PP_2058	transporter YraO/porin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,342,103	0	T	G	8.6%	186.4 / 3.4	64	intergenic (+18/‑46)	yraO/PP_2058	transporter YraO/porin
Reads supporting (aligned to +/- strand):  ref base T (27/31);  new base G (2/4);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCTGGTGGGCATGGTCGACCGCGATATCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTGCCGCCGGTTTCATCGAAGACAGC  >  NC_002947/2341957‑2342238                                                                                                                                                    |                                                                                                                                        ctgctgGTGGGCATGGTCGACCGCGATATCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCg                                                                                                                                       <  1:297001/149‑1 (MQ=255)    ctgGTGGGCATGGTCGACCGCGATATCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGgcg                                                                                                                                    >  7:198013/1‑149 (MQ=255)           gCATGGTCGACCGCGATATCGGTGATTTGCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCt                                                                                                                             >  1:207256/1‑149 (MQ=255)           gCATGGTCGACCGCGATATCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCt                                                                                                                             >  2:180739/1‑149 (MQ=255)              ttgTCGACCGCTATATCGGTGATTTCCAGAAAGCCACCGTTAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGAc                                                                                                                          <  8:337302/147‑1 (MQ=255)               ggTCGACCGCGATATCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACt                                                                                                                         >  2:274562/1‑149 (MQ=255)                gTCCACCGCGACATCGGGGATTTCCAGAACCCCACCGGGAAATCGGCCGTGCTCACCTCTCTCGTGATCACTGCGCTGCCACTGCAGACCGGCGCCTTG‑TCCTTCATGGGTTGAAGCACCCAGGGCGCCGCGCGGCGCCCCGCTGACTc                                                                                                                        <  4:30354/149‑1 (MQ=255)                 tCGACCGCGATATCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCa                                                                                                                       <  5:225995/149‑1 (MQ=255)                    aCCGCGATATCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCAcc                                                                                                                     >  5:9295/1‑148 (MQ=255)                      cgcgATATCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCaa                                                                                                                  <  1:180739/149‑1 (MQ=255)                       gcgATATCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGTTCACTGCGCTGGCATTTTTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAg                                                                                                                 <  1:298985/149‑1 (MQ=255)                        cgATATCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGa                                                                                                                <  2:207256/149‑1 (MQ=255)                            atCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGaacaa                                                                                                            <  6:9295/149‑1 (MQ=255)                             tCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGaacaa                                                                                                            <  7:1719/148‑1 (MQ=255)                             tCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCGCACCGCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGaacaac                                                                                                           <  7:149392/149‑1 (MQ=255)                                 tGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGaacaacaa                                                                                                         >  2:196277/1‑147 (MQ=255)                                 tGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAacaacaaca                                                                                                       <  1:274562/149‑1 (MQ=255)                                 tGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAacaacaaca                                                                                                       <  7:336834/149‑1 (MQ=255)                                      tCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGt                                                                                                   >  2:9883/1‑148 (MQ=255)                                        cAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAAc                                                                                                >  4:284105/1‑149 (MQ=255)                                          gAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACca                                                                                              <  7:248748/149‑1 (MQ=255)                                            aaGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACcatc                                                                                            <  1:161998/149‑1 (MQ=255)                                            aaGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACcatc                                                                                            <  3:254249/149‑1 (MQ=255)                                            aaGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTACTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACcatc                                                                                            <  2:325847/149‑1 (MQ=255)                                             aGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACcatca                                                                                           >  3:265569/1‑149 (MQ=255)                                                 aCCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTc                                                                                       >  3:107590/1‑149 (MQ=255)                                                              ccGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTc                                                                          >  2:358698/1‑149 (MQ=255)                                                                  gCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCCCCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTg                                                                      >  2:109290/1‑149 (MQ=255)                                                                     cACCTCTCTGGTGATCACTGCTCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCGGGGCGCCGAGCGGCGccct                                                                                                                                <  8:89038/88‑2 (MQ=255)                                                                     cACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCGGGGCGCCGAGCGGCGccct                                                                                                                                >  7:89038/1‑87 (MQ=255)                                                                      aCCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGc                                                                  >  5:288709/1‑149 (MQ=255)                                                                               gTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTTCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTtgtg                                                          >  1:66552/1‑149 (MQ=255)                                                                               gTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTgg                                                         <  1:9883/149‑1 (MQ=255)                                                                               gTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCCGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGAtt                                                                                 <  5:131356/125‑1 (MQ=255)                                                                               gTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCCGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGAtt                                                                                 >  6:131356/1‑125 (MQ=255)                                                                                 gATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCtttt                                                                           <  6:374456/129‑1 (MQ=255)                                                                                 gATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCtttt                                                                           >  5:374456/1‑129 (MQ=255)                                                                                  aTCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGcc                                                       <  6:332011/148‑1 (MQ=255)                                                                                   tCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCctg                                                     <  4:107590/149‑1 (MQ=255)                                                                                    cACTGCGCTGGCACTGCGGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGc                                                                  <  8:131839/135‑1 (MQ=255)                                                                                       tGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACACAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCtgt                                                 >  1:300918/1‑149 (MQ=255)                                                                                         cgcTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCtgt                                                 <  3:324424/147‑1 (MQ=255)                                                                                               cACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCggctgg                                         <  3:284105/149‑1 (MQ=255)                                                                                                aCTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCggctggc                                        >  2:174396/1‑149 (MQ=255)                                                                                                  tgctgACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCTCCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGgctggctg                                      >  8:146634/1‑149 (MQ=255)                                                                                                       aCGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCt                                   <  1:358698/147‑1 (MQ=255)                                                                                                         cggCGCCTTG‑TCCTTCATGGCTTGAAGCCCCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAACAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCa                                                                <  8:317382/115‑1 (MQ=255)                                                                                                                      cTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTGccgccg                   >  6:360546/1‑149 (MQ=255)                                                                                                                      cTTCATGGTTTGAAGCACCCAGGGCGCCGCGCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTgccgcc                    <  2:152695/148‑1 (MQ=255)                                                                                                                       ttCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGa                                                                                                                >  7:305455/1‑55 (MQ=255)                                                                                                                       ttCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGa                                                                                                                <  8:305455/55‑1 (MQ=255)                                                                                                                       ttCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTGCCGCCgg                  >  2:329856/1‑149 (MQ=255)                                                                                                                         cATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCg                            <  7:292240/137‑1 (MQ=255)                                                                                                                         cATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCg                            >  8:292240/1‑137 (MQ=255)                                                                                                                          aTGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTGCCGCCGGttt               <  5:360546/149‑1 (MQ=255)                                                                                                                           tGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTGCCGCCGGttt               >  8:179736/1‑148 (MQ=255)                                                                                                                           tGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTGCCGCCGGtt                <  2:216592/147‑1 (MQ=255)                                                                                                                           tGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTGCCGCCGGtt                >  1:216592/1‑147 (MQ=255)                                                                                                                           tGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTGCCGCCGGTTTc              <  1:174396/149‑1 (MQ=255)                                                                                                                             gTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTGCCGCCGGTTTCAt            <  1:196277/149‑1 (MQ=255)                                                                                                                                   aGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTGCCGCCGGTTTCATCGAAGa      >  8:224900/1‑149 (MQ=255)                                                                                                                                        ccAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGGTGGCTGCGTGTGGCCGGTGGGGGCGGTTTGATCGAAGACAgg  >  6:199578/1‑147 (MQ=255)                                                                                                                                          ggggCGCCGAGCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTGCCGCCGGTTTCATCGaa        <  8:118629/139‑1 (MQ=255)                                                                                                                                          ggggCGCCGAGCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTGCCGCCGGTTTCATCGaa        >  7:118629/2‑140 (MQ=255)                                                                                                                                                    |                                                                                                                                        CTGCTGGTGGGCATGGTCGACCGCGATATCGGTGATTTCCAGAAAGCCACCGTGAAATGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTG‑TCCTTCATGGTTTGAAGCACCCAGGGCGCCGCTCGGCGCCCCGCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGCCACTTGTGGCCTGCTGTGCGGCTGGCTGCCTGTGGCCGGTGCCGCCGGTTTCATCGAAGACAGC  >  NC_002947/2341957‑2342238

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