Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,181,413	A→C	E331D (GAA→GAC)	ygiC →	ATP‑Grasp family ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,181,413	0	A	C	98.7%	281.4 / ‑4.4	79	E331D (GAA→GAC)	ygiC	ATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (43/35);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTT  >  NC_000913/3181266‑3181558                                                                                                                                                    |                                                                                                                                                  ccTGCTGCCCGCTTATGTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAggg                                                                                                                                                <  2:225996/151‑1 (MQ=255)  ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACggg                                                                                                                                                >  2:388629/1‑150 (MQ=255)              tATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAg                                                                                                                                    <  2:531160/150‑1 (MQ=255)              tATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGc                                                                                                                                   >  1:483532/1‑151 (MQ=255)               aTTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGc                                                                                                                                   >  2:875218/1‑150 (MQ=255)               aTTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGc                                                                                                                                   <  1:947059/150‑1 (MQ=255)               aTTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCa                                                                                                                                  <  1:760598/151‑1 (MQ=255)                 tttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCaa                                                                                                                                 <  1:184610/150‑1 (MQ=255)                    gCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTcc                                                                                                                             <  2:76525/151‑1 (MQ=255)                        aaGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGGTCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTc                                                                                                                              <  1:418383/146‑1 (MQ=255)                        aaGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCAccc                                                                                                                         <  1:3409/151‑1 (MQ=255)                             gatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTAc                                                                                                                    >  1:461661/1‑151 (MQ=255)                               tcatcCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCg                                                                                                                  >  1:103765/1‑151 (MQ=255)                                catcCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGa                                                                                                                 <  2:1060842/151‑1 (MQ=255)                                  tcCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGaa                                                                                                                >  2:559144/1‑150 (MQ=255)                                  tcCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGaa                                                                                                                >  2:288367/1‑150 (MQ=255)                                   ccGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAAt                                                                                                              <  2:898275/151‑1 (MQ=255)                                    cGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAAtt                                                                                                             <  2:1039298/151‑1 (MQ=255)                                    cGCAAATGGAAAAATATGTGGTTAAAACGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGaaa                                                                                                               <  2:665418/149‑1 (MQ=255)                                      cAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCg                                                                                                           >  1:761049/1‑151 (MQ=255)                                       aaaTGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCgg                                                                                                          <  1:507370/151‑1 (MQ=255)                                         aTGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCCTTGAAGCAGCGGAAGGTCCGCATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCg                                                                                                        >  1:248339/1‑151 (MQ=255)                                         aTGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCg                                                                                                        <  2:752006/151‑1 (MQ=255)                                           ggAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGa                                                                                                       <  1:942015/150‑1 (MQ=255)                                              aaaaTATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGc                                                                                                   >  1:566087/1‑151 (MQ=255)                                               aaaTATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCt                                                                                                  >  1:807034/1‑151 (MQ=255)                                               aaaTATGGGGTTAAACCGATCTTCTCCCGGGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTACGGCGAAGACGGGATGATTGTTCAGCAATTCCACCTGTTACCGAAATTCGGCGagag                                                                                                    >  2:850262/1‑149 (MQ=255)                                                  tatGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCtat                                                                                                >  1:354090/1‑150 (MQ=255)                                                   atGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCtata                                                                                               <  2:583072/150‑1 (MQ=255)                                                     gtgGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATg                                                                                             <  2:1011926/150‑1 (MQ=255)                                                           aaaCCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTg                                                                                      >  2:459851/1‑151 (MQ=255)                                                            aaCCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCACTTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTgg                                                                                     >  1:630632/1‑151 (MQ=255)                                                            aaCCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTgg                                                                                     >  1:723738/1‑151 (MQ=255)                                                                 aTCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAgct                                                                                >  1:830922/1‑151 (MQ=255)                                                                  tcttctCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAgctg                                                                               >  1:303776/1‑151 (MQ=255)                                                                   cttctCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAgctgg                                                                              >  2:783835/1‑151 (MQ=255)                                                                   cttctCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAgctgg                                                                              >  2:602330/1‑151 (MQ=255)                                                                   cttctCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAgctgg                                                                              >  1:1015100/1‑151 (MQ=255)                                                                     tctcCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAgctg                                                                               <  2:536463/148‑1 (MQ=255)                                                                     tctcCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAgctg                                                                               <  1:555629/148‑1 (MQ=255)                                                                     tctcCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAg                                                                                  <  1:785220/145‑1 (MQ=255)                                                                        cccGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAgctggctg                                                                           <  2:1165084/149‑1 (MQ=255)                                                                        cccGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGctggctgg                                                                          >  1:1009137/1‑150 (MQ=255)                                                                           gTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGa                                                                       >  1:748788/1‑150 (MQ=255)                                                                             gAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACg                                                                    <  1:19034/151‑1 (MQ=255)                                                                                  cgcAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCa                                                                >  2:530426/1‑150 (MQ=255)                                                                                    cAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAAc                                                              >  1:657407/1‑150 (MQ=255)                                                                                     aaaCGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAAccc                                                            >  2:867236/1‑151 (MQ=255)                                                                                     aaaCGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAAccc                                                            <  1:115393/151‑1 (MQ=255)                                                                                      aaCGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACccg                                                           <  2:906791/151‑1 (MQ=255)                                                                                        cGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTAcc                                                                                                                   <  2:660838/93‑1 (MQ=255)                                                                                        cGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTAcc                                                                                                                   >  1:660870/1‑93 (MQ=255)                                                                                        cGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACccgcc                                                         >  1:229788/1‑151 (MQ=255)                                                                                        cGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACccgcc                                                         >  2:760844/1‑151 (MQ=255)                                                                                              gATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAAtt                                                   >  2:693228/1‑151 (MQ=255)                                                                                               aTCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTg                                                  <  1:836356/151‑1 (MQ=255)                                                                                               aTCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTg                                                  <  1:663326/151‑1 (MQ=255)                                                                                                 cATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCCCCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTgg                                                 <  1:122373/150‑1 (MQ=255)                                                                                                     gagaACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTc                                            >  1:767596/1‑151 (MQ=255)                                                                                                        aaCGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCTGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTg                                         >  2:545695/1‑151 (MQ=255)                                                                                                             cAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGAc                                    <  1:379248/151‑1 (MQ=255)                                                                                                                 aCCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTg                                <  1:795288/151‑1 (MQ=255)                                                                                                                 aCCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTg                                >  2:676027/1‑151 (MQ=255)                                                                                                                 aCCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTg                                <  2:286054/151‑1 (MQ=255)                                                                                                                   cATTGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGc                               >  2:482015/1‑150 (MQ=255)                                                                                                                     ttGAAGAAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCAt                             >  2:7931/1‑150 (MQ=255)                                                                                                                      tGAAGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCAtt                            <  1:48980/150‑1 (MQ=255)                                                                                                                        aaGCAGCGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGAt                         >  2:377661/1‑151 (MQ=255)                                                                                                                             gcGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCAccc                    >  2:1191588/1‑151 (MQ=255)                                                                                                                             gcGGAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCAccc                    >  2:395062/1‑151 (MQ=255)                                                                                                                               ggAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGACATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCAccc                    >  2:1172635/1‑149 (MQ=255)                                                                                                                                gAAGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAgg                 >  1:578663/1‑151 (MQ=255)                                                                                                                                  aGGTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGat               >  2:552763/1‑151 (MQ=255)                                                                                                                                   ggTCCGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCCGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGata              <  1:162597/151‑1 (MQ=255)                                                                                                                                       cGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGtc          <  2:139367/151‑1 (MQ=255)                                                                                                                                       cGTATGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGt           >  2:685050/1‑150 (MQ=255)                                                                                                                                         tATGGCGAAGACGGGATGCTTGTTCGGCAATTCCACCCGGTACCTAAATTCGGCGACAGCTATATCCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCGGGAATTGGCATTCGTGGAGACCATGCATTGATCACCCAGGATCTGtatc        >  2:76387/1‑151 (MQ=255)                                                                                                                                          aTGGCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCg       <  1:388646/151‑1 (MQ=255)                                                                                                                                            ggCGAAGACGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCgg      <  2:920401/150‑1 (MQ=255)                                                                                                                                               gaagacGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCTCCCAGGATATGTCTCGGtttt  <  1:90304/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTT  >  NC_000913/3181266‑3181558

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