Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,386,872	C→T	A833A (GCG→GCA)	mscM ←	mechanosensitive channel protein, miniconductance

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,386,872	0	C	T	100.0%	243.3 / NA	74	A833A (GCG→GCA)	mscM	mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/39);  total (35/39)

TTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGG  >  NC_000913/4386737‑4387001                                                                                                                                          |                                                                                                                                  tttGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAgcactgca                                                                                                                                <  1:695722/139‑1 (MQ=255)    ggtgatggtggtgatggCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGcactgcaccg                                                                                                                             >  2:51183/1‑139 (MQ=255)            tggtgatggCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAgg                                                                                                                          <  2:259842/134‑1 (MQ=255)            tggtgatggCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAgg                                                                                                                          >  1:259842/1‑134 (MQ=255)                  tggCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTc                                                                                                               <  1:1112020/139‑1 (MQ=255)                     cGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTcc                                                                                                              <  1:1481615/139‑1 (MQ=255)                     cGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAg                                                                                                            <  2:1461462/139‑1 (MQ=255)                         aCCCGTACCCGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACtt                                                                                                        >  2:1341482/1‑139 (MQ=255)                            cGTACCCGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                     <  2:1227143/139‑1 (MQ=255)                            cGTACCCGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                     <  2:463833/139‑1 (MQ=255)                            cGTACCCGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                     <  1:215440/139‑1 (MQ=255)                            cGTACCCGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGAAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCt                                                                                                     <  2:1417095/139‑1 (MQ=255)                                 ccGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGcc                                                                                                <  1:831179/139‑1 (MQ=255)                                   ggCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCt                                                                                              >  1:415560/1‑139 (MQ=255)                                      gTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTa                                                                                           <  2:635897/139‑1 (MQ=255)                                          aaTCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                                                                                       >  2:132043/1‑139 (MQ=255)                                             ccAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGa                                                                                    <  2:1506944/139‑1 (MQ=255)                                               aGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGcc                                                                                                <  1:99708/125‑1 (MQ=255)                                               aGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGcc                                                                                                >  2:99708/1‑125 (MQ=255)                                                    ctgc‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACa                                                                             <  1:132043/139‑1 (MQ=255)                                                      gc‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTgg                                                                                     <  1:504648/129‑1 (MQ=255)                                                      gc‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTgg                                                                                     >  2:504648/1‑129 (MQ=255)                                                       c‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTa                                                                                           <  1:648083/122‑1 (MQ=255)                                                       c‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTa                                                                                           >  2:648083/1‑122 (MQ=255)                                                                     ttCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCACCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAata                                                              >  1:901059/1‑139 (MQ=255)                                                                        cagcagCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATAttt                                                           >  1:851716/1‑139 (MQ=255)                                                                              cgcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAAt                                                                                                                     <  2:1465710/75‑1 (MQ=255)                                                                              cgcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAAt                                                                                                                     >  1:1465710/1‑75 (MQ=255)                                                                              cgcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGAc                                                                                                          >  2:445779/1‑86 (MQ=255)                                                                              cgcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGAc                                                                                                          <  1:445779/86‑1 (MQ=255)                                                                              cgcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAgg                                                     >  2:370381/1‑139 (MQ=255)                                                                               gcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGcactgcaccg                                                                                                                             >  1:573588/1‑66 (MQ=255)                                                                               gcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGcactgcaccg                                                                                                                             <  2:573588/66‑1 (MQ=255)                                                                               gcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTgg                                                                                     >  2:517673/1‑106 (MQ=255)                                                                               gcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTgg                                                                                     <  1:517673/106‑1 (MQ=255)                                                                               gcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGa                                                    >  2:819151/1‑139 (MQ=255)                                                                                cgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTTATTGGCTCCAGACTTTCTa                                                                                                    <  1:1103904/90‑1 (MQ=255)                                                                                cgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCt                                                                                                                >  2:735338/1‑78 (MQ=255)                                                                                cgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCt                                                                                                                <  1:735338/78‑1 (MQ=255)                                                                                cgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTGTa                                                                                                    >  2:1103904/1‑90 (MQ=255)                                                                                cgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGaa                                                   >  2:13645/1‑139 (MQ=255)                                                                                   gCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGcc                                                <  2:1547351/139‑1 (MQ=255)                                                                                             cACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGaa                                      >  1:110981/1‑139 (MQ=255)                                                                                                 aGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCACTGTACCGTGGGGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGa                                  <  2:1266589/139‑1 (MQ=255)                                                                                                      cGTGGTGATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCt                             <  2:1033106/139‑1 (MQ=255)                                                                                                         ggtgATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGAc                          <  2:1137432/139‑1 (MQ=255)                                                                                                         ggtgATGATAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGAc                          <  1:1391716/139‑1 (MQ=255)                                                                                                           tgatgaTAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACca                        >  2:313534/1‑139 (MQ=255)                                                                                                             atgatAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGAc                          >  2:812533/1‑135 (MQ=255)                                                                                                             atgatAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGAc                          <  1:812533/135‑1 (MQ=255)                                                                                                               gatAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAgg                                                                                  >  2:189/1‑77 (MQ=255)                                                                                                               gatAAACACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAgg                                                                                  <  1:189/77‑1 (MQ=255)                                                                                                                  aaaCACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATg                                    <  2:307455/120‑1 (MQ=255)                                                                                                                  aaaCACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATg                                    >  1:307455/1‑120 (MQ=255)                                                                                                                  aaaCACCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGcaca                 >  2:698695/1‑139 (MQ=255)                                                                                                                    acacCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGAcc                         >  1:421369/1‑129 (MQ=255)                                                                                                                    acacCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGAcc                         <  2:421369/129‑1 (MQ=255)                                                                                                                     cacCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACg                                                                                                  <  1:1053679/55‑1 (MQ=255)                                                                                                                     cacCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACg                                                                                                  >  2:1053679/1‑55 (MQ=255)                                                                                                                     cacCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCt                             <  1:587165/124‑1 (MQ=255)                                                                                                                     cacCAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCt                             >  2:587165/1‑124 (MQ=255)                                                                                                                     cacCAGAATGGCAATCAGCACTACACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATg              >  2:459825/1‑139 (MQ=255)                                                                                                                        cAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTg                                                                                             <  2:483329/57‑1 (MQ=255)                                                                                                                        cAGAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTg                                                                                             >  1:483329/1‑57 (MQ=255)                                                                                                                          gAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                                                                                       >  2:836102/1‑61 (MQ=255)                                                                                                                          gAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                                                                                       <  1:836102/61‑1 (MQ=255)                                                                                                                          gAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTg                                                                                <  1:197995/68‑1 (MQ=255)                                                                                                                          gAATGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTg                                                                                >  2:197995/1‑68 (MQ=255)                                                                                                                             tGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGAc                          >  1:805789/1‑119 (MQ=255)                                                                                                                             tGGCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGAc                          <  2:805789/119‑1 (MQ=255)                                                                                                                               gCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACg                                                                                                  <  1:1623230/45‑1 (MQ=255)                                                                                                                               gCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACg                                                                                                  >  2:1623230/1‑45 (MQ=255)                                                                                                                               gCAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTCACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCa    >  1:1309479/1‑139 (MQ=255)                                                                                                                                cAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCg                                                                                        >  1:1633722/1‑54 (MQ=255)                                                                                                                                cAATCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCg                                                                                        <  2:1633722/54‑1 (MQ=255)                                                                                                                                 aaTCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAgg  <  1:51183/139‑1 (MQ=255)                                                                                                                                  aTCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAgg  <  2:243373/138‑1 (MQ=255)                                                                                                                                  aTCAGCACTGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAgg  >  1:243373/1‑138 (MQ=255)                                                                                                                                          |                                                                                                                                  TTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGC‑‑AAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGG  >  NC_000913/4386737‑4387001

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