Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,389,043	G→A	A28A (GCC→GCT)	aaeA ←	p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,389,043	0	G	A	100.0%	238.8 / NA	74	A28A (GCC→GCT)	aaeA	p‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/37);  total (37/37)

TTTTTCACCAGCTGGTTATCATGAACATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAGGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACA  >  NC_000913/3388915‑3389165                                                                                                                                 |                                                                                                                           tttttCACCAGCTGGTTATCATGAACATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTaaaaa                                                                                                                  >  2:994422/1‑139 (MQ=255) tttttCACCAGCTGGTTATCATGAACATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTaaaaa                                                                                                                  >  1:1611951/1‑139 (MQ=255)       aCCAGCTGGTTATCATGAACATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGa                                                                                                            >  2:1049014/1‑139 (MQ=255)         cAGCTGGTTATCATGAACATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATg                                                                                                          <  1:52053/139‑1 (MQ=255)         cAGCTGGTTATCATGAACATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATg                                                                                                          <  2:626103/139‑1 (MQ=255)              ggTTATCATGAACATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAgg                                                                                                      >  2:1448850/1‑138 (MQ=255)              ggTTATCATGAACATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAgg                                                                                                      <  1:1448850/138‑1 (MQ=255)              ggTTATCATGAACATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGc                                                                                                     <  1:1049014/139‑1 (MQ=255)                  aTCATGAACATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGa                                                                                                 <  1:1346755/139‑1 (MQ=255)                       gAACATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGAc                                                                                            <  1:1437073/139‑1 (MQ=255)                         aCATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGc                                                                                                     >  1:1578200/1‑128 (MQ=255)                         aCATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGc                                                                                                     <  2:1578200/128‑1 (MQ=255)                          cATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGc                                                                                                     <  2:1353383/127‑1 (MQ=255)                          cATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGc                                                                                                     >  1:1353383/1‑127 (MQ=255)                             tCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGAt                                                                                                           >  1:501376/1‑118 (MQ=255)                             tCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGAt                                                                                                           <  2:501376/118‑1 (MQ=255)                              cACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACg                                                                                     <  2:364821/139‑1 (MQ=255)                              cACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACg                                                                                     <  1:994422/139‑1 (MQ=255)                                    ggTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGTGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGAATAATACGAcctgg                                                                               >  1:1188010/1‑136 (MQ=255)                                    ggTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATg                                                                                                          <  1:995219/112‑1 (MQ=255)                                    ggTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATg                                                                                                          >  2:995219/1‑112 (MQ=255)                                    ggTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATg                                                                                                          >  2:1337560/1‑112 (MQ=255)                                    ggTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATg                                                                                                          <  1:1337560/112‑1 (MQ=255)                                     gTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGa                                                                              >  2:266378/1‑139 (MQ=255)                                     gTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGa                                                                              >  1:330483/1‑139 (MQ=255)                                         tGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGc                                                                          >  2:1694388/1‑139 (MQ=255)                                         tGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGc                                                                          >  2:638085/1‑139 (MQ=255)                                         tGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGc                                                                          >  2:1055311/1‑139 (MQ=255)                                              ccAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGaa                                                                                                <  2:1605115/112‑1 (MQ=255)                                              ccAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGaa                                                                                                >  1:1605115/1‑112 (MQ=255)                                              ccAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGAc                                                                                            <  1:1393467/116‑1 (MQ=255)                                              ccAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGAc                                                                                            >  2:1393467/1‑116 (MQ=255)                                              ccAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTAc                                                                     >  2:781756/1‑139 (MQ=255)                                              ccAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTAc                                                                     >  1:953879/1‑139 (MQ=255)                                              ccAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTAc                                                                     >  2:1472321/1‑139 (MQ=255)                                              ccAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTAc                                                                     >  1:1688957/1‑139 (MQ=255)                                               cAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACg                                                                    >  2:1751772/1‑139 (MQ=255)                                                      gTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAAt                                                             <  1:1472321/139‑1 (MQ=255)                                                       tCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAg                                                                                                       <  1:30012/96‑1 (MQ=255)                                                       tCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAg                                                                                                       >  2:30012/1‑96 (MQ=255)                                                       tCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGa                                                                                                 >  1:1181031/1‑102 (MQ=255)                                                       tCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGa                                                                                                 <  2:1181031/102‑1 (MQ=255)                                                       tCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAAtt                                                            <  1:434021/139‑1 (MQ=255)                                                       tCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAAtt                                                            >  2:1406387/1‑139 (MQ=255)                                                        ccGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAAttt                                                           <  2:1611951/139‑1 (MQ=255)                                                           gcgcGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCt                                                        <  1:781756/139‑1 (MQ=255)                                                                 tCGCAACGACGTCTGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAAttt                                                           >  1:1203945/1‑130 (MQ=255)                                                                 tCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAAttt                                                           <  2:1203945/130‑1 (MQ=255)                                                                            tCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACtt                                       <  2:953879/139‑1 (MQ=255)                                                                            tCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACtt                                       <  2:1688957/139‑1 (MQ=255)                                                                                        cgcgcgTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACg                           <  1:1694388/139‑1 (MQ=255)                                                                                          cgcgTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGaa                         >  1:177153/1‑139 (MQ=255)                                                                                          cgcgTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGaa                         <  1:1406387/139‑1 (MQ=255)                                                                                          cgcgTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGaa                         <  2:177153/139‑1 (MQ=255)                                                                                          cgcgTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGaa                         <  1:1751772/139‑1 (MQ=255)                                                                                           gcgTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACtt                                       >  2:1465770/1‑124 (MQ=255)                                                                                           gcgTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACtt                                       <  1:1465770/124‑1 (MQ=255)                                                                                            cgTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTg                                                                                                               >  1:1077723/1‑51 (MQ=255)                                                                                            cgTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTg                                                                                                               <  2:1077723/51‑1 (MQ=255)                                                                                            cgTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAAttt                                                           >  2:1612932/1‑103 (MQ=255)                                                                                            cgTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAAttt                                                           <  1:1612932/103‑1 (MQ=255)                                                                                                aCGCGTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGTAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTcg                   <  2:1188010/139‑1 (MQ=255)                                                                                                 cgcgTCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTcgc                  >  1:1708381/1‑139 (MQ=255)                                                                                                     tCCAGGGGGATTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGata              <  2:330483/139‑1 (MQ=255)                                                                                                              aTTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGt                                                                                >  2:48411/1‑64 (MQ=255)                                                                                                              aTTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGt                                                                                <  1:48411/64‑1 (MQ=255)                                                                                                              aTTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATaaa     >  2:1064620/1‑139 (MQ=255)                                                                                                              aTTCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATaaa     >  2:630331/1‑139 (MQ=255)                                                                                                               ttCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATaaaa    >  2:652529/1‑139 (MQ=255)                                                                                                                tCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGa                                                                                                         <  1:1361878/37‑1 (MQ=255)                                                                                                                tCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGa                                                                                                         >  2:1361878/1‑37 (MQ=255)                                                                                                                tCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACtt                                       >  2:1082247/1‑103 (MQ=255)                                                                                                                tCGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACtt                                       <  1:1082247/103‑1 (MQ=255)                                                                                                                 cGGTGTAATAGACCCAAGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACa  <  1:638085/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                           TTTTTCACCAGCTGGTTATCATGAACATTCACCTGGGTAATGAGTCCAGAAACGTCCGGCGCGATCGCAACGACGTCAGCGCTAAAGCGCGCGTCACGCGTCCAGGGGGATTCGGTGTAATAGACCCAGGCATTAAAAATTGCGATGAAGGCCAGAATGACTAATACGACCGTGATGGCCGTACGGGAGAATTTTCTTATTAGTGTTTTCACTTCAACCTCAAACGAACAGTCGCGATATCAAATAAAACA  >  NC_000913/3388915‑3389165

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