Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,858,929	A→C	T613P (ACC→CCC)	mutS →	methyl‑directed mismatch repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,858,929	0	A	C	100.0%	355.7 / NA	93	T613P (ACC→CCC)	mutS	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (38/55);  total (38/55)

TACCCTGAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCACCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGG  >  NC_000913/2858781‑2859065                                                                                                                                                     |                                                                                                                                         taCCCTGAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATcccc                                                                                                                                        >  1:280357/1‑151 (MQ=255)  aCCCTGAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCg                                                                                                                                       >  1:727805/1‑151 (MQ=255)       gAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCg                                                                                                                                  <  1:704142/151‑1 (MQ=255)            acacCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAAc                                                                                                                               <  1:265390/149‑1 (MQ=255)            acacCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACa                                                                                                                              <  1:561061/150‑1 (MQ=255)            acacCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACAt                                                                                                                             <  4:168168/151‑1 (MQ=255)            acacCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACAt                                                                                                                             <  3:59740/151‑1 (MQ=255)              acCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCGTCCCCGGTCCGAACATgg                                                                                                                           >  2:837773/1‑151 (MQ=255)               ccTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATgg                                                                                                                           <  1:494445/150‑1 (MQ=255)                 tGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCg                                                                                                                        >  2:340131/1‑151 (MQ=255)                  gCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGc                                                                                                                         <  1:192372/149‑1 (MQ=255)                  gCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCg                                                                                                                        <  1:646339/150‑1 (MQ=255)                   cccGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGt                                                                                                                      <  3:270217/151‑1 (MQ=255)                      gACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTaa                                                                                                                    >  1:1082023/1‑150 (MQ=255)                       aCCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTaaaa                                                                                                                  <  1:170664/151‑1 (MQ=255)                       aCCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTaaa                                                                                                                   >  2:1185369/1‑150 (MQ=255)                       aCCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTaaa                                                                                                                   >  4:457907/1‑150 (MQ=255)                         cTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGt                                                                                                                >  1:212397/1‑151 (MQ=255)                           tCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTa                                                                                                               <  1:482296/150‑1 (MQ=255)                            cATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTAcc                                                                                                             >  1:1208367/1‑151 (MQ=255)                            cATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTAc                                                                                                              <  4:95684/150‑1 (MQ=255)                             aTTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCt                                                                                                            >  2:641345/1‑151 (MQ=255)                               tGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCtat                                                                                                          >  2:401135/1‑151 (MQ=255)                                   aaaCCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATgc                                                                                                      <  3:313940/151‑1 (MQ=255)                                     aCCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATgcgc                                                                                                    <  4:247954/151‑1 (MQ=255)                                          gCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGAc                                                                                               <  2:1085377/151‑1 (MQ=255)                                          gCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGAc                                                                                               >  1:811090/1‑151 (MQ=255)                                            aTTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCg                                                                                             >  2:429968/1‑151 (MQ=255)                                              tCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGc                                                                                            >  4:447253/1‑150 (MQ=255)                                                gCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACt                                                                                         <  1:1167319/151‑1 (MQ=255)                                                gCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACt                                                                                         >  1:24172/1‑151 (MQ=255)                                                 cATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTg                                                                                        <  3:58256/151‑1 (MQ=255)                                                     aCCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTg                                                                                    >  2:107793/1‑151 (MQ=255)                                                      ccGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTgc                                                                                   <  4:34547/151‑1 (MQ=255)                                                       cGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTgcg                                                                                  >  1:837630/1‑151 (MQ=255)                                                       cGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTgcg                                                                                  >  4:425029/1‑151 (MQ=255)                                                          aGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTg                                                                               <  2:1125101/151‑1 (MQ=255)                                                               gCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGc                                                                          >  1:374448/1‑151 (MQ=255)                                                                 cATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCt                                                                        <  3:64886/151‑1 (MQ=255)                                                                     cGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACAt                                                                    >  4:185293/1‑151 (MQ=255)                                                                        gaGTCGAAGAAGTACTGAATGTGCAATTTATCGCCAACCCGCTGAATCTGTCGCGGCAGCGCCGCATGTTGATCATCCCCGGACCGAACATGGGCGGTAAAAGTACCTATATTCACAAGACCGCACTGATTGCGCTGATGGCCTACATCgg                                                                 >  3:281311/2‑151 (MQ=255)                                                                         aGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCgg                                                                 <  1:271528/150‑1 (MQ=255)                                                                         aGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGc                                                                <  1:31034/151‑1 (MQ=255)                                                                            tGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAg                                                              <  1:384472/150‑1 (MQ=255)                                                                            tGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAg                                                              <  1:427731/150‑1 (MQ=255)                                                                             gAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGc                                                             <  2:1235251/150‑1 (MQ=255)                                                                             gAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCt                                                            <  3:278391/151‑1 (MQ=255)                                                                             gAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCt                                                            <  2:35242/151‑1 (MQ=255)                                                                             gAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCt                                                            >  2:926665/1‑151 (MQ=255)                                                                               aCAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTAt                                                          <  2:988903/151‑1 (MQ=255)                                                                               aCAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTAt                                                          <  1:1031384/151‑1 (MQ=255)                                                                               aCAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTAt                                                          >  1:1251227/1‑151 (MQ=255)                                                                                 aaGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGt                                                        >  2:222120/1‑151 (MQ=255)                                                                                  aGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGt                                                        >  2:105708/1‑150 (MQ=255)                                                                                     aCTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCg                                                    >  2:390078/1‑151 (MQ=255)                                                                                       tGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGc                                                  >  2:502274/1‑151 (MQ=255)                                                                                        gAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGc                                                  <  2:185125/150‑1 (MQ=255)                                                                                        gAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGc                                                  >  1:61128/1‑150 (MQ=255)                                                                                        gAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGc                                                  <  1:1056072/150‑1 (MQ=255)                                                                                         aaTGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGcac                                                >  2:343286/1‑151 (MQ=255)                                                                                          aTGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGcaca                                               <  1:340189/151‑1 (MQ=255)                                                                                            gAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACaa                                              >  2:426047/1‑150 (MQ=255)                                                                                             aGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACaaaa                                            <  3:231238/151‑1 (MQ=255)                                                                                             aGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACaaaa                                            <  2:1259706/151‑1 (MQ=255)                                                                                                cATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGt                                         >  4:157679/1‑151 (MQ=255)                                                                                                   ttATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCg                                       <  2:1107907/150‑1 (MQ=255)                                                                                                    tATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCga                                      <  1:1068602/150‑1 (MQ=255)                                                                                                    tATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCg                                       >  3:103888/1‑149 (MQ=255)                                                                                                      tCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCg                                       >  1:466099/1‑147 (MQ=255)                                                                                                       cGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGAtt                                  <  1:1233905/151‑1 (MQ=255)                                                                                                        gCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTg                                 <  3:363901/151‑1 (MQ=255)                                                                                                         ccAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTgg                                >  4:96122/1‑151 (MQ=255)                                                                                                          cAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGa                               <  2:363628/151‑1 (MQ=255)                                                                                                           aaCCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGAc                              <  1:535383/151‑1 (MQ=255)                                                                                                           aaCCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGAc                              <  1:727783/151‑1 (MQ=255)                                                                                                             cccGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCt                            >  3:198148/1‑151 (MQ=255)                                                                                                             cccGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGAACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCt                            >  4:181934/1‑151 (MQ=255)                                                                                                              ccGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTa                           <  4:419788/151‑1 (MQ=255)                                                                                                              ccGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTa                           <  1:677015/151‑1 (MQ=255)                                                                                                               cGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCt                            >  3:135479/1‑149 (MQ=255)                                                                                                                 cTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTatcg                        >  4:378492/1‑151 (MQ=255)                                                                                                                   gAATCTGTCTGCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTCCCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACAACGGCAAATATGTTCCGGCACAAAAAGTCGAGGTTGGagcaaac                         <  2:218662/148‑7 (MQ=255)                                                                                                                   gAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTatcg                        <  4:375421/149‑1 (MQ=255)                                                                                                                    aaTCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTatcgatc                     <  2:1208151/151‑1 (MQ=255)                                                                                                                         gTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGc                   <  1:814764/148‑1 (MQ=255)                                                                                                                         gTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCa                  <  2:106050/149‑1 (MQ=255)                                                                                                                           cgccgcAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCtt              <  3:105351/151‑1 (MQ=255)                                                                                                                           cgccgcAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCt               <  1:1185585/150‑1 (MQ=255)                                                                                                                              cgcAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTAc           <  1:752181/151‑1 (MQ=255)                                                                                                                                cAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTAccc         >  3:19927/1‑151 (MQ=255)                                                                                                                                 aGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCcg        <  3:115654/151‑1 (MQ=255)                                                                                                                                    gccgcATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGt     <  1:890711/151‑1 (MQ=255)                                                                                                                                     ccgcATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTa    <  4:81074/151‑1 (MQ=255)                                                                                                                                       gcATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAgg  <  1:243004/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                         TACCCTGAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCACCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGG  >  NC_000913/2858781‑2859065

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