Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,396,097	A→G	R8G (AGA→GGA)	insH‑4 →	IS5 transposase and trans‑activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,396,097	0	A	G	100.0%	220.9 / NA	68	R8G (AGA→GGA)	insH‑4	IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (40/28);  total (40/28)

CGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAG  >  NC_000913/1395960‑1396177                                                                                                                                          |                                                                                 cGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGa                                                                                <  2:1760072/140‑1 (MQ=255)  ggAAACCCCTGGCGGGGCTATTTTAGCGCCACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGaa                                                                               <  1:1102919/140‑1 (MQ=255)  ggAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGaa                                                                               <  2:288523/140‑1 (MQ=255)  ggAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGaa                                                                               <  1:93267/140‑1 (MQ=255)  ggAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGaa                                                                               <  1:2466216/140‑1 (MQ=255)   gAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAAc                                                                              <  2:1790341/140‑1 (MQ=255)       ccccTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAggg                                                                          <  2:619777/140‑1 (MQ=255)       ccccTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAggg                                                                          <  2:1272595/140‑1 (MQ=255)       ccccTGGCGGGGCTACTTTACCTCGTCAATTCATTCAAATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAggg                                                                          <  1:1169537/140‑1 (MQ=25)           tGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTtca                                                                      >  1:1299181/1‑140 (MQ=255)             gCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCTTGATATGAGATCGTGTTTGTCATCTTGAGCAATGGAACAGGGGtcatc                                                                    >  2:876653/1‑140 (MQ=37)               ggggCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTcatcat                                                                  <  2:271723/140‑1 (MQ=255)               ggggCTATTTTACCGAGACAATCCACTCAGATCATCACTAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTcatcat                                                                  <  1:787605/140‑1 (MQ=25)                  gCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAg                                                               >  1:1576880/1‑140 (MQ=255)                  gCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAg                                                               >  2:2191662/1‑140 (MQ=255)                      ttttACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGtcat                                                           <  2:345453/140‑1 (MQ=255)                           ccGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACt                                                      >  1:2427342/1‑140 (MQ=255)                            cgcgACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACtt                                                     >  1:1372590/1‑140 (MQ=255)                            cgcgACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACtt                                                     >  1:498697/1‑140 (MQ=255)                             gcgACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATAAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTa                                                    >  1:1398434/1‑140 (MQ=255)                                aCAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCt                                                 >  2:1663256/1‑140 (MQ=255)                                aCAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCt                                                 >  1:1994780/1‑140 (MQ=255)                                aCAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCt                                                 >  2:1358743/1‑140 (MQ=255)                                    ttcattcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGc                                             <  2:2355971/140‑1 (MQ=255)                                    ttcattcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGc                                             <  2:2017552/140‑1 (MQ=255)                                      cattcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCg                                           <  1:2191662/140‑1 (MQ=255)                                         tcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGTGCCATGGAACAGCGTTGATCATGAGACATCAACTTACCTTCGCCGACa                                        >  2:990500/1‑140 (MQ=18)                                         tcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGACCAGAGTTCATCATGAGTCATCACCTTACCTTCGCCGACa                                        >  2:1060671/1‑140 (MQ=255)                                         tcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACa                                        <  2:2427342/140‑1 (MQ=255)                                           aGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGt                                      <  1:1337616/140‑1 (MQ=255)                                           aGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGt                                      >  1:1235584/1‑140 (MQ=255)                                              tcatcaATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGaa                                   <  2:1372590/140‑1 (MQ=255)                                                atcaATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGGCATGAAGTTACCTGCGCCGACAGTGAAtt                                 >  1:649433/1‑140 (MQ=18)                                                 tcaATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTc                                >  1:436540/1‑140 (MQ=255)                                                 tcaATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTAATCGTGAATCATCAACGTACCGTCGCCGACAGTGAATTc                                >  2:1126734/1‑140 (MQ=17)                                                      aGTCAGGGAAGGAAGTAGCAACATTGGCTGAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTcagcag                           <  1:876653/140‑1 (MQ=255)                                                      aGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTcagcag                           <  1:69374/140‑1 (MQ=255)                                                       gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGTAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGt                          >  1:1832095/1‑140 (MQ=255)                                                       gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAATTTAACTGCGCCGAGAGTGAATTCAGCAGt                          >  1:647797/1‑140 (MQ=17)                                                       gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGt                          >  1:2090355/1‑140 (MQ=255)                                                       gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGt                          >  1:1241903/1‑140 (MQ=255)                                                       gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGt                          >  1:1372510/1‑140 (MQ=255)                                                       gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGt                          >  2:1974719/1‑140 (MQ=255)                                                       gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGAAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGt                          >  1:1690144/1‑140 (MQ=255)                                                         cAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGaa                                   <  1:154219/129‑1 (MQ=255)                                                         cAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGaa                                   >  2:154219/1‑129 (MQ=255)                                                          aGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCGTCAACTTACCTACGCCGACAGCGAATTCAGCAGTAAg                       >  1:1045887/1‑140 (MQ=255)                                                          aGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAg                       >  1:613913/1‑140 (MQ=255)                                                           gggAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTAATCCAATTACCTTCGCGGACATAGAATTCAGCAGTAAgc                      >  2:1209256/1‑140 (MQ=14)                                                            ggaaggaaGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCCTCAACTTACCTTCGCCTACAGCGAATTCAGCAGTAAgcg                     >  2:783410/1‑140 (MQ=255)                                                            ggaaggaaGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAgcg                     >  2:1360515/1‑140 (MQ=255)                                                               aggaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAAt                                  <  1:1662547/124‑1 (MQ=255)                                                               aggaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAAt                                  >  2:1662547/1‑124 (MQ=255)                                                                ggaagTAGCAACATTAGCTAAGCAACGTGTGAACAAGGCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGt                 <  1:1209256/140‑1 (MQ=14)                                                                 gaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTc                >  1:330765/1‑140 (MQ=255)                                                                 gaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCACCAGTAAGCGCCGTc                >  1:125077/1‑140 (MQ=255)                                                                  aagTAGCAACATTAGTTAAGGACGGCGCGCACAAGTCGCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCa               <  1:783410/140‑1 (MQ=14)                                                                  aagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTGCCTTCGCAGACTGTGAATTCATCAGTAAGCGCCGTCa               >  1:1169745/1‑140 (MQ=14)                                                                    gTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGa             <  2:1690144/140‑1 (MQ=255)                                                                     tAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGAc            <  2:431459/140‑1 (MQ=255)                                                                      aGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGAcc           >  1:1885760/1‑140 (MQ=38)                                                                      aGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGAcc           >  1:2081954/1‑140 (MQ=38)                                                                      aGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCAGCAACTTACGTTCGCCGGCAGTGAATCCAGCGGTAAGCGCCGTCAGAcc           >  2:775144/1‑140 (MQ=11)                                                                       gCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCa          >  2:2350384/1‑140 (MQ=37)                                                                       gCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCa          <  1:1663256/140‑1 (MQ=37)                                                                          aCATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGaa       <  2:1372510/140‑1 (MQ=21)                                                                           cATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGaaa      >  2:2413689/1‑140 (MQ=18)                                                                               aGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATGGAACAGGGTTCATCATGAGTCATCAACTTACTTTCGCCGACAGTTAATTCAGCAGTAGGCGCAGTCAGACCAGAAAagag  >  2:1048210/1‑140 (MQ=11)                                                                                                                                          |                                                                                 CGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAG  >  NC_000913/1395960‑1396177

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