Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,545,565	G→A	*228* (TAG→TAA)	gntX →	DNA catabolic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,545,565	0	G	A	100.0%	293.0 / NA	88	*228* (TAG→TAA)	gntX	DNA catabolic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (41/47);  total (41/47)

CCGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCC  >  NC_000913/3545432‑3545699                                                                                                                                      |                                                                                                                                       ccGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCt                                                                                                                                   <  1:731543/139‑1 (MQ=255)     gCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGTCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgat                                                                                                                               >  1:1841493/1‑139 (MQ=255)       aaGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                             <  1:1542746/139‑1 (MQ=255)          gTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGcc                                                                                                                                    <  2:668402/129‑1 (MQ=255)          gTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGcc                                                                                                                                    >  1:668402/1‑129 (MQ=255)            cGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGTCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                        <  2:1841493/139‑1 (MQ=255)            cGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                        <  2:1734010/139‑1 (MQ=255)            cGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                        <  2:241279/139‑1 (MQ=255)            cGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                        <  1:1247105/139‑1 (MQ=255)                     gTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTTCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATaa                                                                                                               >  1:354627/1‑139 (MQ=255)                     gTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATaa                                                                                                               <  2:773764/139‑1 (MQ=255)                         ttGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCaa                                                                                                           >  1:119152/1‑139 (MQ=255)                          tgtgGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAAc                                                                                                          <  2:1772606/139‑1 (MQ=255)                           gtgGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACt                                                                                                         >  2:623079/1‑139 (MQ=255)                            tgGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTa                                                                                                        <  1:706731/139‑1 (MQ=255)                               atgatgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatgat                                                                                                                            <  1:217895/116‑1 (MQ=255)                               atgatgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatgat                                                                                                                            >  2:217895/1‑116 (MQ=255)                               atgatgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTaaaa                                                                                                     <  1:1862057/139‑1 (MQ=255)                                   tgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCt                                                                                                                                   >  1:898619/1‑105 (MQ=255)                                   tgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCt                                                                                                                                   <  2:898619/105‑1 (MQ=255)                                   tgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGt                                                                                                 >  1:384094/1‑139 (MQ=255)                                       gTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAAc                                                                                             >  2:216137/1‑139 (MQ=255)                                        ttACGACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACt                                                                                            >  1:1091440/1‑139 (MQ=255)                                        ttACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                        <  2:36496/111‑1 (MQ=255)                                        ttACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                        >  1:36496/1‑111 (MQ=255)                                          accaccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                          <  2:119152/139‑1 (MQ=255)                                             accGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                        <  2:1695669/106‑1 (MQ=255)                                             accGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                        >  1:1695669/1‑106 (MQ=255)                                               cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                >  2:731838/1‑128 (MQ=255)                                               cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                <  1:731838/128‑1 (MQ=255)                                               cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                     <  2:1091440/139‑1 (MQ=255)                                               cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                     <  1:47860/139‑1 (MQ=255)                                               cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                     >  2:1772336/1‑139 (MQ=255)                                               cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                     <  1:1772336/139‑1 (MQ=255)                                                 gAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                   <  2:910484/139‑1 (MQ=255)                                                         gTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATg                                                                           <  2:383722/139‑1 (MQ=255)                                                            gCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATc                                                                        >  1:29476/1‑139 (MQ=255)                                                               gagaTTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                             >  2:993840/1‑83 (MQ=255)                                                               gagaTTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                             <  1:993840/83‑1 (MQ=255)                                                                    tGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                <  1:335858/107‑1 (MQ=255)                                                                    tGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                >  2:335858/1‑107 (MQ=255)                                                                      cgcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                        <  1:1744648/81‑1 (MQ=255)                                                                      cgcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                        >  2:1744648/1‑81 (MQ=255)                                                                            tGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                        <  1:188132/107‑1 (MQ=255)                                                                            tGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                        >  2:188132/1‑107 (MQ=255)                                                                                  tACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAg                                                  <  2:29476/139‑1 (MQ=255)                                                                                   aCGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAgc                                                 >  2:192986/1‑139 (MQ=255)                                                                                         tGGTGCTGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                     <  2:1785138/97‑1 (MQ=255)                                                                                         tGGTGCTGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                     >  1:1785138/1‑97 (MQ=255)                                                                                         tGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                     <  2:22913/97‑1 (MQ=255)                                                                                         tGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                     >  1:22913/1‑97 (MQ=255)                                                                                                      gTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                >  2:675561/1‑73 (MQ=255)                                                                                                      gTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                <  1:675561/73‑1 (MQ=255)                                                                                                        ccAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATaa                                                                                                               >  2:1107844/1‑56 (MQ=255)                                                                                                        ccAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATaa                                                                                                               <  1:1107844/56‑1 (MQ=255)                                                                                                             tCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                   >  2:185449/1‑79 (MQ=255)                                                                                                             tCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                   <  1:1012922/79‑1 (MQ=255)                                                                                                             tCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                   >  2:1012922/1‑79 (MQ=255)                                                                                                             tCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                   <  1:185449/79‑1 (MQ=255)                                                                                                                 gtgtCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                         <  2:5081/37‑1 (MQ=255)                                                                                                                 gtgtCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                         >  1:5081/1‑37 (MQ=255)                                                                                                                 gTGCCTTTGTCGTACCTTGTAAAGCCTCGATGATggg                                                                                                                         <  2:1145802/37‑1 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                         >  1:1145802/1‑37 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                <  2:13462/62‑1 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                >  1:13462/1‑62 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                <  2:197518/62‑1 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                >  1:197518/1‑62 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCa                                                                                               <  2:316358/63‑1 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCa                                                                                               >  2:999914/1‑63 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCa                                                                                               >  1:316358/1‑63 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCa                                                                                               <  1:999914/63‑1 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                        <  1:587851/70‑1 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                        >  2:587851/1‑70 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                        >  1:1342716/1‑70 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                        <  2:1342716/70‑1 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                   >  1:1831790/1‑75 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                   <  2:1831790/75‑1 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTc                                                                >  2:1136752/1‑94 (MQ=255)                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTc                                                                <  1:1136752/94‑1 (MQ=255)                                                                                                                    ccTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGgaaga                <  2:354627/139‑1 (MQ=255)                                                                                                                      tttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                   >  1:270852/1‑70 (MQ=255)                                                                                                                      tttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                   <  2:270852/70‑1 (MQ=255)                                                                                                                       ttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCATCTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGga                   <  2:1132582/133‑1 (MQ=255)                                                                                                                       ttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGga                   >  1:1132582/1‑133 (MQ=255)                                                                                                                                  tGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGgaaga                <  1:1466485/125‑1 (MQ=255)                                                                                                                                  tGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGgaaga                >  2:1466485/1‑125 (MQ=255)                                                                                                                                  tGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATcc  >  2:844937/1‑139 (MQ=255)                                                                                                                                  tGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATcc  >  2:1557614/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                       CCGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCC  >  NC_000913/3545432‑3545699

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