Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,897,050	+G	coding (277/666 nt)	yieH →	phosphoenolpyruvate and 6‑phosphogluconate phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,897,046	1	.	G	95.7%	237.9 / 1.8	69	coding (273/666 nt)	yieH	phosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base G (27/39);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 7.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTTCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑G‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCA  >  NC_000913/3896908‑3897173                                                                                                                                            |                                                                                                                                   tttCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                     <  2:747016/139‑1 (MQ=255) tttCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                     >  1:905591/1‑139 (MQ=255) tttCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTGCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGGAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                     <  1:743559/139‑1 (MQ=255)      aaGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑g‑g‑ggg                                                                                                                             <  1:1044588/139‑2 (MQ=255)           gtAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGa                                                                                                                        <  2:395822/139‑1 (MQ=255)             aaaaCTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAgc                                                                                                                      <  1:428990/139‑1 (MQ=255)                aCTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCt                                                                                                                   <  2:44621/139‑1 (MQ=255)                  tGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTcc                                                                                                                 <  1:553280/139‑1 (MQ=255)                  tGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTcc                                                                                                                 <  2:1086125/139‑1 (MQ=255)                   gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCt                                                                                                                >  2:804279/1‑139 (MQ=255)                    tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAATCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTg                                                                                                               >  1:81604/1‑139 (MQ=255)                    tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGACGCCATAGAAGG‑G‑GG‑CTGGAGCGCTCCTg                                                                                                               >  2:283248/1‑139 (MQ=255)                    tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTg                                                                                                               >  2:175908/1‑139 (MQ=255)                    tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCACAACTGGAAGCCATCGAA‑G‑GAGG‑CTGGAGCGCTCCTg                                                                                                               >  2:722821/1‑139 (MQ=255)                            tATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGAt                                                                                                       <  1:804279/139‑1 (MQ=255)                            tATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGAt                                                                                                       <  2:35121/139‑1 (MQ=255)                              tCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCa                                                                                                     <  1:718216/139‑1 (MQ=255)                              tCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCa                                                                                                     <  1:217375/139‑1 (MQ=255)                                   atTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTg                                                                                                               >  2:296091/1‑124 (MQ=255)                                   atTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTg                                                                                                               <  1:296091/124‑1 (MQ=255)                                     tGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGcc                                                                                              >  2:696045/1‑139 (MQ=255)                                      gTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                             <  1:912143/139‑1 (MQ=255)                                        ttCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAt                                                                                           <  1:696045/139‑1 (MQ=255)                                         tCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtg                                                                                          <  1:283248/139‑1 (MQ=255)                                          cccTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgt                                                                                         <  2:1096473/139‑1 (MQ=255)                                             ttGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                      <  2:256279/139‑1 (MQ=255)                                              tGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtgt                                                                                     >  2:927532/1‑139 (MQ=255)                                              tGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAt                                                                                           >  2:603843/1‑133 (MQ=255)                                              tGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAt                                                                                           <  1:603843/133‑1 (MQ=255)                                               gAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtgtg                                                                                    >  1:851063/1‑139 (MQ=255)                                                 aCATGGTGTTACGTTAGCGAAATCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGt                                                                                  <  2:81604/139‑1 (MQ=255)                                                     ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                      >  2:25860/1‑131 (MQ=255)                                                     ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                      <  1:25860/131‑1 (MQ=255)                                                     ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCt                                                                              >  2:751109/1‑139 (MQ=255)                                                     ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCt                                                                              >  1:419947/1‑139 (MQ=255)                                                          tACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACgg                                                                         >  1:616563/1‑139 (MQ=255)                                                          tACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACgg                                                                         >  1:270241/1‑139 (MQ=255)                                                           aCGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGc                                                                        <  1:689622/139‑1 (MQ=255)                                                            cGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGcc                                                                       >  2:320583/1‑139 (MQ=255)                                                              ttAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCa                                                                     <  2:269573/139‑1 (MQ=255)                                                                aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCaaa                                                                   >  1:334428/1‑139 (MQ=255)                                                                 gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAAt                                                                  >  2:1076827/1‑139 (MQ=255)                                                                 gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAAt                                                                  >  2:408307/1‑139 (MQ=255)                                                                  cGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATa                                                                 <  1:1076827/139‑1 (MQ=255)                                                                      aaCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACaa                                                             >  2:508018/1‑139 (MQ=255)                                                                       aCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACaaa                                                            <  1:508018/139‑1 (MQ=255)                                                                              cTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagc                                                     <  2:419947/139‑1 (MQ=255)                                                                              cTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagc                                                     <  1:330420/139‑1 (MQ=255)                                                                               tGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagca                                                    <  1:656881/139‑1 (MQ=255)                                                                                        ttACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑GAG‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATggg                                           >  2:815935/1‑139 (MQ=255)                                                                                         tACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                          <  1:357514/139‑1 (MQ=255)                                                                                         tACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                          <  2:1032626/139‑1 (MQ=255)                                                                                         tACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                          <  1:175908/139‑1 (MQ=255)                                                                                                   aaGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAata                                <  2:951063/139‑1 (MQ=255)                                                                                                     gTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATg                              <  2:616563/139‑1 (MQ=255)                                                                                                        gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAgc                                                                                                                      <  1:607570/48‑1 (MQ=39)                                                                                                        gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAgc                                                                                                                      >  2:607570/1‑48 (MQ=255)                                                                                                         cTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGc                          >  1:694940/1‑139 (MQ=255)                                                                                                             gCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCacta                      >  1:103574/1‑139 (MQ=255)                                                                                                                  tCGATTCAGAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTcc                 <  1:320583/139‑1 (MQ=255)                                                                                                                          gAACTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATaaa         <  2:592443/139‑1 (MQ=255)                                                                                                                           aaCTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCg                                                                                                         <  1:653113/42‑1 (MQ=255)                                                                                                                           aaCTGGAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCg                                                                                                         >  2:653113/1‑42 (MQ=38)                                                                                                                           aaCTGGAAGCCATCGAA‑G‑G‑GGACTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAAc        >  2:256740/1‑139 (MQ=255)                                                                                                                               ggAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGtt    <  2:851063/139‑1 (MQ=255)                                                                                                                                gAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCg                                                                                                         <  2:982194/37‑1 (MQ=255)                                                                                                                                gAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCg                                                                                                         >  1:982194/1‑37 (MQ=255)                                                                                                                                gAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGc                                                                        >  2:1045148/1‑70 (MQ=255)                                                                                                                                gAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGc                                                                        <  1:1045148/70‑1 (MQ=255)                                                                                                                                gAAGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTc   >  2:741255/1‑139 (MQ=255)                                                                                                                                 aaGCCATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCa  <  2:103574/139‑1 (MQ=255)                                                                                                                                     cATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAata                                >  2:220537/1‑105 (MQ=255)                                                                                                                                     cATCGAAGG‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAata                                <  1:220537/105‑1 (MQ=255)                                                                                                                                            |                                                                                                                                   TTTCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑G‑G‑GG‑CTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCA  >  NC_000913/3896908‑3897173

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