Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,246,015	T→C	A454A (GCA→GCG)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,246,015	0	T	C	98.9%	313.6 / ‑4.3	87	A454A (GCA→GCG)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (48/38);  minor base A (1/0);  total (49/38)
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.74e-01

GGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGG  >  NC_000913/1245871‑1246162                                                                                                                                                 |                                                                                                                                                    ggtgCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCa                                                                                                                                                <  2:581040/150‑1 (MQ=255)   tgCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGcc                                                                                                                                             >  1:1028706/1‑151 (MQ=255)    gCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCa                                                                                                                                            <  2:43480/151‑1 (MQ=255)         catcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCCCCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGGCGCGACCCACTGTAACGGCGCCCAGCCATTTgg                                                                                                                                       <  2:798528/151‑1 (MQ=255)         catcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTgg                                                                                                                                       >  2:1010278/1‑151 (MQ=255)          atcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTgg                                                                                                                                       >  2:1073061/1‑150 (MQ=255)          atcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGgc                                                                                                                                      >  2:1058025/1‑151 (MQ=255)                tttttCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATcc                                                                                                                                >  1:505704/1‑151 (MQ=255)                   ttCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCa                                                                                                                              <  1:284854/150‑1 (MQ=255)                   ttCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCAt                                                                                                                             >  2:1434041/1‑151 (MQ=255)                    tCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCAtt                                                                                                                            >  1:1389776/1‑151 (MQ=255)                    tCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCAtt                                                                                                                            >  1:345625/1‑151 (MQ=255)                         caGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGc                                                                                                                       <  1:769492/151‑1 (MQ=255)                          aGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGcc                                                                                                                      >  2:903540/1‑151 (MQ=255)                                 tctcCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACtttt                                                                                                               >  2:114663/1‑151 (MQ=255)                                 tctcCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACttt                                                                                                                <  1:299319/150‑1 (MQ=255)                                     ccGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCAcc                                                                                                           >  2:1176508/1‑151 (MQ=255)                                      cGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCg                                                                                                          <  1:854600/151‑1 (MQ=255)                                        gTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGAc                                                                                                        >  1:950798/1‑151 (MQ=255)                                          cATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACg                                                                                                       <  2:369178/150‑1 (MQ=255)                                           aTAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTGTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtg                                                                                                     <  2:883047/151‑1 (MQ=255)                                            tAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtgg                                                                                                    >  1:1086889/1‑151 (MQ=255)                                            tAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtgg                                                                                                    >  1:689874/1‑151 (MQ=255)                                            tAGGGGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtgg                                                                                                    >  2:753669/1‑151 (MQ=255)                                             aGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtggt                                                                                                   >  1:318970/1‑151 (MQ=255)                                              ggTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtggtg                                                                                                  <  2:121488/151‑1 (MQ=255)                                               gtgtGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGgtgt                                                                                                 >  1:1115882/1‑151 (MQ=255)                                               gtgtGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGgtgt                                                                                                 >  1:1112867/1‑151 (MQ=255)                                                 gtgCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGtt                                                                                                <  1:1387883/150‑1 (MQ=255)                                                      gAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGgc                                                                                           <  2:1059263/150‑1 (MQ=255)                                                        aCATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGgccgc                                                                                        >  1:179027/1‑151 (MQ=255)                                                                cAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgc                                                                                >  2:715172/1‑151 (MQ=255)                                                                 aGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgc                                                                                >  2:338627/1‑150 (MQ=255)                                                                  ggAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgca                                                                               <  2:179029/150‑1 (MQ=255)                                                                           aGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATgtgt                                                                     >  1:1470502/1‑151 (MQ=255)                                                                              tAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGtttt                                                                  >  2:424333/1‑151 (MQ=255)                                                                               aaTGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGtttt                                                                  <  1:412156/150‑1 (MQ=255)                                                                               aaTGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGtttt                                                                  <  2:579836/150‑1 (MQ=255)                                                                                aTGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGTCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGtttt                                                                  <  2:522028/149‑1 (MQ=255)                                                                                aTGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTtcg                                                                <  2:80369/151‑1 (MQ=255)                                                                                  gTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTtcgtc                                                              <  1:1016770/151‑1 (MQ=255)                                                                                      aTCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTCTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccg                                                          >  1:596045/1‑151 (MQ=255)                                                                                          ccACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccgtcg                                                       >  1:284370/1‑150 (MQ=255)                                                                                           cACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCgtcgccgtcgcc                                                     >  1:1005331/1‑151 (MQ=255)                                                                                            aCCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCgtcgccgtcgcc                                                     >  2:1307510/1‑150 (MQ=255)                                                                                             ccTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCAt                                                   <  1:860248/151‑1 (MQ=255)                                                                                                  ttttGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCAttt                                                 <  2:1028608/148‑1 (MQ=255)                                                                                                  ttttGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttgt                                              >  1:466872/1‑151 (MQ=255)                                                                                                   tttGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttgtt                                             >  2:1430106/1‑151 (MQ=255)                                                                                                   tttGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttgt                                              <  1:1287429/150‑1 (MQ=255)                                                                                                   tttGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttgt                                              <  2:1347684/150‑1 (MQ=255)                                                                                                      gCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTgg                                           <  1:457386/150‑1 (MQ=255)                                                                                                      gCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGgc                                          >  1:1241284/1‑151 (MQ=255)                                                                                                      gCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGgc                                          >  2:274893/1‑151 (MQ=255)                                                                                                       cccGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGgcg                                         <  1:555788/151‑1 (MQ=255)                                                                                                       cccGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGgcg                                         >  1:1088404/1‑151 (MQ=255)                                                                                                       cccGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGgca                                         <  2:290610/151‑2 (MQ=255)                                                                                                          gTAGTTTTGTAATCCTTCTGTCGCGACCCTCTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCGCCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGa                                      <  1:86102/151‑1 (MQ=255)                                                                                                          gTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGcgcg                                       >  2:837842/1‑150 (MQ=255)                                                                                                          gTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGAGCCAAGCCATCTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGcgcg                                       >  1:1052351/1‑150 (MQ=255)                                                                                                           tAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGAt                                     >  2:882467/1‑151 (MQ=255)                                                                                                            aGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGAt                                     <  2:1186555/150‑1 (MQ=255)                                                                                                            aGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATc                                    >  2:1224778/1‑151 (MQ=255)                                                                                                            aGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATc                                    <  1:808991/151‑1 (MQ=255)                                                                                                                ttGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCttt                                 <  1:1434146/150‑1 (MQ=255)                                                                                                                  gTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGc                              >  1:151349/1‑151 (MQ=255)                                                                                                                    aaTCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCt                             <  2:730684/150‑1 (MQ=255)                                                                                                                     aTCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGc                           >  1:1199288/1‑151 (MQ=255)                                                                                                                      tCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGcc                          >  2:123077/1‑151 (MQ=255)                                                                                                                        cTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGc                        >  1:933106/1‑151 (MQ=255)                                                                                                                          tCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCa                       <  1:110324/150‑1 (MQ=255)                                                                                                                          tCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCAt                      >  1:1132335/1‑151 (MQ=255)                                                                                                                             gTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTg                    <  1:1058121/150‑1 (MQ=255)                                                                                                                              tCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGAc                  >  1:1383562/1‑151 (MQ=255)                                                                                                                              tCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGAc                  >  1:420305/1‑151 (MQ=255)                                                                                                                                  gACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTAc              >  2:136430/1‑151 (MQ=255)                                                                                                                                    cccACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAg            <  1:594438/151‑1 (MQ=255)                                                                                                                                    cccACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAg            <  2:127566/151‑1 (MQ=255)                                                                                                                                     ccACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGc           <  1:1010373/151‑1 (MQ=255)                                                                                                                                       aCTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGGCGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGTCGCATTGACGTa               <  2:16490/145‑1 (MQ=255)                                                                                                                                       aCTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCgg         >  1:1361973/1‑151 (MQ=255)                                                                                                                                        cTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCggg        >  2:438359/1‑151 (MQ=255)                                                                                                                                          gTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCggg        <  2:103536/149‑1 (MQ=255)                                                                                                                                           tAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAAc     >  2:412141/1‑151 (MQ=255)                                                                                                                                            aaCGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAAc     >  2:627271/1‑150 (MQ=255)                                                                                                                                             aCGGCGCCCAGCCATTTGGCGCATCCCATTGGTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGaa      >  2:824615/1‑148 (MQ=255)                                                                                                                                              cGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAgg  <  2:1112768/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGG  >  NC_000913/1245871‑1246162

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