Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,381,789	G→A	pseudogene (13/126 nt)	ycjV →	pseudogene, putative ATP‑binding component of a transport system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,381,789	0	G	A	100.0%	234.0 / NA	72	pseudogene (13/126 nt)	ycjV	pseudogene, putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/38);  total (34/38)

CCATCGTGATGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCC  >  NC_000913/1381653‑1381924                                                                                                                                         |                                                                                                                                        ccATCGTGATGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATg                                                                                                                                       <  2:1478554/139‑1 (MQ=255)         aTGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGc                                                                                                                               >  2:570295/1‑139 (MQ=255)           gggAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTg                                                                                                                             >  1:1211987/1‑139 (MQ=255)                  cGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCg                                                                                                                        >  2:430180/1‑137 (MQ=255)                  cGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCg                                                                                                                        <  1:430180/137‑1 (MQ=255)                     ccGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTaaa                                                                                                                   <  1:34076/139‑1 (MQ=255)                        gAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                                     >  2:372473/1‑118 (MQ=255)                        gAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                                     <  1:372473/118‑1 (MQ=255)                              atTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCa                                                                                                          >  2:1343203/1‑139 (MQ=255)                                tCATCCGGACGCGCATGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTTCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATg                                                                                                        <  2:1077349/139‑1 (MQ=255)                                 cATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGc                                                                                                       <  2:1211987/139‑1 (MQ=255)                                   tCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAg                                                                                                     >  1:1323935/1‑139 (MQ=255)                                       gACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAGTTAACCGGTGCGGAAGTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGgaga                                                                                                 >  2:1024151/1‑139 (MQ=255)                                       gACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTaaa                                                                                                                   >  1:5302/1‑121 (MQ=255)                                       gACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTaaa                                                                                                                   <  2:5302/121‑1 (MQ=255)                                           cgcAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                                     >  2:1547149/1‑99 (MQ=255)                                           cgcAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                                     <  1:1547149/99‑1 (MQ=255)                                            gcAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                                     >  2:269358/1‑98 (MQ=255)                                            gcAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                                     <  1:269358/98‑1 (MQ=255)                                                  ggAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTa                                                                                      <  1:1113299/139‑1 (MQ=255)                                                   gAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTAt                                                                                     <  2:515547/139‑1 (MQ=255)                                                    agaaTAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTAGTGCGTTAAATGATTATCATg                                                                                                        >  1:965579/3‑119 (MQ=255)                                                    aaaaTAACGTTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATg                                                                                                        <  2:965579/119‑1 (MQ=255)                                                           cATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCAc                                                                                        <  2:1235145/128‑1 (MQ=255)                                                           cATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCAc                                                                                        >  1:1235145/1‑128 (MQ=255)                                                              ttCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGata                                                                          >  2:1170547/1‑139 (MQ=255)                                                               tCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGatat                                                                         >  1:170949/1‑139 (MQ=255)                                                                  gCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                                     >  1:237293/1‑76 (MQ=255)                                                                  gCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                                     <  2:237293/76‑1 (MQ=255)                                                                   ccAAAATTAGCGTGGCAGAATTAACCGGGGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAACTATCCCTATTCATTTTGATATGACg                                                                     >  2:1582783/1‑139 (MQ=255)                                                                            gCGTGGCAGAATTAACCGGTGCGGGATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCAt                                                            <  2:515941/139‑1 (MQ=255)                                                                                 gCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGa                                                                    >  2:1054966/1‑126 (MQ=255)                                                                                 gCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGa                                                                    <  1:1054966/126‑1 (MQ=255)                                                                                 gCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCtt                                                       >  2:117002/1‑139 (MQ=255)                                                                                       ttAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGc                                                 <  1:117002/139‑1 (MQ=255)                                                                                          aCCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTa                                                                                                             <  1:1564503/76‑1 (MQ=255)                                                                                          aCCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTa                                                                                                             >  2:1564503/1‑76 (MQ=255)                                                                                               tGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTcc                                                                                                                                   <  2:1283596/49‑1 (MQ=255)                                                                                               tGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTcc                                                                                                                                   >  1:1283596/1‑49 (MQ=255)                                                                                               tGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATg                                                                                                        <  2:1599965/76‑1 (MQ=255)                                                                                               tGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATg                                                                                                        >  1:1599965/1‑76 (MQ=255)                                                                                                 cGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCa                                                                                         <  1:1331728/89‑1 (MQ=255)                                                                                                 cGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCa                                                                                         >  2:1331728/1‑89 (MQ=255)                                                                                                  ggAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGaaa                                      <  1:570295/139‑1 (MQ=255)                                                                                                              ctACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCAt                                                            >  2:770023/1‑105 (MQ=255)                                                                                                              ctACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCAt                                                            <  1:770023/105‑1 (MQ=255)                                                                                                              ctACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCt                          <  2:29052/139‑1 (MQ=255)                                                                                                              ctACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCt                          <  1:1343203/139‑1 (MQ=255)                                                                                                                acacCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTa                                                                                                                     >  1:1423642/1‑46 (MQ=255)                                                                                                                acacCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTa                                                                                                                     <  2:1423642/46‑1 (MQ=255)                                                                                                                acacCACGGTTGGGGGCACGAGTTAATGGTCCGTGCCGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTaa                        >  2:99063/1‑139 (MQ=255)                                                                                                                 caccacGGTTGGGGGCACGAGTTAATGGTCCGTGCTgg                                                                                                                            >  2:774071/1‑38 (MQ=255)                                                                                                                 caccacGGTTGGGGGCACGAGTTAATGGTCCGTGCTgg                                                                                                                            <  1:774071/38‑1 (MQ=255)                                                                                                                    cacGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGAt                                                                                                               <  2:1284734/48‑1 (MQ=255)                                                                                                                    cacGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGAt                                                                                                               >  1:1284734/1‑48 (MQ=255)                                                                                                                      cGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGaaa                                                                                               >  1:731027/1‑62 (MQ=255)                                                                                                                      cGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGaaa                                                                                               <  2:731027/62‑1 (MQ=255)                                                                                                                         ttGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACg                                                                     <  2:865596/85‑1 (MQ=255)                                                                                                                         ttGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACg                                                                     >  1:865596/1‑85 (MQ=255)                                                                                                                           gggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAAt                                                                                                                  <  1:669470/38‑1 (MQ=255)                                                                                                                           gggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAAt                                                                                                                  >  2:669470/1‑38 (MQ=255)                                                                                                                           gggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGga                                                                                                   <  1:1231110/53‑1 (MQ=255)                                                                                                                           gggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGga                                                                                                   >  2:1231110/1‑53 (MQ=255)                                                                                                                           gggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGa             >  2:622797/1‑139 (MQ=255)                                                                                                                            ggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGaaaa                                                                                              >  2:60856/1‑57 (MQ=255)                                                                                                                            ggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGaaaa                                                                                              <  1:60856/57‑1 (MQ=255)                                                                                                                                cACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCa        <  2:850617/139‑1 (MQ=255)                                                                                                                                 aCGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCAt       <  2:170949/139‑1 (MQ=255)                                                                                                                                   gAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTc     >  1:1118754/1‑139 (MQ=255)                                                                                                                                      ttAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTcccc  <  2:1323935/139‑1 (MQ=255)                                                                                                                                      ttAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTcccc  <  2:858684/139‑1 (MQ=255)                                                                                                                                      ttAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTcccc  <  2:1118754/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        CCATCGTGATGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCC  >  NC_000913/1381653‑1381924

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