Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,295,021	T→A	11.1%	intergenic (+37/+33)	slyB → / ← nhaA‑I	outer membrane lipoprotein/Na(+)/H(+) antiporter NhaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,295,021	0	T	A	11.1%	162.3 / 9.4	63	intergenic (+37/+33)	slyB/nhaA‑I	outer membrane lipoprotein/Na(+)/H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base A (4/3);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GATGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTG  >  NC_002947/1294874‑1295165                                                                                                                                                    |                                                                                                                                                  gatgGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGctgctg                                                                                                                                                  >  3:325003/1‑149 (MQ=255)  atgGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCtgctgc                                                                                                                                                 <  5:292372/149‑1 (MQ=255)          aTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCggg                                                                                                                                        <  6:245248/149‑1 (MQ=255)           tGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGcc                                                                                                                                                <  7:185647/140‑1 (MQ=255)           tGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGcc                                                                                                                                                >  8:185647/1‑140 (MQ=255)              gcgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCGGGGttt                                                                                                                                    >  8:71980/1‑149 (MQ=39)                 ccTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCGGGGtttttt                                                                                                                                 <  7:71980/149‑1 (MQ=39)                  cTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGttgttt                                                                                                                                 >  2:246215/1‑148 (MQ=255)                        gCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGAGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGAGCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCt                                                                                                                          >  8:169135/1‑149 (MQ=255)                        gCAGGCCGTGCAGAAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCt                                                                                                                          >  1:200017/1‑149 (MQ=255)                         cAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGTCTGTAGACGGTACCAGTCGCGTTTCGCACTGAACGTAGTGCAGTAATAAAAAACCCCGACATGTCTTCTGCCA‑GTCGGGGTTTTTTCATGCCTc                                                                                                                         <  8:317105/149‑1 (MQ=255)                          aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCa                                                                                                                        <  1:251860/149‑1 (MQ=255)                          aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCa                                                                                                                        >  2:251860/1‑149 (MQ=255)                            gCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGc                                                                                                                      >  2:49445/1‑149 (MQ=255)                             ccGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGcc                                                                                                                     <  5:244978/149‑1 (MQ=255)                                     gagaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAgg                                                                                                             <  8:305201/149‑1 (MQ=255)                                     gagaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAgg                                                                                                             >  1:265155/1‑149 (MQ=255)                                     gagaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTAGCACTGACCGTAGTGCAGTAATAAAAAACCCCCACATGGCTGCTGCCA‑GTCGGGGTTTTTTAATTCCTCAGCCAATTCAgg                                                                                                             >  2:317763/1‑149 (MQ=255)                                      agaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGc                                                                                                            <  4:325003/149‑1 (MQ=255)                                           gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGttt                                                                                                                                    >  2:134752/1‑120 (MQ=255)                                           gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGttt                                                                                                                                    <  1:134752/120‑1 (MQ=255)                                                 ttCCGTTTTGGAGACCGACTGCCCATCATCACAGTATACGGTACAAGTCGCGATTCGCAATTACCGTAGTGCAGTAATAAAAAACCCCGAAATGAATGCTGACA‑GTCAGGGTTTTTTCATGCCTCAGCCAAGTCATGCAATGCTGGttt                                                                                                 >  4:186770/1‑149 (MQ=255)                                                     gtgtTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCg                                                                                             >  4:144035/1‑149 (MQ=255)                                                      tgtTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCgg                                                                                            <  1:49445/149‑1 (MQ=255)                                                           gCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTgg                                                                                        <  1:246215/148‑1 (MQ=255)                                                                 gAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGc                                                                                 >  2:128291/1‑149 (MQ=255)                                                                    tGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGt                                                                              <  5:112381/149‑1 (MQ=255)                                                                         atcatGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACgg                                                                         >  3:9355/1‑149 (MQ=255)                                                                          tcatGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGc                                                                        >  2:183045/1‑149 (MQ=255)                                                                           catGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTg                                                                                  <  3:25689/138‑1 (MQ=255)                                                                           catGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTg                                                                                  >  4:25689/1‑138 (MQ=255)                                                                              gACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTc                                                                                                                                <  7:21878/89‑1 (MQ=255)                                                                              gACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTc                                                                                                                                >  8:21878/1‑89 (MQ=255)                                                                              gACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAg                                                                    >  6:229600/1‑149 (MQ=255)                                                                               aCTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAg                                                                    <  5:30/148‑1 (MQ=255)                                                                                   tAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCg                                                                 <  2:265155/147‑1 (MQ=255)                                                                                     gACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCa                                                             >  7:130917/1‑149 (MQ=255)                                                                                      aCGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAc                                                            >  1:50512/1‑149 (MQ=255)                                                                                       cGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAcc                                                           <  2:50512/149‑1 (MQ=255)                                                                                            ccAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCgcag                                                       <  1:128291/148‑1 (MQ=255)                                                                                             cAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCgcagc                                                      >  6:244248/1‑148 (MQ=255)                                                                                                   cgTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaac                                               >  7:300364/1‑149 (MQ=255)                                                                                                    gTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCaa                                                                                                                   >  2:261279/1‑80 (MQ=255)                                                                                                    gTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCaa                                                                                                                   <  1:261279/80‑1 (MQ=255)                                                                                                       tCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaa                                            <  5:244248/148‑1 (MQ=255)                                                                                                               aCCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAg                                    <  1:317763/148‑1 (MQ=255)                                                                                                                   tAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGc                               <  8:130917/149‑1 (MQ=255)                                                                                                                   tAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGc                               >  3:222001/1‑149 (MQ=255)                                                                                                                        cAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCAt                           <  3:144035/148‑1 (MQ=255)                                                                                                                         aGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGGTGGTTTTGCGGATGGCCATTGCGGTGACGGAGTATCCGATCACAGCAGCAAAGAACGAACCAGTAAGAATGCCCATGc                         >  1:115902/1‑149 (MQ=255)                                                                                                                              taaAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                           >  8:110498/1‑126 (MQ=255)                                                                                                                              taaAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                           <  7:110498/126‑1 (MQ=255)                                                                                                                              gaaaaaaCCCCGAC‑TGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTc                     >  8:305233/2‑148 (MQ=38)                                                                                                                               aaaaaaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCa                   >  4:88924/1‑149 (MQ=255)                                                                                                                                aaaaaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCAt                  <  4:9355/149‑1 (MQ=255)                                                                                                                                 aaaaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTGTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATg                 >  3:212316/1‑149 (MQ=255)                                                                                                                                   aaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAAt                                 <  7:255708/131‑1 (MQ=255)                                                                                                                                   aaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAAt                                 >  8:255708/1‑131 (MQ=255)                                                                                                                                     ccccGAC‑TGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCgg             >  5:79437/1‑149 (MQ=255)                                                                                                                                      cccGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCAt                  >  6:210835/1‑143 (MQ=255)                                                                                                                                      cccGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCAt                  <  5:210835/143‑1 (MQ=255)                                                                                                                                      cccGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGc            >  8:311143/1‑149 (MQ=255)                                                                                                                                         cgactGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGt          <  7:305233/145‑1 (MQ=255)                                                                                                                                           cATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGc            <  6:354217/144‑1 (MQ=255)                                                                                                                                           cATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGc            >  5:354217/1‑144 (MQ=255)                                                                                                                                               gctgctGCCA‑GTCGTGTTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTTGTTTTGCGTCTGGCCATTGCGGTGACGGCTTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCAc    <  2:115902/148‑1 (MQ=255)                                                                                                                                               gctgctGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCAc    <  4:222001/148‑1 (MQ=255)                                                                                                                                               gctgctGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACt   >  5:54007/1‑149 (MQ=255)                                                                                                                                               gctgctGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACt   <  1:120163/149‑1 (MQ=255)                                                                                                                                                ctgctgCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTg  >  3:299498/1‑149 (MQ=255)                                                                                                                                                 tgctgcCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGc            <  7:343/138‑1 (MQ=255)                                                                                                                                                 tgctgcCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGc            >  8:343/1‑138 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GATGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTG  >  NC_002947/1294874‑1295165

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