Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,747	0	T	G	45.3%	100.8 / 20.4	64	T255P (ACC→CCC)	citG	2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (4/31);  new base G (29/0);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 2.98e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAC  >  NC_000913/646599‑646895                                                                                                                                                     |                                                                                                                                                      aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGtt                                                                                                                                                      <  1:909679/150‑1 (MQ=255)  atatTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCg                                                                                                                                                    >  1:176666/1‑151 (MQ=255)      ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAt                                                                                                                                                 <  2:56213/150‑1 (MQ=255)       tCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAt                                                                                                                                                 <  1:551005/149‑1 (MQ=255)       tCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATg                                                                                                                                                <  2:207143/150‑1 (MQ=255)         cAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATG‑ccc                                                                                                                                            >  1:513270/1‑151 (MQ=255)          aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑cccc                                                                                                                                           >  1:435226/1‑151 (MQ=255)                     gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑CCCCCcttttt                                                                                                                                    >  2:27196/1‑146 (MQ=255)                       gCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑CCCCCCTTTTTCaata                                                                                                                               >  1:191122/1‑150 (MQ=255)                          tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATg                                                                                                                            <  2:608698/150‑1 (MQ=255)                           aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGtt                                                                                                                          <  2:603173/151‑1 (MQ=255)                            aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAATG‑CCCCCCTTTTTTaataat                                                                                                                             >  1:500937/1‑147 (MQ=255)                             ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCCTTTTTCAATAATGtttt                                                                                                                        >  2:1145584/1‑150 (MQ=255)                               aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTgcg                                                                                                                      <  1:33739/151‑1 (MQ=255)                               aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCCTTTTTCAATAATGttt                                                                                                                         >  2:894461/1‑148 (MQ=255)                                   aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCATG‑GCCCCCCTTTTCAATAATTTTTTCGCCTc                                                                                                                  >  1:1188954/1‑151 (MQ=255)                                   aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGAGGGGTGGAAAA‑TCCCCGTTTTTTAATAATGttttggcc                                                                                                                    >  2:176793/1‑144 (MQ=255)                                    aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATT‑CCCCCTTTTTTTTATAATTTTTTCGccccc                                                                                                                 >  2:257293/1‑148 (MQ=255)                                     aTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGCTGTGACCCCCTTTTTGAAATAAGGtttcgccgccc                                                                                                                 >  1:831289/1‑143 (MQ=255)                                         gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTcgcg                                                                                                               <  1:784208/148‑1 (MQ=255)                                         gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTcgc                                                                                                                <  2:340641/147‑1 (MQ=255)                                         gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCt                                                                                                             <  1:120180/150‑1 (MQ=255)                                         gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCCTTTTTTAATAATATTTTTGCCCCgcgc                                                                                                              >  1:287890/1‑149 (MQ=255)                                             tAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCCTTTTTCCATAAAGTTTTTGCCCCCCGCTGTAg                                                                                                         >  1:938307/1‑150 (MQ=255)                                              aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAg                                                                                                         <  2:99567/149‑1 (MQ=255)                                              aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAg                                                                                                         <  1:106654/149‑1 (MQ=255)                                                 aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCTTTTTTCAAAAATGtttttgccccccgccg                                                                                                            >  2:259124/1‑130 (MQ=255)                                                       gTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGc                                                                                              <  1:762585/151‑1 (MQ=255)                                                         aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAG‑CCCCCCTTTTTCAATAATTTTTTCGCCCCCCCCTGTAGCCCGCGCAGGccc                                                                                            >  2:1161481/1‑151 (MQ=255)                                                          aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATG‑CCCCCCTTTTTCAAAAATGTTTGCGCCCCGCCCTGTGGGCAGGGCACGcccc                                                                                           >  1:52363/1‑151 (MQ=255)                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                       <  2:563174/150‑1 (MQ=255)                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAG‑CCCCCTTTTTTCAATAATGTTTTTGCCCCGCGCCGTAGGCAGCGCGGGgccccccc                                                                                       >  1:77955/1‑148 (MQ=255)                                                                       cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcg                                                                              <  1:1203971/151‑1 (MQ=255)                                                                       cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcg                                                                              <  1:296839/151‑1 (MQ=255)                                                                        aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑GCCCCCTTTTTTAATAATCTTTGCCACCCCCCCTGAAGCCCGCGCAAGcccccccccccccccc                                                                               >  2:374328/1‑139 (MQ=255)                                                                           aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCCTTTTTCAATAAAGTTTTCGCCCCCCGCCGTAGCCCGCGCCGGgccccccccccccc                                                                                 >  1:696626/1‑136 (MQ=255)                                                                              gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGc                                                                        <  1:938582/150‑1 (MQ=255)                                                                                   cTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATG‑CCCCCCTTTTTTAATAATGTTTTCGCCTCGCGCCGGAGCCAGCGCAGGCCCCCCCCCCCACCCCCCGATGCAACgct                                                                  >  1:166105/1‑149 (MQ=255)                                                                                   cTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATG‑CCCCCCTTTTTTAAAAAAATTTTGGCCTCGCGCTGTAGCCAGCGCAGGGCCCCCCCCCCCCCCCGCCAAGCGACggt                                                                  >  1:212533/1‑149 (MQ=255)                                                                                   cTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAT‑CCCCCCTTTTTTAATAATGTTTTTGCCCCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCAcccc                                                                               >  2:570183/1‑136 (MQ=255)                                                                                      gACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGt                                                               <  1:28030/151‑1 (MQ=255)                                                                                        cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑CCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCCGTAGCCAGCGCAGGCCCCCCCCGCCCCCCCCCGCTGCAAccgtgg                                                                >  1:38655/1‑143 (MQ=255)                                                                                          gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTAt                                                             <  1:771123/149‑1 (MQ=255)                                                                                            gaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGcc                                                         <  2:347186/151‑1 (MQ=255)                                                                                                 cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTg                                                     <  1:516303/150‑1 (MQ=255)                                                                                                     cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCg                                                 <  1:329165/150‑1 (MQ=255)                                                                                                      gATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTTGCATG‑CCCCCTTTTTTTAATAATTTTTTTGCCTCGCCCTGTTGCCAGCGCAGGGCCCCCTCGCCCCCCCGCGACGCAACGGTGGTATCGCCGTTGATCGcc                                               >  2:749386/1‑151 (MQ=255)                                                                                                         acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAt                                             <  2:287892/150‑1 (MQ=255)                                                                                                          cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATG‑CCCCCCTTTTGCAATAATGTTTTCGCCCCCCGCTGGAGCCACCGCAGGCCCCCCCCGCCCCCCCGCGGTGCAACGGTGGGATTTCCGTTGATCGCCacc                                            >  2:357960/1‑148 (MQ=255)                                                                                                             tCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCTTTTTTCCATAATGTTTTCGCCCCCCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGATGCAACGTTGGGATCGCCGTTGATCGCCATCAGt                                         >  2:1114871/1‑150 (MQ=255)                                                                                                               ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAgg                                      <  2:212534/151‑1 (MQ=255)                                                                                                               ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAgg                                      <  1:83167/151‑1 (MQ=255)                                                                                                                  gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTTGAATG‑CCCCCTTTTTTGCATAATAATTTCTCCCCCCCCCGTAGCCAGCGCCGGACCCCCCCCCCCCCCCGCGGGGCAAAGGTGGGAATGCCGGTGATCGCCAACCGTAgg                                      >  2:84921/1‑148 (MQ=255)                                                                                                                  gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAg                                    <  2:251241/150‑1 (MQ=255)                                                                                                                       aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGt                              <  2:700845/151‑1 (MQ=255)                                                                                                                        cTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATG‑CCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCCATGCAACGTTGGTATCGCCGGTGATCGCCCTCAGTAGGAg                                    >  2:241019/1‑144 (MQ=255)                                                                                                                           ccGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCTTTTTTGAATAATGTTTTCGCCCCCCGCCGGAGCCCGCGCCGGCCCCCCCCCCCCCCCCGCCGCGCAACCTTGGTATCGCCGTTGATCGCCCTCCGTAGGAGCaaaggat                            >  1:110883/1‑144 (MQ=255)                                                                                                                              gagaTAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                       <  2:856957/151‑1 (MQ=255)                                                                                                                                   aaTCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATg                   <  1:1119677/150‑1 (MQ=255)                                                                                                                                    aTCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGcc                 <  1:760487/151‑1 (MQ=255)                                                                                                                                           tCGGCGGGGGGTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGGATCGAGCAAAGCCAGTTc            >  2:925842/1‑149 (MQ=255)                                                                                                                                             ggcggGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGa        <  1:360152/151‑1 (MQ=255)                                                                                                                                                gggggTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt     >  1:268907/1‑151 (MQ=255)                                                                                                                                                gggggTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt     >  1:1175040/1‑151 (MQ=255)                                                                                                                                                  gggTTCGAATG‑‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAc  <  1:1100667/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAC  >  NC_000913/646599‑646895

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