Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,082,192	T→C	Y94H (TAT→CAT)	loiP →	metalloprotease LoiP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,082,192	0	T	C	100.0%	166.8 / NA	47	Y94H (TAT→CAT)	loiP	metalloprotease LoiP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/25);  total (22/25)

TGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCACGTGGCGTTAGG  >  NC_000913/3082055‑3082325                                                                                                                                          |                                                                                                                                      tGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGcat                                                                                                                                     >  2:519107/1‑140 (MQ=255) tGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGcat                                                                                                                                     >  1:669915/1‑140 (MQ=255)  gAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGcata                                                                                                                                    >  2:172109/1‑140 (MQ=255)       aTCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCg                                                                                                                               <  2:641778/140‑1 (MQ=255)       aTCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCg                                                                                                                               <  2:445762/140‑1 (MQ=255)                       gaTGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCt                                                                                                               <  2:109300/140‑1 (MQ=255)                        aTGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCtt                                                                                                              <  2:466746/140‑1 (MQ=255)                          ggACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGAtgtg                                                                                                                      <  2:94498/130‑1 (MQ=255)                          ggACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGAtgtg                                                                                                                      >  1:94498/1‑130 (MQ=255)                           gACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGc                                                                                                           >  1:665032/1‑140 (MQ=255)                             cAGCAAGGCGACGATTGTGCCACCCAATAGCGAATTCGCGAAACGTCTTACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCaa                                                                                                         <  2:199232/140‑1 (MQ=255)                               gCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATg                                                                                                       <  1:630250/140‑1 (MQ=255)                                cAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATgg                                                                                                      >  1:63702/1‑140 (MQ=255)                                              cgcCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATc                                                                                        >  2:670648/1‑140 (MQ=255)                                               gccagccaATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATcc                                                                                       <  2:63702/140‑1 (MQ=255)                                                   gccaATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATcc                                                                                       <  1:515816/136‑1 (MQ=255)                                                   gccaATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATcc                                                                                       >  2:515816/1‑136 (MQ=255)                                                   gccaATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCAGCGt                                                                                   <  2:80379/140‑1 (MQ=255)                                                             aaTACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCggg                                                                         <  2:520932/140‑1 (MQ=255)                                                             aaTACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCggg                                                                         >  1:520932/1‑140 (MQ=255)                                                               tACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCt                                                                       <  2:665032/140‑1 (MQ=255)                                                                aCGCTAAACGTCTGACATCTATTGCCAATGCTCTATGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTg                                                                      <  2:273096/140‑1 (MQ=255)                                                                aCGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTg                                                                      <  1:438570/140‑1 (MQ=255)                                                                aCGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTg                                                                      <  1:670648/140‑1 (MQ=255)                                                                    tAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAg                                                                             <  1:655963/129‑1 (MQ=255)                                                                    tAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAg                                                                             >  2:655963/1‑129 (MQ=255)                                                                       aCGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGata                                                               <  1:365962/140‑1 (MQ=255)                                                                           cTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATatgat                                                           <  2:635743/140‑1 (MQ=255)                                                                            tGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATatgatg                                                          <  1:519107/140‑1 (MQ=255)                                                                            tGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATatgatg                                                          >  1:397180/1‑140 (MQ=255)                                                                            tGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATatgatg                                                          >  2:72616/1‑140 (MQ=255)                                                                               cAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACg                                                       <  2:397180/140‑1 (MQ=255)                                                                                            tGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGTTGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCg                                          >  2:611192/1‑140 (MQ=255)                                                                                                        caaTATCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCg                              >  2:615033/1‑140 (MQ=255)                                                                                                             tCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCAc                         >  2:649944/1‑140 (MQ=255)                                                                                                             tCAACGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATAGGATGACGAATAACGGAGTCGAGGCGGTGATCGGTCAc                         >  1:267273/1‑140 (MQ=255)                                                                                                               aaCGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGa                       >  2:604549/1‑140 (MQ=255)                                                                                                               aaCGGTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATACGATCAAGGATAACGAACTCGAACCGGTGATCGGTCACGa                       >  2:317886/1‑140 (MQ=255)                                                                                                                 cGGTCAGCCGGTAAACTACAAACTGTATATGGCGAAGAATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGaaa                     <  1:345089/140‑1 (MQ=255)                                                                                                                  ggTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGACGGATATGACGACAGATAACGAAGTCGAAGCAGTGATCGGCCACGAAAt                    >  2:345089/1‑140 (MQ=255)                                                                                                                   gTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGATGGATAACGAAGGCGGAGCGGTGATCGGTCACGAAATg                   >  1:259457/1‑140 (MQ=255)                                                                                                                   gTCAGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATg                   >  1:383793/1‑140 (MQ=255)                                                                                                                      aGCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATgggg                >  2:574312/1‑140 (MQ=255)                                                                                                                       gCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGc               >  1:406453/1‑140 (MQ=255)                                                                                                                       gCCGGTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGAGAACGAAGTCGAGGCTGTGATCAGTCACGAAATGGGGc               >  1:226362/1‑140 (MQ=255)                                                                                                                          ggTAAATTACAAAGTGCTTATGCCGAAGGATTTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCACg            <  2:259457/140‑1 (MQ=255)                                                                                                                          ggTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCACg            <  1:649944/140‑1 (MQ=255)                                                                                                                          ggTAAATTACAAAGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCCATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCACg            <  2:413878/140‑1 (MQ=255)                                                                                                                                    aaaGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCACGTGGCGTTAgg  <  1:172109/140‑1 (MQ=255)                                                                                                                                      aGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATgggg                <  1:139522/124‑1 (MQ=255)                                                                                                                                      aGTGCATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATgggg                >  2:139522/1‑124 (MQ=255)                                                                                                                                          |                                                                                                                                      TGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCACGTGGCGTTAGG  >  NC_000913/3082055‑3082325

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