Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,227,040	A→G	A1071A (GCA→GCG)	metH →	cobalamin‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,227,040	0	A	G	100.0%	327.7 / NA	93	A1071A (GCA→GCG)	metH	cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (41/52);  total (41/52)

GCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGCC  >  NC_000913/4226904‑4227174                                                                                                                                         |                                                                                                                                       gCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTg                                                                                                                                     >  2:1625909/1‑140 (MQ=255) gCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTg                                                                                                                                     >  1:187350/1‑140 (MQ=255)  ccaccaTCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCGTTTGCCGAGACGGGAGGGCTGAAAGAGGACGCGCTgg                                                                                                                                    >  1:511113/1‑140 (MQ=255)  ccaccaTCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTgg                                                                                                                                    >  1:244685/1‑140 (MQ=255)      caTCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGa                                                                                                                                >  1:1846007/1‑140 (MQ=255)      caTCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGa                                                                                                                                <  2:1846007/140‑1 (MQ=255)              tCAACAGACCGAAAAAATAGGCATCGCTAACGACTCTCTTGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTg                                                                                                                        <  2:511113/140‑1 (MQ=255)              tCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTg                                                                                                                        <  2:1506012/140‑1 (MQ=255)                 aCAGACCGAAAAAACTGGCTTCGCTGACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAg                                                                                                                     <  1:882715/140‑1 (MQ=255)                 aCAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGCACGCGCTGGCTGATGCCTTTGAAg                                                                                                                     <  2:1261789/140‑1 (MQ=255)                   aGACCGAAAAAACAGGCCTCGCAAACTACTGCCTCTCTGTCTTCGTTGAGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAgcg                                                                                                                   <  2:756638/140‑1 (MQ=255)                    gACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAgcgc                                                                                                                  >  2:1616357/1‑140 (MQ=255)                    gACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGCCTGGAAGTGGATGCTCTGGCTGATGCCTTTGAAgcgc                                                                                                                  >  1:785273/1‑140 (MQ=255)                      ccGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCgcag                                                                                                                >  2:1667367/1‑140 (MQ=255)                      ccGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCgcag                                                                                                                >  1:1211719/1‑140 (MQ=255)                       cGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGa                                                                                                                                <  1:813821/123‑1 (MQ=255)                       cGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGa                                                                                                                                >  2:813821/1‑123 (MQ=255)                       cGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTg                                                                                                                        <  2:1189611/131‑1 (MQ=255)                       cGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTg                                                                                                                        >  1:1189611/1‑131 (MQ=255)                       cGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCgcagc                                                                                                               <  2:1702495/140‑1 (MQ=255)                          agagaCAGGCTGCGATAACTGCTGAATCACTGACTTCGTTTCGCGGAAGCTTTCTGCTTAAGCAGATTCCCTCGGAACCTTCGCCGTTACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCCCAGCACg                                                                                                            <  2:1260298/136‑1 (MQ=255)                          aaaaaCAGGCTTCGCTACCTACCGTCTCGCTAACTCCGTCGCGCCGAACCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACg                                                                                                            <  2:558046/140‑1 (MQ=255)                          aaaaaCAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGTCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACg                                                                                                            <  2:174572/140‑1 (MQ=255)                          aaaaaCAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACg                                                                                                            <  2:1211719/140‑1 (MQ=255)                          aaaaaCAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACg                                                                                                            <  1:169462/140‑1 (MQ=255)                          aaaaaCAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACg                                                                                                            <  2:935183/140‑1 (MQ=255)                                 ggCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACgatgat                                                                                                       >  2:1492881/1‑138 (MQ=255)                                 ggCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACgatgat                                                                                                       <  1:1492881/138‑1 (MQ=255)                                 ggCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTa                                                                                                     <  1:1667367/140‑1 (MQ=255)                                    ttCGCTAACTACTGTCTCGCTGCCTACGTCGCGCCGACGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTacaa                                                                                                  <  1:515781/140‑1 (MQ=255)                                    ttCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTacaa                                                                                                  <  2:39791/140‑1 (MQ=255)                                          aaCTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAAt                                                                                            <  1:1324892/140‑1 (MQ=255)                                           actactGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATc                                                                                           <  1:222020/140‑1 (MQ=255)                                           actactGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATc                                                                                           <  1:422057/140‑1 (MQ=255)                                            ctactGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACg                                                                                                            >  1:432296/1‑122 (MQ=255)                                            ctactGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACg                                                                                                            <  2:432296/122‑1 (MQ=255)                                             tactGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCAt                                                                                         <  1:1616357/140‑1 (MQ=255)                                                tGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGt                                                                                      >  1:1853505/1‑140 (MQ=255)                                                tGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAAAAAAATCATgga                                                                                      >  1:454960/1‑139 (MQ=255)                                                tGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCAGGCTGATGACTTGGAAGCGCAGCACGATGAATACAACAAAATCATGGt                                                                                      >  2:544845/1‑140 (MQ=255)                                                 gTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTg                                                                                     >  1:1906362/1‑140 (MQ=255)                                                  tctcGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCACAGCACGATGATTACAACAAAATCATGGTGa                                                                                    >  2:1188642/1‑140 (MQ=255)                                                    tcGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGaaa                                                                                  >  2:27115/1‑140 (MQ=255)                                                     cGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAg                                                                                 >  2:346873/1‑140 (MQ=255)                                                     cGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAg                                                                                 <  2:1618205/140‑1 (MQ=255)                                                     cGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAg                                                                                 >  2:1677314/1‑140 (MQ=255)                                                     cGCTGACTGCGTTGCGCGGAAGCTTTCTGGGGAAGCAGACTACATCGGTGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAg                                                                                 <  2:889177/140‑1 (MQ=255)                                                       cTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGAGCAGCACCATGATTACAACAAAATCATGGTGAAAgcg                                                                               >  1:756899/1‑140 (MQ=255)                                                            ttCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTCTGAAGCGCAGCAGGATGATTACAACAAAATGATGGTGAAAGCGCTTGc                                                                          >  2:641363/1‑140 (MQ=255)                                                             tCGTTGCGCCGAAGCTTTCTGGCAATGCAGATTACATGGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGcc                                                                         <  1:641363/140‑1 (MQ=255)                                                              cGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCTGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGa                                                                                    <  1:1423936/128‑1 (MQ=255)                                                              cGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCTGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGa                                                                                    >  2:1423936/1‑128 (MQ=255)                                                              cGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGGACGATGATTACAACAAAATCATGGTGAAAGCGCTTGccc                                                                        >  1:1831754/1‑139 (MQ=255)                                                              cGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGccc                                                                        >  2:1324200/1‑139 (MQ=255)                                                                    gcCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGAATACAACAAAATTATCGTGAAAGCGCTTGCCGACCGtt                                                                  >  1:931821/1‑140 (MQ=255)                                                                        aaGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATTATGGTGGAAGCGCCTGCCGACCGTTTAGc                                                              >  1:705237/1‑140 (MQ=255)                                                                         aGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGcc                                                             >  1:79998/1‑140 (MQ=255)                                                                           cTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGa                                                           <  2:1227461/140‑1 (MQ=255)                                                                                   tAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTg                                                   <  1:1079329/140‑1 (MQ=255)                                                                                   tAAAGCAGATTACATCGGCGCAGTTGCCGTGACTGGCGGTCTGGAAGATGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTg                                                   <  2:785273/140‑1 (MQ=255)                                                                                   tAAAACATATTTCAGCGGCGCATTGGCCCTGACTGGCGGGCTGTAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTg                                                   <  2:1831754/140‑1 (MQ=255)                                                                                     aaGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCg                                                 <  1:1677314/140‑1 (MQ=255)                                                                                       gCAGATTATATCGGCGACTTTGCCGTGACTGCCCGGCTGTAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGa                                               <  2:931821/140‑1 (MQ=255)                                                                                       gCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGa                                               <  1:1625909/140‑1 (MQ=255)                                                                                       caggATTACATCGGCGCATTTGGCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGa                                               <  1:653063/137‑1 (MQ=255)                                                                                         aGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGt                                             <  1:1016341/140‑1 (MQ=255)                                                                                         aGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGt                                             <  1:1324200/140‑1 (MQ=255)                                                                                         aGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGt                                             <  1:305396/140‑1 (MQ=255)                                                                                         aGATTACATCGGCGCATTTGCCGTGACTGGCCGGCAGGATGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGt                                             <  2:624563/140‑1 (MQ=255)                                                                                             tACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCg                                                 >  1:1341647/1‑132 (MQ=255)                                                                                             tACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCg                                                 <  2:1341647/132‑1 (MQ=255)                                                                                             tACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTAtct                                         <  2:79998/140‑1 (MQ=255)                                                                                                 tCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTCTAGCCGCAGCCTTGGCGGAGTATATCCAt                                     >  1:624563/1‑140 (MQ=255)                                                                                                  cGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGa                                               <  1:1774629/129‑1 (MQ=255)                                                                                                  cGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGa                                               >  2:1774629/1‑129 (MQ=255)                                                                                                     cgcATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCTTTTAGCCGAAGCATTTGAGGAGTATCTCCATgctg                                 >  2:653063/1‑137 (MQ=255)                                                                                                     cgcATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTg                                                   <  2:300614/122‑1 (MQ=255)                                                                                                     cgcATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGt                                             <  1:1019374/128‑1 (MQ=255)                                                                                                     cgcATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGt                                             >  2:1019374/1‑128 (MQ=255)                                                                                                     cgcATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGc                                 >  1:272233/1‑140 (MQ=255)                                                                                                     cgcATTTGCCGTGACTGGCGGGCTGGAAGAGGACGAGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTg                                                   >  1:300614/1‑122 (MQ=255)                                                                                                       cATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCgt                               <  1:1155966/140‑1 (MQ=255)                                                                                                       cATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCgt                               <  2:1853505/140‑1 (MQ=255)                                                                                                       cATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCgt                               >  1:1488466/1‑140 (MQ=255)                                                                                                       cATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCgt                               >  2:1000825/1‑140 (MQ=255)                                                                                                                 gACTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAg                     <  2:272233/140‑1 (MQ=255)                                                                                                                  aCTGGCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGt                    <  2:1906362/140‑1 (MQ=255)                                                                                                                      gCGGGCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTAc                <  1:27115/140‑1 (MQ=255)                                                                                                                      gCGGGCTGGAAGAGGACCCGCTTGCTGATGCCTTTCAAGAGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTAc                <  2:756899/140‑1 (MQ=255)                                                                                                                        gggCTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTg              <  2:244685/140‑1 (MQ=255)                                                                                                                           cTGGAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTgggg           >  2:60357/1‑140 (MQ=255)                                                                                                                             ggAAGAGGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCt         <  2:1488466/140‑1 (MQ=255)                                                                                                                                  agGACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATgcg    <  1:60357/140‑1 (MQ=255)                                                                                                                                    gACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGGCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGcc  >  2:940127/1‑140 (MQ=255)                                                                                                                                    gACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGcc  >  2:1823856/1‑140 (MQ=255)                                                                                                                                    gACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGcc  >  2:443225/1‑140 (MQ=255)                                                                                                                                    gACGCGCTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGcc  >  1:1111667/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                       GCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGCC  >  NC_000913/4226904‑4227174

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