Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,575,956	0	A	C	40.8%	118.3 / 21.3	71	V369G (GTG→GGG)	gntU	gluconate transporter, low affinity GNT 1 system
Reads supporting (aligned to +/- strand):  ref base A (28/14);  new base C (0/29);  total (28/43)
Fisher's exact test for biased strand distribution p-value = 1.44e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGCTTCGGTCGCGCCGGTAAATTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGACCTACGCCAGAGTCAACC  >  NC_000913/3575809‑3576102                                                                                                                                                     |                                                                                                                                                   tCGGCTTCGGTCGCGCCGGTAAATTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCAcc                                                                                                                                                  <  1:881387/150‑1 (MQ=255) tCGGCTTCGGTCGCGCCGGTAAATTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCAcc                                                                                                                                                  >  2:720134/1‑150 (MQ=255)  cGGCTTCGGTCGCGCCGGTAAATTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAg                                                                                                                                                <  2:1065750/151‑1 (MQ=255)    gCTTCGGTCGCGCCGGTAAATTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTc                                                                                                                                              <  2:366083/151‑1 (MQ=255)          gTCGCGCCGGTAAATTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTAccg                                                                                                                                        <  2:492021/151‑1 (MQ=255)            cgcgCCGGTAAATTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTAccgcc                                                                                                                                      <  2:334845/151‑1 (MQ=255)                     aaaTTTACCAAACAACCAAAACCCGGCGTCGTTACCGGGGCTGACAACAATCGAACCACCAGCGATACAAACCGACAGCGCCGCCATTGGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCCCCAGTCCTACCGCCGTTAAACAg                                                                                                                             <  2:1029577/151‑1 (MQ=255)                      aaTTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAgg                                                                                                                            >  2:346715/1‑151 (MQ=255)                        tttACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCt                                                                                                                          >  1:1042678/1‑151 (MQ=255)                              gAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTg                                                                                                                    >  1:209390/1‑151 (MQ=255)                               aacaccCAGAACCCGGCGTCGTTAACGTGGCTGACAACAATCGACCCACCAGCGATAAAAATCGACAGCGCCGCCATTGGCGCACCGGAG‑TAGTTCAGTTTTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTgg                                                                                                                   <  1:396969/151‑1 (MQ=255)                               aacaacCAAAACCGGGCGTCGTTAACGTGGCTGACAAAAACCGAACCCCCAGCGATACAAACCGACAGCGCCCCCATTTGCGCCCGGGAG‑TAGTTCATTTTTTCAATACCCGGCATCCCCAGCCCTACCGCCGTTAAACAGGCTACGGTg                                                                                                                    <  2:861449/150‑1 (MQ=255)                               aaaaacCAGAAACCGGCGTCGTTAACGGGGCTGACAAAAATCGACCCCCCAGCGATACAAACCGACACCGCCCCCATTGGCGCACCGGAG‑TAGTTCAGTTGTTCAATACCCGGCATCCCCAGCCCTACCGCCGTTAAACAGGCTACGGTg                                                                                                                    <  2:471658/147‑1 (MQ=255)                                acaacCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGc                                                                                                                  >  1:246419/1‑151 (MQ=255)                                acaacCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGc                                                                                                                  >  1:1080328/1‑151 (MQ=255)                                      aaaacccGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCAATACAAATCGACAGCGCCGCCATTTGCGCCCCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAAc                                                                                                             <  1:764131/145‑1 (MQ=255)                                         accgggCGTCGTTAACGGGGCTGACAAAAACCGACCCCCCAGCAATAAAAATCGACAGCGCCCCCATTGGCGCACCGGAG‑TATTTCAGTTTTTCAATAACCGGCATCCCAATTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTg                                                                                                         <  2:1041812/147‑1 (MQ=255)                                          aCCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTg                                                                                                         >  1:459663/1‑150 (MQ=255)                                           ccGGCGTCTTAACCGTGGCTGACAACAATCGACCCCCCAGCGATAAAAGTCGACAGCGCCGCCATTTGCGCCCCGGAT‑TATTTCAGTTGTTCAATACCCGGCATCCCCATTCCTCCCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGa                                                                                                        <  1:217379/150‑1 (MQ=255)                                           ccGGCGTCGTTAACGTGGCTAACAACAACCGACCCCCCAGCAAAAAAAACCGACAGCCCCGCCATTTGCGCCCCGGAG‑TAGTTCAGTTGTTCAATACCCGCCACCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGaa                                                                                                       <  2:1019072/151‑1 (MQ=255)                                            cGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAat                                                                                                      >  2:595066/1‑151 (MQ=255)                                               cgtcgtTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAatgat                                                                                                   >  1:1018204/1‑151 (MQ=255)                                                gtcgtTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAatgatg                                                                                                  >  1:131457/1‑151 (MQ=255)                                                gtcgtTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAatgatg                                                                                                  >  1:774975/1‑151 (MQ=255)                                                              taaaaaaaaccGAACCCCCGGCAATACAAATCAACAGCCCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATACCGGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGct                                                                                         <  1:441859/136‑2 (MQ=255)                                                                  aacaaTCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGc                                                                                >  1:460888/1‑151 (MQ=255)                                                                  aacaaACGAACCACGAGCGAAAAAAATCGACGGCGCCGCCTTTTGCGCACGGGAG‑TATTGAAGTTGTTCAAGACCCGGCATCACAGGTCCTGCCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAg                                                                                 <  1:391080/150‑1 (MQ=255)                                                                     aaccaaaCCACCAGCGAAAAAAATCGACAGCCCCGCCATTTGCGCA‑CGGAGTTAGTTCAGTTTTTAAATACCCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCAc                                                                              <  2:1008198/145‑1 (MQ=255)                                                                     aaTCGAACCCCCAGCAAAAAAAACCGCCAGCCCCCCCATTGCCGCACCGGAG‑TATTCCAGTTTTTCAAAAACCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCa                                                                               <  2:271769/149‑1 (MQ=255)                                                                     aaTCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACg                                                                             >  1:824116/1‑151 (MQ=255)                                                                        cgcaccgccATCGATACAAACCGACAGCGCCCCTTTTTGCGCACCGGAG‑TATTTCAGTTTTTAAATACCCGGCTCCCCCTTCCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAg                                                                          <  2:839510/148‑1 (MQ=255)                                                                        cGACCCCCCAGCGATAAAAATCGACAGCGCCGCCATTTGCGCCCCGGAG‑TAGTTCAGTTGTTCAATACCCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGaa                                                                           <  1:346708/150‑1 (MQ=255)                                                                        cGACCCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCCCCGGAG‑TAGTTCAGTTGTTCAATACCCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAg                                                                          <  1:741254/151‑1 (MQ=255)                                                                          accccccaGCAAAAAAAATCGACGGCGCCGCCTTTGGCGCCCCGGAG‑TATTTCTGTTGTTCAATACCCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGaa                                                                           <  2:807130/143‑1 (MQ=255)                                                                           accaccAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGc                                                                         >  1:622892/1‑149 (MQ=255)                                                                              accaGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTg                                                                    >  2:823031/1‑151 (MQ=255)                                                                               ccaGCGATACAAACCGACGGCGCCGCCTTTGGCGCACCGGAG‑TAGTTCATTTGTTCAATACCCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGa                                                                   <  2:169706/151‑1 (MQ=255)                                                                                  gCGAAACAAACCAACAGCCCCCCCATTTGCGCCCCGGAG‑TAGTTCAGTTGTTCAATACCCGGCACCCCCATTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATa                                                                 <  2:159389/150‑1 (MQ=255)                                                                                   caaTACAAATCACCAGCCCCCCCATTTGCGCCCGGGAG‑TAGTTAGGTTGTAAAATACCCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATa                                                                 <  2:694894/147‑1 (MQ=255)                                                                                     aTACAAATCGACAGCGCCGCCATTTGCGCCCCGGAG‑TAGTTCAGTTGTTAAATACCCGGCACCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCaa                                                             <  1:869985/151‑1 (MQ=255)                                                                                       aCAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATc                                                           >  1:1024804/1‑151 (MQ=255)                                                                                         aaaTCCACGCCGCCCCCATTTGCGCACGGGAG‑TTGTTCATTTTTACAAAAACCGACACCCCCGGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATa                                                                 <  2:285535/143‑1 (MQ=255)                                                                                         aaaGCGCCAGCGCCGCCTTTTGCCCCCCGGGG‑TAGTTCATTTTTCAAATACCCGGCACCCCCGGTCCTCCCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATa                                                                 <  2:822614/143‑1 (MQ=255)                                                                                          aaccGACCGCGCCCCCATTTGCGCACCGGAG‑TATTTCATTTGTTCAATACCCGCCATCCCCATTCCTCCCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGc                                                        <  2:104649/148‑1 (MQ=255)                                                                                          aaTCAACAGCCCCCCCATTCGCGCACGGGAG‑TAGTTCAGTTGTTAAATACCCGACACCACAAGTCCTCCCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCg                                                          <  2:262765/149‑1 (MQ=255)                                                                                             cGACAGCGCCGCCTTTGGCGCACCGGAG‑TGGTTCATTTGTTAAAAAACCGGCACCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAg                                                      <  1:716287/150‑1 (MQ=255)                                                                                                   cccccccttttGCGCCCGGGAG‑TATTTCATTTTTTAAAAAACCGCCACCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAgg                                                     <  2:516330/137‑1 (MQ=255)                                                                                                     ccgccATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGcc                                             >  1:475215/1‑151 (MQ=255)                                                                                                     cccccATTTGCGCCCGGGAG‑TAGTTCATTTTTCAAATACCCGGCATCCCCAGCCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGcc                                             <  2:417830/148‑1 (MQ=255)                                                                                                        ccATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCgg                                           >  2:452862/1‑150 (MQ=255)                                                                                                           tttGCGCCCGGGAT‑TATTTCAGTTGTTCAATAACCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGtt                                         <  2:1069596/149‑1 (MQ=255)                                                                                                               cgcccgggAG‑TATTTCGGTTTTTCAATAACCGGCTTCCCCAGCCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAAc                                      <  2:840444/142‑1 (MQ=255)                                                                                                               cccACCGGAG‑AAGTTCAGTTTTTCAATACCCGGCACCCCCAGCCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCt                                   <  2:1080372/149‑1 (MQ=255)                                                                                                                gcACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCtt                                  >  1:670313/1‑151 (MQ=255)                                                                                                                 ccccGGAG‑TAGTTCAGTTGTTAAATACCCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCtt                                  <  1:939308/148‑1 (MQ=255)                                                                                                                 cACGGGAG‑TATTTCATTTGTTAAAAAACCGGCATCCCCGTCCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTc                                 <  2:416230/151‑1 (MQ=255)                                                                                                                  aCCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACGGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTc                                 >  1:1115699/1‑150 (MQ=255)                                                                                                                     ggAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGccc                             >  1:989609/1‑151 (MQ=255)                                                                                                                          tagtTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTg                         >  1:52954/1‑151 (MQ=255)                                                                                                                           agtTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGC‑GAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCCTCGCCCa                            >  2:836781/1‑146 (MQ=255)                                                                                                                               cAGTTGTTCAATAACCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCgg                     <  2:246435/150‑1 (MQ=255)                                                                                                                                 gttttccAATACCCGCTACCACCAGCCCAACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGAc                   <  2:638102/144‑1 (MQ=255)                                                                                                                                 gttgttCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGAcc                  >  1:410164/1‑151 (MQ=255)                                                                                                                                 gttgttCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGAcc                  >  1:453014/1‑151 (MQ=255)                                                                                                                                       cAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGACCTACGcc            >  1:275884/1‑151 (MQ=255)                                                                                                                                        aaTACCCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGACCTACGCCa           <  1:374273/151‑1 (MQ=255)                                                                                                                                        aaTACCCGCCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGACCTACGCCa           <  1:1102930/151‑1 (MQ=255)                                                                                                                                          tAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGACCTACGcc            >  1:52/1‑148 (MQ=255)                                                                                                                                               ggCACCCCCATTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGACCTACGCCAGAGTCaa    <  2:1042717/151‑1 (MQ=255)                                                                                                                                               ggCACCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGACCTACGCCAGAGTCa     <  2:459661/150‑1 (MQ=255)                                                                                                                                                 cATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGACCTACGCCAGAGTCAAcc  >  2:390328/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TCGGCTTCGGTCGCGCCGGTAAATTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAG‑TAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGACCTACGCCAGAGTCAACC  >  NC_000913/3575809‑3576102

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