Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,438,060	C→T	R13C (CGC→TGC)	mscL →	mechanosensitive channel protein, high conductance

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,438,060	0	C→T	100.0%	142.9	40	R13C (CGC→TGC)	mscL	mechanosensitive channel protein, high conductance
Reads supporting (aligned to +/- strand):  new base (18/22):  ref base (0/0):  total (18/22)

GAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGCGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCAGGGGGATATCCCTGCTGTTGTGATGCATTACGGTGTCTTCATTCAAAAC  >  NC_000913/3437859‑3438266                                                                                                                                                                                                          |                                                                                                                                                                                                               gAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAAc                                                                                                                                                                                                        >  2:429051/1‑210 (MQ=255)        tctTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGtggtgg                                                                                                                                                                                                 <  2:126013/210‑1 (MQ=255)               gCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGAGGTGCGGGAACGTGGTGGATTTGGc                                                                                                                                                                                          >  1:417291/1‑210 (MQ=255)                 cAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCgg                                                                                                                                                                                        >  2:514264/1‑210 (MQ=255)                       cTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTTGATTTGGCggtgggtg                                                                                                                                                                                  >  2:308606/1‑210 (MQ=255)                       cTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCggtgggtg                                                                                                                                                                                  <  1:226494/210‑1 (MQ=255)                           cttcttTTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTTAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCAt                                                                                                                                                                              <  2:69764/210‑1 (MQ=255)                             tcttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTa                                                                                                                                                                            <  1:227196/210‑1 (MQ=255)                                  taattaatAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGt                                                                                                                                                                       >  1:53470/1‑210 (MQ=255)                                  gaaataataaGGCGTCATTAGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGt                                                                                                                                                                       >  1:544546/5‑210 (MQ=255)                                    attaatAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTTATTATCGGTGc                                                                                                                                                                     >  1:131142/1‑210 (MQ=255)                                                    aTGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTTCGGCATTCGGGAAGAtt                                                                                                                                                     >  1:206327/1‑210 (MQ=255)                                                               ttattTCCCTTTGATTATCAAATATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGAATTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACt                                                                                                                                          <  1:320051/210‑1 (MQ=255)                                                                  tttCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGt                                                                                                                                       <  2:68892/210‑1 (MQ=255)                                                                     cccTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGc                                                                                                                                    <  2:53471/210‑1 (MQ=255)                                                                      ccTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGcc                                                                                                                                   <  2:241916/210‑1 (MQ=255)                                                                          tGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGata                                                                                                                               <  1:2284/210‑1 (MQ=255)                                                                                tCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATAtcatca                                                                                                                         >  2:188383/1‑210 (MQ=255)                                                                                                       cTGGCAGGAAAATGGCTTAACATTTGTGAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAAtt                                                                                                  <  1:135324/210‑1 (MQ=255)                                                                                                         ggCAGGAAAATGGCTTAACATTTGTCAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCATTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTgg                                                                                                >  2:21940/1‑210 (MQ=255)                                                                                                          gCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGc                                                                                               <  2:76278/210‑1 (MQ=255)                                                                                                           cAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCg                                                                                              >  1:538465/1‑210 (MQ=255)                                                                                                               aaaaTGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATGGGCGGGAt                                                                                          >  2:19585/1‑210 (MQ=255)                                                                                                                         aaCATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTaaa                                                                                <  1:429028/210‑1 (MQ=255)                                                                                                                          aCATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAAc                                                                               <  2:243591/210‑1 (MQ=255)                                                                                                                          aCATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAAc                                                                               <  2:196755/210‑1 (MQ=255)                                                                                                                                   aGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTg                                                                      <  2:370503/210‑1 (MQ=255)                                                                                                                                                      tCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCa                                                   >  2:518282/1‑210 (MQ=255)                                                                                                                                                       cATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGTGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCAg                                                  >  1:507141/1‑210 (MQ=255)                                                                                                                                                           gggAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGTCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTTTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCAggggg                                              <  1:115696/210‑1 (MQ=255)                                                                                                                                                           gggAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCAggggg                                              >  2:155111/1‑210 (MQ=255)                                                                                                                                                            ggAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCAGGGGGa                                             >  2:497109/1‑210 (MQ=255)                                                                                                                                                                        gAGCATTATTAAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGGGATTCGCAGGGGGCTATCCCGGCtgt                                 >  1:12703/1‑210 (MQ=255)                                                                                                                                                                        gAGCATTATTAAAGAATGTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCAGGGATATCCCCctgcgca                                 <  2:53616/210‑4 (MQ=255)                                                                                                                                                                             ttattaAAGAATTTCGCGAATTTGGGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCAGGGGGATATCCCTGCTGTTGTGa                            <  1:270499/210‑1 (MQ=255)                                                                                                                                                                                 taAAGAATTTCGCGAATTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCAGGGGGATATCCCTGCTGTTGTGATGCa                        <  2:16095/210‑1 (MQ=255)                                                                                                                                                                                            gcgAGTTTGCGATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCAGGGGGATATCCCTGCTGTTGTGATGCATTACGGTGTCt             <  1:350590/210‑1 (MQ=255)                                                                                                                                                                                                      taTGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCAGGGGGATATCCCTGCTGTTGTGATGCATTACGGTGTCTTCATTCaaaa   <  1:160999/209‑1 (MQ=255)                                                                                                                                                                                                      gATGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCAGGGGGATATCCCTGCTGTTGTGATGCATTACGGTGTCTTCATTCaaaa   <  2:131149/210‑1 (MQ=255)                                                                                                                                                                                                       aTGTGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCAGGGGGATATCCCTGCTGTTGTGATGCATTACGGTGTCTTCATTCaaaaa  >  2:161007/1‑209 (MQ=255)                                                                                                                                                                                                          |                                                                                                                                                                                                               GAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGCGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTGGGCTTATTAATTGGCGGGATCGATTTTAAACAGTTTGCTGTCACGCTACGCGATGCGCAGGGGGATATCCCTGCTGTTGTGATGCATTACGGTGTCTTCATTCAAAAC  >  NC_000913/3437859‑3438266

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.