Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,094,234	T→C	T283A (ACA→GCA)	yggR ←	putative PilT family AAA+ ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,094,234	0	T	C	100.0%	268.6 / NA	74	T283A (ACA→GCA)	yggR	putative PilT family AAA+ ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/39);  total (35/39)

GTCACAAACGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGTGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACGGGGTCTTTTTCCTGCGCCGGAAATGAATCCAC  >  NC_000913/3094099‑3094369                                                                                                                                        |                                                                                                                                        gTCACAAACGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGtt                                                                                                                                      >  1:882332/1‑139 (MQ=255)      aaaCGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                                 <  1:811818/139‑1 (MQ=255)      aaaCGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                                 <  2:1798698/139‑1 (MQ=255)            ccTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAAtt                                                                                                                           >  2:529181/1‑139 (MQ=255)                cccccACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCaaa                                                                                                                       <  2:848466/139‑1 (MQ=255)                    cACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAgc                                                                                                                   <  1:493821/139‑1 (MQ=255)                     aCCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAgcg                                                                                                                  >  1:1489361/1‑139 (MQ=255)                          gTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGcc                                                                                                                >  2:665027/1‑136 (MQ=255)                          gTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGcc                                                                                                                <  1:665027/136‑1 (MQ=255)                           tGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcgc                                                                                                            <  2:882332/139‑1 (MQ=255)                           tGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcgc                                                                                                            >  2:1043776/1‑139 (MQ=255)                            gCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcgcg                                                                                                           <  1:1557866/139‑1 (MQ=255)                                aTAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTTGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtcct                                                                                                       >  2:100252/1‑139 (MQ=255)                                aTAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                                 <  1:498463/113‑1 (MQ=255)                                aTAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                                 >  2:498463/1‑113 (MQ=255)                                aTAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtcct                                                                                                       >  2:2032313/1‑139 (MQ=255)                                aTAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtcct                                                                                                       <  1:529181/139‑1 (MQ=255)                                aTAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtcc                                                                                                        >  2:248808/1‑138 (MQ=255)                                aTAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtcc                                                                                                        <  1:248808/138‑1 (MQ=255)                                 tAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtcctt                                                                                                      <  1:1629551/139‑1 (MQ=255)                                       tgctgAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGtt                                                                                                >  1:709112/1‑139 (MQ=255)                                        gctgAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGttt                                                                                               >  1:48996/1‑139 (MQ=255)                                           gAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATc                                                                                            <  2:48996/139‑1 (MQ=255)                                            aaaCGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                                 >  1:33241/1‑101 (MQ=255)                                            aaaCGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                                 <  2:33241/101‑1 (MQ=255)                                            aaaCGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATcc                                                                                           >  1:628660/1‑139 (MQ=255)                                            aaaCGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATcc                                                                                           <  1:1517857/139‑1 (MQ=255)                                               cGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACt                                                                                        >  2:987659/1‑139 (MQ=255)                                                   atcatcCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCa                                                                                    <  2:1489361/139‑1 (MQ=255)                                                   atcatcCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCa                                                                                    >  1:931320/1‑139 (MQ=255)                                                        ccccACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCttt                                                                               <  1:987659/139‑1 (MQ=255)                                                         cccACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCTCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCtttt                                                                              >  2:297287/1‑139 (MQ=255)                                                         cccACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcg                                                                                                             <  1:636283/108‑1 (MQ=255)                                                         cccACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcg                                                                                                             >  2:636283/1‑108 (MQ=255)                                                         cccACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCtttt                                                                              >  2:420519/1‑139 (MQ=255)                                                                gctgcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAg                                                                       <  2:628660/139‑1 (MQ=255)                                                                gctgcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAg                                                                       >  2:431170/1‑139 (MQ=255)                                                                  tgcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCTCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCa                                                                     <  1:297287/139‑1 (MQ=255)                                                                   gcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCAc                                                                    <  1:1043776/139‑1 (MQ=255)                                                                   gcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCAc                                                                    <  1:420519/139‑1 (MQ=255)                                                                                   aTGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtcctt                                                                                                      <  2:127747/89‑1 (MQ=255)                                                                                   aTGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtcctt                                                                                                      >  2:249578/1‑89 (MQ=255)                                                                                   aTGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtcctt                                                                                                      <  1:249578/89‑1 (MQ=255)                                                                                   aTGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtcctt                                                                                                      >  1:127747/1‑89 (MQ=255)                                                                                   aTGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCt                                                                                 <  1:1205242/110‑1 (MQ=255)                                                                                   aTGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCt                                                                                 >  2:1205242/1‑110 (MQ=255)                                                                                     gCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGc                                                  <  1:100252/139‑1 (MQ=255)                                                                                      cGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGcc                                                 <  1:2032313/139‑1 (MQ=255)                                                                                        gTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGt                                                                                                          <  1:984405/80‑1 (MQ=255)                                                                                        gTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGt                                                                                                          >  2:984405/1‑80 (MQ=255)                                                                                           aCTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCAc                                                                    <  2:1627592/115‑1 (MQ=255)                                                                                           aCTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCAc                                                                    >  1:1627592/1‑115 (MQ=255)                                                                                                 gggTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGt                                              <  1:418434/131‑1 (MQ=255)                                                                                                 gggTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGt                                              >  2:418434/1‑131 (MQ=255)                                                                                                    ttttCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCAc                                   <  2:931320/139‑1 (MQ=255)                                                                                                     tttCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCttt                                                                               >  1:691166/1‑94 (MQ=255)                                                                                                     tttCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCttt                                                                               <  2:691166/94‑1 (MQ=255)                                                                                                     tttCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCAc                                                                    <  2:1864500/105‑1 (MQ=255)                                                                                                     tttCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCAc                                                                    >  1:1864500/1‑105 (MQ=255)                                                                                                       tCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCa                                                                     <  1:2023187/102‑1 (MQ=255)                                                                                                       tCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCa                                                                     >  2:2023187/1‑102 (MQ=255)                                                                                                       tCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCAc                                                                    <  1:1772612/103‑1 (MQ=255)                                                                                                       tCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCAc                                                                    >  2:1772612/1‑103 (MQ=255)                                                                                                       tCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACggg                                >  1:1369466/1‑139 (MQ=255)                                                                                                              gcgAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcgcg                                                                                                           <  1:1166193/57‑1 (MQ=255)                                                                                                              gcgAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcgcg                                                                                                           >  2:1166193/1‑57 (MQ=255)                                                                                                                      aaTTCCCCACCGTGTGTGCGTTAATCAGCAATTCAAATAGCGGCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACGTGGTCTTTTTCCTgcgc                 <  1:766891/139‑1 (MQ=255)                                                                                                                            ccACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACgg                                 >  1:431933/1‑117 (MQ=255)                                                                                                                            ccACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACgg                                 <  2:431933/117‑1 (MQ=255)                                                                                                                             cACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACGGGGTCTTTTTCCTGCGCCGGAAAt          >  1:782192/1‑139 (MQ=255)                                                                                                                             cACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACGGGGTCTTTTTCCTGCGCCGGAAAt          >  1:1037743/1‑139 (MQ=255)                                                                                                                                cgcgGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTg                                                   <  1:1592938/95‑1 (MQ=255)                                                                                                                                cgcgGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTg                                                   >  2:1592938/1‑95 (MQ=255)                                                                                                                                cgcgGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACGGGGTCTTTTTCCTGCGCCGGAAATGaa       <  2:1037743/139‑1 (MQ=255)                                                                                                                                     gtgcgtTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACGGGGTCTTTTTCCTGCGCCGGAAATGAATccac  >  1:652966/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        GTCACAAACGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGTGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACGGGGTCTTTTTCCTGCGCCGGAAATGAATCCAC  >  NC_000913/3094099‑3094369

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