Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	835,912	0	C	G	40.3%	92.1 / 14.4	67	A222G (GCG→GGG)	ybiB	putative family 3 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (13/27);  new base G (27/0);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 3.33e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGATGGGCGTGCGTAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGG  >  NC_000913/835767‑836061                                                                                                                                                  |                                                                                                                                                      cGGATGGGCGTGCGTAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGctac                                                                                                                                                  <  2:759677/151‑1 (MQ=255)  ggATGGGCGTGCGTAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGctac                                                                                                                                                  >  2:168955/1‑150 (MQ=255)  ggATGGGCGTGCGTAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGctac                                                                                                                                                  <  1:1071595/150‑1 (MQ=255)     tGGGCGTGCGTAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGccacta                                                                                                                                                >  1:1255316/1‑149 (MQ=255)     tGGGCGTGCGTAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCGACTaa                                                                                                                                               >  1:944075/1‑150 (MQ=255)     tGGGCGTGCGTAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCCCCTAAt                                                                                                                                              >  2:571935/1‑151 (MQ=255)      gggCGTGCGTAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAAt                                                                                                                                              <  1:608357/150‑1 (MQ=255)      gggCGTGCGTAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAAt                                                                                                                                              <  2:790585/150‑1 (MQ=255)      gggCGTGCGTAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATg                                                                                                                                             <  1:428629/151‑1 (MQ=255)        gcgtgcgtAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCa                                                                                                                                           <  2:1090286/151‑1 (MQ=255)          gtgcgtAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGGGGGGGCTACTAAAGCAt                                                                                                                                          >  2:63257/1‑150 (MQ=255)           tgcgtAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGGGGGCGCGTCTAATGCATgg                                                                                                                                        >  2:442573/1‑151 (MQ=255)            gcgtAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGGTGCTAATgaa                                                                                                                                           >  1:722268/1‑145 (MQ=255)            gcgtAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATgg                                                                                                                                        <  2:1136113/150‑1 (MQ=255)             cgtAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGGGGGCGCCACTAATGCATg                                                                                                                                         >  1:1244894/1‑148 (MQ=255)                aaCAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTg                                                                                                                                   <  1:34383/151‑1 (MQ=255)                 aCAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGa                                                                                                                                  <  2:468440/151‑1 (MQ=255)                  cAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCCACTTATGCATGGTAc                                                                                                                                     >  1:132494/1‑147 (MQ=255)                   aGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGccacttaagacatg                                                                                                                                        >  1:479287/1‑134 (MQ=255)                   aGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCGGCCCAAGCATGGTACTg                                                                                                                                   >  1:695971/1‑148 (MQ=255)                       cacaTACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAgggg                                                                                                                             >  2:1328036/1‑148 (MQ=255)                             cccTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCCACCAATGCATGGTACTGAAgggggag                                                                                                                          >  1:154710/1‑142 (MQ=255)                              ccTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGGGGGGGCTACTAATGCCTGGTACTGAAgggg                                                                                                                             >  2:541484/1‑141 (MQ=255)                                          ggCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCTACTAATGCATGGGACTGAAGGGGGAGTGTATGCTAATcc                                                                                                             >  2:898343/1‑147 (MQ=255)                                          ggCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGc                                                                                                           <  2:961162/149‑1 (MQ=255)                                               cGCCATTTGCCGAAGGGGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCGACTTATTGATGGTCCTGAGGGGGAAGGGTTTGCTGAtttgccacg                                                                                                       >  1:977385/1‑140 (MQ=255)                                                    tttGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCTACTAATGCATGGTACTGGAGGGGGAGGGTATGCTAATCCCCAGCGCTGCCCGc                                                                                               >  2:178678/1‑151 (MQ=255)                                                       gCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCCACTAATGCATGGTACTGAAAGGGAAGGGTATGCTAATCCCCCGCGCTGccccc                                                                                               >  2:871461/1‑146 (MQ=255)                                                        ccGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGGGGGGGCGACTAATGGATGGGAATGAAGGGGGAGTGGATGCTAATCCCCCGCGCCGccc                                                                                                 >  1:748361/1‑145 (MQ=255)                                                        ccGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGAt                                                                                           <  1:633402/151‑1 (MQ=255)                                                         cGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGGTTCTAATGCATGGTACTGAAGGGGAAGGGGATGCTAATCCGCCGCGCTGCCCGc                                                                                               >  2:549615/1‑146 (MQ=255)                                                           aaGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGGGGGGGCTACTAATGCCTGGGACTGAAGGGGGAGTGGATGCTAATACGCCGCGCTGCCCGCAGATCa                                                                                         >  2:732209/1‑150 (MQ=255)                                                            aGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCCACTTATGCATGGTACTTAAGGGGAAGTGGAGGCTAATCCCCAGCGCCGCCCCCCGATcca                                                                                        >  1:884966/1‑148 (MQ=255)                                                               tGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCTACTAAAGCATGGGACTGTAGGGGAAGGGGATGCTAATTCGCCGCGCTGCCCGCAGATCAAtct                                                                                     >  2:447401/1‑150 (MQ=255)                                                                 aGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCAt                                                                                  <  2:1076487/151‑1 (MQ=255)                                                                 aGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCAt                                                                                  <  1:125210/151‑1 (MQ=255)                                                                            cTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATTCCTGGGACTGAAGGGGGAGGGGATGCTAATCCCCCGCGCTGCCCGCCGATCAAACTCATTGTCCGTg                                                                          >  1:1078061/1‑148 (MQ=255)                                                                                ccAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGAt                                                                   <  2:402047/151‑1 (MQ=255)                                                                                ccAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCGACTAAAGCATGGTACTGAAGGGGGAGTGTATGCTAAACCCCCGCGCCGCCCCCAGAGCAATCTCATTTAACGCGGAgggg                                                                    >  1:992527/1‑149 (MQ=255)                                                                                   gCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCCGGCGGCACCAATACATGGTACTGAAGGTGGAGAGGATGTTAATCCCCAGCGGGGCCCCCACCTCCAACTCCATGTTCGTGAAggggagcg                                                                >  1:1088218/1‑146 (MQ=255)                                                                                        tCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCTACTTATGCATGGGACTGAAGGGGAAGGGGATGCTAATCCCCAGCGCTGCCCCCAGATCAATCTCCTTTACCGTGAAGGGGTGCggggg                                                            >  2:621086/1‑148 (MQ=255)                                                                                           cATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCCACTAAAGCATGGTACTGAAGGTGAAGGGGAAGCTAATCCGCAGCGCCGCCCCCCGATCAATCTCATTTACCCGGAAAGGGAGCggggg                                                            >  1:496871/1‑145 (MQ=255)                                                                                           cATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGt                                                        <  2:317036/151‑1 (MQ=255)                                                                                            aTCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCTACTAATGCATGGTACTGGAGGGGAAGGGGATGCTAATCCGCAGCGCTGCCCCCAGATCAAACTCCTTTACCGTGGAGGGGTGCggggg                                                            >  2:15534/1‑144 (MQ=255)                                                                                              ccGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTAt                                                      <  2:120840/150‑1 (MQ=255)                                                                                               cGGAATACAGTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTTACCGTGAAGGGATGCGGGTGCTGTATg                                                     <  1:421759/150‑1 (MQ=255)                                                                                                 gAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGAAGGCGCTCCTAAAAAATGGTACGTAAGGTGAAGGGGGTGCTAATCCGCAAGGCTGCCCGCAGGTCAATCTCATTGACCGTGGAAGGATTCggggggggtctgaa                                                   >  1:91797/1‑139 (MQ=255)                                                                                                  aaTACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTa                                                       <  1:1069502/145‑1 (MQ=255)                                                                                                       aTTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAgg                                            <  2:274096/151‑1 (MQ=255)                                                                                                        ttGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTACGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAgg                                            <  1:69514/150‑1 (MQ=255)                                                                                                         tGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGa                                           <  1:184821/150‑1 (MQ=255)                                                                                                         tGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACCAAAGCATGGTACTGAAGGGGGAGTGTGTGCTAAATCGCCGCGCTGCCCCCAGAACAATCTCATTGACCGTGGAGGGATGCGGGGGCTGTATGAAAAACAAGa                                           >  2:877168/1‑150 (MQ=255)                                                                                                          ggACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGaca                                         <  2:302214/151‑1 (MQ=255)                                                                                                               cgTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCTACTAATGCATGGTACTGAAGGGGAAGGGGATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACActgc                                     >  2:168427/1‑150 (MQ=255)                                                                                                               cgTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACActgct                                    <  1:1186473/151‑1 (MQ=255)                                                                                                                gTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGGTACTACTCCAGGGTACTGCTGGGGAAAGGGAGGGTTATATGCAGGGCTGCCCCCAGATCCATCTCATTGGCCGGTAAGGGGGGGGGGGGGTGTATGAAAAACa                                              >  1:233274/1‑140 (MQ=255)                                                                                                                  cgcgAAGTTCTTTAGCGATATCGGCGGGCGGGGGCCACTAAAGGATGGGACTGAAGGGGGAGGGGATGCTAATCCGCCGCGCTGCCCGCCGACCAATCTCATTTACCCCGAAGGGAAGGGGGGGGTGGATGAAAAAAAGGacaacgc                                     >  2:378326/1‑143 (MQ=255)                                                                                                                   gcgAAGTTCTTTAGCGATATCGGCGGGCGGGGGCCACTAATTCATGGTTCTGAAGGGGAAGGGTATGCTGAATCGCCGCCCTGCCCGCAGATCAATCTCACTGAGCGTGAAGGGATGCGGGTGGTGTATGaaaaaaaagaaacctctgg                                  >  1:672608/1‑135 (MQ=255)                                                                                                                     gAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGc                              <  2:951063/151‑1 (MQ=255)                                                                                                                     gAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGc                              <  1:864129/151‑1 (MQ=255)                                                                                                                     gAAGTTCTTTAGCGATATCGGCGGGCGGGCGCCACTAATGCATGGGACTGAAGGGGAAGTGTATGCTAAACCGCAGCGCTGCCCGCAGATCAAACTCATTTACCGTGGAGGGGTGCGGGGGCTGTATGGAAAACAAGACACTGCTgg                                  >  2:381378/1‑147 (MQ=255)                                                                                                                      aaGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCg                             <  2:1283707/151‑1 (MQ=255)                                                                                                                        gTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTTAAGGGGAAGTGTATGCTAATCCCCAGCGCCGCCCCCCGATCAATCTCCTTTACCGTGGAGGGGTGCGGGGGGTGTATGAAAAAAAAGACACCGCTGGGAg                               >  2:209689/1‑147 (MQ=255)                                                                                                                           cTTTAGCGATATCGGCGCGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGt                          <  1:97336/149‑1 (MQ=255)                                                                                                                             ttAGCGATATCGGCGGGCGGGGGCTACTAATGCATGGTACTGAAGGTGGAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATAAATCTCATTGACCGTGAAGGGGTGCGGGGGCTGTATGAAAAAAAGGACACTGCTGGTAGCGAGGTACt                      >  2:736924/1‑151 (MQ=255)                                                                                                                                  gATATCGGCGGGCGGGGGCTAATAATGCATGGTACTGAAGGGGAAGGGTATGCTAATCCGCAGCGCTGCCCGCACATCAATCTCATTGTCCGCGAAGGGAGGCGGGGGCTGTTTGAAAAAAAGGAAACTGCTGGGAGCGAGGTACTg                     >  1:859951/1‑147 (MQ=255)                                                                                                                                       cggcggGCCGGCGCTACCAATGCATGGTACTGGAGGTGCAGTGCATGCTAATACGCAGCGCTGCCCGCAGATCAATCTCAATGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTCCCGGTAGCGCGTTACTGCCACaa               >  1:107550/1‑148 (MQ=255)                                                                                                                                        ggcgggcgggcgCTACTAATGCATGGTACTGAAGGCGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGTCACTGCTGGTAGCGAGTTACTGCCACAAGCa            >  2:1211587/1‑150 (MQ=255)                                                                                                                                           gggcgggcgCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAg        <  1:1335414/151‑1 (MQ=255)                                                                                                                                                ggcgCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCg   >  1:51405/1‑151 (MQ=255)                                                                                                                                                 gcgcTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCgg  >  2:975391/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                      CGGATGGGCGTGCGTAACAGTGCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGCATGGTACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGG  >  NC_000913/835767‑836061

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