Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,163,435	T→A	6.3%	I99N (ATT→AAT)	PP_5482 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,163,435	0	T	A	6.3%	191.6 / 3.1	63	I99N (ATT→AAT)	PP_5482	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/27);  new base A (2/2);  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 = 9.48e-01

TAGTGGGGCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCAGTAAATTTTG  >  NC_002947/2163291‑2163579                                                                                                                                                  |                                                                                                                                                  tAGTGGGGCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGAt                                                                                                                                               >  3:117088/1‑149 (MQ=255)    tGGGGCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGcc                                                                                                                                            <  1:417278/149‑1 (MQ=255)    tGGGGCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGcc                                                                                                                                            >  3:421108/1‑149 (MQ=255)        gCAACTGGTGCCTGGACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCttt                                                                                                                                         >  7:368484/1‑148 (MQ=255)        gCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCttt                                                                                                                                         >  2:264605/1‑148 (MQ=255)        gCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCttt                                                                                                                                         >  2:270837/1‑148 (MQ=255)          aaCTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCAATGATGCCttcata                                                                                                                                      <  7:60210/149‑5 (MQ=255)          aaCTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTAGAGTATGAAGGCATCATTGATGCCTTTCAg                                                                                                                                      <  4:263283/149‑1 (MQ=255)             tGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCta                                                                                                                                   >  6:396161/1‑149 (MQ=255)                tGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTatat                                                                                                                                >  2:356277/1‑149 (MQ=255)                      cACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTgg                                                                                                                            <  3:59244/147‑1 (MQ=255)                      cACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTgg                                                                                                                            <  4:117088/147‑1 (MQ=255)                       aCTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGAcc                                                                                                                         <  5:272350/149‑1 (MQ=255)                        cTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCa                                                                                                                        <  1:264605/149‑1 (MQ=255)                        cTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCa                                                                                                                        <  1:270837/149‑1 (MQ=255)                         tCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCAt                                                                                                                       >  5:99803/1‑149 (MQ=255)                           gAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATgg                                                                                                                     >  2:349146/1‑149 (MQ=255)                           gAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATg                                                                                                                      <  6:99803/148‑1 (MQ=255)                                cTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGAc                                                                                                                 <  1:49264/148‑1 (MQ=255)                                cTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTCGAACATTGCgcc                                                                                                                 >  4:359997/1‑146 (MQ=255)                                 tttAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATGGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTc                                                                                                               >  6:50139/1‑149 (MQ=255)                                         ggTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCtat                                                                                                                                  >  4:389484/1‑122 (MQ=255)                                         ggTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCtat                                                                                                                                  <  3:389484/122‑1 (MQ=255)                                         ggTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGt                                                                                                       <  5:217236/149‑1 (MQ=255)                                               cATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTg                                                                                                  <  1:216203/148‑1 (MQ=255)                                                aTTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATGATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGa                                                                                                >  6:269373/1‑149 (MQ=255)                                                aTTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGa                                                                                                >  7:321255/1‑149 (MQ=255)                                                   gAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGATTTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATAc                                                                                             >  6:316450/1‑149 (MQ=255)                                                    aGCGAAATAACGATGAGCGTGTTATCGCAGCCATCTAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCGTTTCAGCTATATGTGGACCATGGCGACTCGCCGTCGTTGCTGGATAc                                                                                             >  6:245148/1‑148 (MQ=255)                                                    aGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATAc                                                                                             <  8:368484/148‑1 (MQ=255)                                                     gCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTATGGCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATAc                                                                                             >  1:84154/1‑148 (MQ=255)                                                     gCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCa                                                                                           <  8:321255/149‑1 (MQ=255)                                                       gAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAgg                                                                                         <  2:84154/149‑1 (MQ=255)                                                           tAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGt                                                                                     >  5:387541/1‑149 (MQ=255)                                                             aCGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGATCTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGtag                                                                                   >  5:401766/1‑147 (MQ=255)                                                              cGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTtgt                                                                                  <  7:300310/149‑1 (MQ=255)                                                               gATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGTCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTCCTGGATACCAGTTCGTtgtg                                                                                 >  3:24921/1‑149 (MQ=255)                                                                     cGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAgg                                                                                         >  8:242672/1‑135 (MQ=255)                                                                            tCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGc                                                                     <  5:375873/148‑1 (MQ=255)                                                                            tCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCa                                                                    >  8:64514/1‑149 (MQ=255)                                                                                       aaaaTTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCg                                                          >  3:204009/1‑148 (MQ=255)                                                                                        aaaTTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTaa                                                               <  3:215987/142‑1 (MQ=255)                                                                                        aaaTTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTaa                                                               >  4:215987/1‑142 (MQ=255)                                                                                        aaaTTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCg                                                          >  5:282039/1‑147 (MQ=255)                                                                                         caTTAT‑GCTTGACGTGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAATGCATCATTGATGCCTTTCATCTATCTGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCTCCGCAGACGCACCTAAATTCGCA‑g                                                      <  8:79172/148‑1 (MQ=255)                                                                                         aaTTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑g                                                      <  5:396161/149‑1 (MQ=255)                                                                                                     aCGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGa                                            >  7:44073/1‑148 (MQ=255)                                                                                                     aCGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGGCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGa                                            <  3:359997/148‑1 (MQ=255)                                                                                                     aCGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCCCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGa                                            <  1:349146/148‑1 (MQ=255)                                                                                                                  gCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGaaa                              <  6:387541/149‑1 (MQ=255)                                                                                                                  gCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGa                                >  3:137762/1‑147 (MQ=255)                                                                                                                  gCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGa                                >  2:263398/1‑147 (MQ=255)                                                                                                                     ggAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGAAAGtt                           >  8:47227/1‑149 (MQ=255)                                                                                                                          tCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGAAAGTTTCAAc                      <  5:50139/149‑1 (MQ=255)                                                                                                                          tCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGAAAGTTTCAAc                      <  7:64514/149‑1 (MQ=255)                                                                                                                          tCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTACTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGAAAGTTTCAAc                      >  2:99115/1‑149 (MQ=255)                                                                                                                            gAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATAcc                                                                                            >  7:307396/1‑78 (MQ=255)                                                                                                                            gAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATAcc                                                                                            <  8:307396/78‑1 (MQ=255)                                                                                                                                tCGAGTATGAAGGCATCATTGATGGCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTa                 >  7:60253/1‑148 (MQ=255)                                                                                                                                    gTATGAAGTCATCATTGATGCCTTTCAGCTATATCTGGACCATGGCGACTCGGATTCGTTGCTGGATACGATGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTacac              <  5:245148/147‑1 (MQ=255)                                                                                                                                      tgaaacgCATCAATGATGCCTTTCATCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCACGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCACGCGGCATTTGAAACTACGGTTGaa                               <  7:266944/123‑1 (MQ=255)                                                                                                                                      tgaaaGGCATCAATGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGaa                               >  8:266944/4‑128 (MQ=255)                                                                                                                                      tgaaaGGCATCAATGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCAg           >  7:294010/4‑148 (MQ=255)                                                                                                                                              atcatTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCAGTAAATTTTg  <  4:204009/149‑1 (MQ=255)                                                                                                                                              atcatTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCAGTAAATTTTg  >  6:51022/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TAGTGGGGCAACTGGTGCCTGCACTCGAATGCTTTAAAACGGTAGTCATTGAGCGAAATAACGATGAGCGTGTTATCGCAGCCATCAAAATTAT‑GCTTGACGAGGCTCGTAGGCAGGAACTCGAATTCGAGTATGAAGGCATCATTGATGCCTTTCAGCTATATGTGGACCATGGCGACTCGCATTCGTTGCTGGATACCAGGTCGTTGTGCACCGCAGCCGCACCTAAATTCGCA‑GCGGCATTTGAAACTACGGTTGAAAGTTTCAACGCTTACACCAGTAAATTTTG  >  NC_002947/2163291‑2163579

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