Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,794,015	G→A	*427* (TAG→TAA)	gabT →	4‑aminobutyrate aminotransferase, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,794,015	0	G	A	100.0%	257.7 / NA	78	*427* (TAG→TAA)	gabT	4‑aminobutyrate aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (42/36);  total (42/36)

ATTCTTCTCTCCTGCGGCCCGTATTACAACGTGCTGCGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGT  >  NC_000913/2793890‑2794142                                                                                                                               |                                                                                                                                aTTCTTCTCTCCTGCGGCCCGTATTACAACGTGCTGCGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATg                                                                                                                    >  2:41076/1‑139 (MQ=255)          tcCTGCGGCCCGTATTACAACGTGCTGCGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCg                                                                                                           >  1:832866/1‑139 (MQ=255)           ccTGCGGCCCGTATTACAACGTGCTGCGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGa                                                                                                          >  1:752252/1‑139 (MQ=255)                             aCGTGCTGCGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCg                                                                                        <  2:1730830/139‑1 (MQ=255)                                tgctgcGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGAATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGc                                                                                      >  2:1718839/1‑139 (MQ=255)                                tgctgcGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGcc                                                                                     >  1:1378868/1‑139 (MQ=255)                                 gctgcGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTc                                                                                               <  1:550119/128‑1 (MQ=255)                                 gctgcGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTc                                                                                               >  2:550119/1‑128 (MQ=255)                                  ctgcGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTa                                                                                   >  2:1473187/1‑139 (MQ=255)                                    gcgcATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACg                                                                                 <  2:1378868/139‑1 (MQ=255)                                    gcgcATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACg                                                                                 <  2:1499443/139‑1 (MQ=255)                                     cgcATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGGCGCTGCGCGTCTTGTCCGGCCTACgg                                                                                >  2:1398789/1‑139 (MQ=255)                                          ccTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATc                                                                           <  2:264182/139‑1 (MQ=255)                                          ccTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATc                                                                           >  1:836946/1‑139 (MQ=255)                                          ccTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTC‑‑GGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGTCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCa                                                                         >  2:1181487/1‑139 (MQ=255)                                               tACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGt                                                                      >  1:852017/1‑139 (MQ=255)                                                    cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAacgc                                                                                                                              <  1:1327828/78‑1 (MQ=255)                                                    cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAacgc                                                                                                                              >  2:1327828/1‑78 (MQ=255)                                                    cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGt                                                                                                >  2:70459/1‑108 (MQ=255)                                                    cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGt                                                                                                >  1:207292/1‑108 (MQ=255)                                                    cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGt                                                                                                <  1:70459/108‑1 (MQ=255)                                                    cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGt                                                                                                <  2:207292/108‑1 (MQ=255)                                                    cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGAt                                                                 <  1:705900/139‑1 (MQ=255)                                                    cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGAt                                                                 >  1:1377768/1‑139 (MQ=255)                                                      cACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTc                                                                                          <  1:608143/112‑1 (MQ=255)                                                      cACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTc                                                                                          >  2:608143/1‑112 (MQ=255)                                                      cACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTAc                                                                                  <  2:1058349/120‑1 (MQ=255)                                                      cACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTAc                                                                                  >  1:1058349/1‑120 (MQ=255)                                                                  cGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgcc                                                   <  2:852017/139‑1 (MQ=255)                                                                  cGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgcc                                                   <  1:1578884/139‑1 (MQ=255)                                                                         aTCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCg                                            <  2:149650/139‑1 (MQ=255)                                                                           ccGTCAGGGTCTGGAGAATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                           <  1:1718839/139‑1 (MQ=255)                                                                            cGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcg                                                     <  2:1737994/127‑1 (MQ=255)                                                                            cGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcg                                                     >  1:1737994/1‑127 (MQ=255)                                                                                 gggTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAgg                                                                       <  2:886072/104‑1 (MQ=255)                                                                                 gggTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAgg                                                                       >  1:886072/1‑104 (MQ=255)                                                                                   gTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTc                                                                                          >  2:652019/1‑83 (MQ=255)                                                                                   gTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTc                                                                                          <  1:652019/83‑1 (MQ=255)                                                                                   gTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATaa                                                               <  2:1512960/110‑1 (MQ=255)                                                                                   gTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATaa                                                               >  1:1512960/1‑110 (MQ=255)                                                                                   gTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                           <  1:349576/130‑1 (MQ=255)                                                                                   gTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                           >  2:349576/1‑130 (MQ=255)                                                                                     cTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATaa                                                               <  1:840766/108‑1 (MQ=255)                                                                                     cTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATaa                                                               >  2:840766/1‑108 (MQ=255)                                                                                        gag‑aTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGa                             >  1:1026786/1‑139 (MQ=255)                                                                                         ag‑aTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGaa                            <  2:752252/139‑1 (MQ=255)                                                                                             tcatcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCg                                                                                        <  2:962715/76‑1 (MQ=255)                                                                                             tcatcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCg                                                                                        >  1:962715/1‑76 (MQ=255)                                                                                             tcatcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACgggg                                                                              >  2:805499/1‑86 (MQ=255)                                                                                             tcatcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACgggg                                                                              <  1:805499/86‑1 (MQ=255)                                                                                             tcatcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAAt                                       >  2:394769/1‑125 (MQ=255)                                                                                             tcatcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAAt                                       <  1:394769/125‑1 (MQ=255)                                                                                              catcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCCCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAgg                        >  1:471646/1‑139 (MQ=255)                                                                                               atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCAt                                               <  2:648359/115‑1 (MQ=255)                                                                                               atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCAt                                               >  1:648359/1‑115 (MQ=255)                                                                                               atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCg                                            <  2:1577203/118‑1 (MQ=255)                                                                                               atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCg                                            >  1:1577203/1‑118 (MQ=255)                                                                                               atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGga                       >  1:749007/1‑139 (MQ=255)                                                                                               atcaGCCAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGga                       >  2:246958/1‑139 (MQ=255)                                                                                                    ccAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                      <  2:504794/119‑1 (MQ=255)                                                                                                    ccAGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                      >  1:504794/1‑119 (MQ=255)                                                                                                      aGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACg                                                                                 <  1:1300699/74‑1 (MQ=255)                                                                                                      aGTGTTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACg                                                                                 >  2:1300699/1‑74 (MQ=255)                                                                                                        tgtTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgcc                                                   >  1:12506/1‑102 (MQ=255)                                                                                                        tgtTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgcc                                                   <  2:12506/102‑1 (MQ=255)                                                                                                        tgtTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTa              >  2:316718/1‑139 (MQ=255)                                                                                                         gtTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGAt                                                                            <  2:446884/76‑1 (MQ=255)                                                                                                         gtTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGAt                                                                            >  1:446884/1‑76 (MQ=255)                                                                                                         gtTTTGATGAGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTAt             <  2:836946/139‑1 (MQ=255)                                                                                                                  aGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                           <  1:694785/100‑1 (MQ=255)                                                                                                                  aGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                           >  2:694785/1‑100 (MQ=255)                                                                                                                  aGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                      <  1:1673122/121‑1 (MQ=255)                                                                                                                  aGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                      >  2:1673122/1‑121 (MQ=255)                                                                                                                  aGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                      <  1:1573500/121‑1 (MQ=255)                                                                                                                  aGGCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                      >  2:1573500/1‑121 (MQ=255)                                                                                                                    gCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGt                                                                                                >  1:1827745/1‑45 (MQ=255)                                                                                                                    gCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGt                                                                                                <  2:1827745/45‑1 (MQ=255)                                                                                                                    gCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcg                                                     <  2:1054693/88‑1 (MQ=255)                                                                                                                    gCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcg                                                     >  1:1054693/1‑88 (MQ=255)                                                                                                                    gCGAAGCAGTAACGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGt  >  2:471737/1‑139 (MQ=255)                                                                                                                               |                                                                                                                                ATTCTTCTCTCCTGCGGCCCGTATTACAACGTGCTGCGCATCCTTGTACCGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAG‑ATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGT  >  NC_000913/2793890‑2794142

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑