Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,953,212	G→A	A171T (GCC→ACC)	cmoA →	carboxy‑SAM synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,953,212	0	G	A	96.0%	259.8 / 1.6	76	A171T (GCC→ACC)	cmoA	carboxy‑SAM synthase
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (37/36);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

GCCATTGAAAACGCATCGATGGTGGTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCAT  >  NC_000913/1953065‑1953358                                                                                                                                                    |                                                                                                                                                   gCCATTGAAAACGCATCGATGGTGGTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATAcc                                                                                                                                                  >  1:920044/1‑150 (MQ=255)  ccATTGAAAACGCATCGATGGTGGTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCaa                                                                                                                                                <  2:984919/151‑1 (MQ=255)  ccATTGAAAACGCATCGATGGTGGTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCa                                                                                                                                                 >  1:643533/1‑150 (MQ=255)     ttGAAAACGCATCGATGGTGGTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGt                                                                                                                                             >  2:344369/1‑151 (MQ=255)       gAAAACGCATCGATGGTGGTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTg                                                                                                                                           >  2:893341/1‑151 (MQ=255)         aaaCGCATCGATGGTGGTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGt                                                                                                                                         <  2:733751/151‑1 (MQ=255)          aaCGCATCGATGGTGGTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTg                                                                                                                                        >  2:56632/1‑151 (MQ=255)          aaCGCATCGATGGTGGTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTg                                                                                                                                        <  2:463555/151‑1 (MQ=255)               aTCGATGGTGGTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACtg                                                                                                                                   >  1:834438/1‑151 (MQ=255)                tCGATGGTGGTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACtg                                                                                                                                   <  1:212695/150‑1 (MQ=255)                 cGATGGTGGTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACtgtt                                                                                                                                 <  2:749456/151‑1 (MQ=255)                    tggtggTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTgttgtt                                                                                                                              >  1:715515/1‑151 (MQ=255)                       tggtgCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCa                                                                                                                            >  1:129177/1‑150 (MQ=255)                         gtgCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACa                                                                                                                         <  1:378701/151‑1 (MQ=255)                                aTTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGcacca                                                                                                                  <  2:513297/151‑1 (MQ=255)                                   ttACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACg                                                                                                                <  1:333309/150‑1 (MQ=255)                                   ttACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACg                                                                                                                >  1:1047733/1‑150 (MQ=255)                                          gCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTaa                                                                                                         <  1:895405/150‑1 (MQ=255)                                          gCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTaa                                                                                                         <  1:1155171/150‑1 (MQ=255)                                           cAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAAc                                                                                                       <  2:771877/151‑1 (MQ=255)                                                ccTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGcc                                                                                                  >  1:195325/1‑151 (MQ=255)                                                      aCCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGt                                                                                            >  1:238921/1‑151 (MQ=255)                                                      aCCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGt                                                                                            >  2:471886/1‑151 (MQ=255)                                                        cTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGtt                                                                                           >  2:750045/1‑150 (MQ=255)                                                        cTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGtt                                                                                           >  2:923937/1‑150 (MQ=255)                                                        cTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGtt                                                                                           >  2:1159895/1‑150 (MQ=255)                                                          tCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTAc                                                                                         >  2:73049/1‑150 (MQ=255)                                                           ccGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTAc                                                                                         >  1:487152/1‑149 (MQ=255)                                                            cGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGc                                                                                      >  2:97609/1‑151 (MQ=255)                                                             gAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCg                                                                                     <  2:1143814/151‑1 (MQ=255)                                                                 gcCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACt                                                                                 >  2:1236355/1‑151 (MQ=255)                                                                 gcCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACt                                                                                 <  1:1108344/151‑1 (MQ=255)                                                                  ccAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTg                                                                                <  2:284277/151‑1 (MQ=255)                                                                    aGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTgg                                                                               >  2:523430/1‑150 (MQ=255)                                                                         ttACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATc                                                                          >  2:648758/1‑150 (MQ=255)                                                                          tACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATc                                                                          >  2:702634/1‑149 (MQ=255)                                                                             tGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGcc                                                                      >  2:253714/1‑150 (MQ=255)                                                                              ggATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGcc                                                                      <  1:693788/149‑1 (MQ=255)                                                                                aTAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCTCGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGaa                                                                  <  2:1113870/151‑1 (MQ=255)                                                                                  aaaaTTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAgc                                                                <  1:950696/151‑1 (MQ=255)                                                                                   aaaTTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAgcg                                                               >  1:154491/1‑151 (MQ=255)                                                                                     aTTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCgca                                                             <  2:311435/151‑1 (MQ=255)                                                                                     aTTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCgca                                                             <  2:1106391/151‑1 (MQ=255)                                                                                      tttATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAgcgc                                                              >  1:441185/1‑149 (MQ=255)                                                                                       ttATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCgcagc                                                           >  2:98781/1‑151 (MQ=255)                                                                                           cAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGc                                                       <  1:890805/151‑1 (MQ=255)                                                                                           cAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGc                                                       >  2:627031/1‑151 (MQ=255)                                                                                            aaGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCt                                                      <  1:130010/151‑1 (MQ=255)                                                                                             aGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAGGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTg                                                     >  2:1023181/1‑151 (MQ=255)                                                                                             aGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTg                                                     <  1:256090/151‑1 (MQ=255)                                                                                                  tGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGaaa                                                 >  2:297586/1‑150 (MQ=255)                                                                                                   gAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGaaaa                                                <  1:98778/150‑1 (MQ=255)                                                                                                        ccGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGAt                                          <  1:173705/151‑1 (MQ=255)                                                                                                        ccGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGAt                                          <  2:589985/151‑1 (MQ=255)                                                                                                         cggcggCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGAt                                          >  1:509794/1‑150 (MQ=255)                                                                                                            cggcgCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTg                                      <  2:1137741/151‑1 (MQ=255)                                                                                                               cgcgcTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGAcc                                   <  1:712626/151‑1 (MQ=255)                                                                                                                 cgcTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCg                                  >  2:153641/1‑150 (MQ=255)                                                                                                                  gcTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGAt                                >  1:797782/1‑151 (MQ=255)                                                                                                                     ggTGCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTcc                             >  2:595091/1‑151 (MQ=255)                                                                                                                        gCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGtt                          <  1:597562/151‑2 (MQ=255)                                                                                                                        gCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGt                           <  1:717384/150‑1 (MQ=255)                                                                                                                        gCTTTCGGAAAAATTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTg                          <  1:88550/151‑1 (MQ=255)                                                                                                                             cGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAAc                     <  1:947631/151‑1 (MQ=255)                                                                                                                                 aaaaTTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCAt                 <  1:794752/151‑1 (MQ=255)                                                                                                                                  aaaTTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCAt                 >  2:1224009/1‑150 (MQ=255)                                                                                                                                   aaTTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATaa               <  1:1165784/151‑1 (MQ=255)                                                                                                                                   aaTTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATa                >  1:582073/1‑150 (MQ=255)                                                                                                                                    aTTCAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATaa               >  1:889280/1‑150 (MQ=255)                                                                                                                                       cAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGc            <  1:704223/150‑1 (MQ=255)                                                                                                                                       cAGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCa           <  1:572071/151‑1 (MQ=255)                                                                                                                                        aGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCAc          <  1:1046053/151‑1 (MQ=255)                                                                                                                                        aGTTTCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCAc          <  1:771305/151‑1 (MQ=255)                                                                                                                                            tCGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCt      >  1:288035/1‑151 (MQ=255)                                                                                                                                             cGAAGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTg     <  1:567478/151‑1 (MQ=255)                                                                                                                                                aGATACCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCAt  <  2:195328/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCCATTGAAAACGCATCGATGGTGGTGCTGAATTTTACCCTGCAATTCCTGGAACCTTCCGAGCGCCAGGCGTTACTGGATAAAATTTATCAAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCAT  >  NC_000913/1953065‑1953358

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