Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,340,588	T→C	L105S (TTA→TCA)	yrbG →	putative calcium/sodium:proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,340,588	0	T	C	100.0%	325.4 / NA	94	L105S (TTA→TCA)	yrbG	putative calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (47/47);  total (47/47)

CGTCTCTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATTACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTCAGGC  >  NC_000913/3340453‑3340723                                                                                                                                        |                                                                                                                                        cgTCTCTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCAcc                                                                                                                                      >  2:148489/1‑139 (MQ=255)   tctcTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCt                                                                                                                                    <  2:109999/139‑1 (MQ=255)      ctGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTaa                                                                                                                                 <  2:312282/139‑1 (MQ=255)      ctGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTaa                                                                                                                                 <  2:1059810/139‑1 (MQ=255)           cGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAAtgttg                                                                                                                            >  2:1858850/1‑139 (MQ=255)                 ccgcgACTTAGCCGTTGGTACAGCCCTCGGTTCAAACATTATCAATATATTGCTGATCCTCGGTCTTGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCATCCTTAATGTTGTTGGTc                                                                                                                      >  1:1810386/2‑139 (MQ=255)                 acGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTc                                                                                                                      >  2:476749/1‑139 (MQ=255)                       tttAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCgtg                                                                                                                <  2:646314/139‑1 (MQ=255)                          aGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCCTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCgtggtg                                                                                                             >  2:1562497/1‑139 (MQ=255)                                    aCAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                                   <  2:1792513/139‑1 (MQ=255)                                             gacTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATg                                                                                          <  2:1716214/137‑1 (MQ=255)                                                tCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACg                                                                                       <  2:465113/139‑1 (MQ=255)                                                 cAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACgg                                                                                      <  1:476749/139‑1 (MQ=255)                                                  aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAAtgttg                                                                                                                            >  2:1792433/1‑100 (MQ=255)                                                  aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAAtgttg                                                                                                                            <  1:1792433/100‑1 (MQ=255)                                                  aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGa                                                                                     >  1:1159478/1‑139 (MQ=255)                                                  aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGa                                                                                     >  1:649648/1‑139 (MQ=255)                                                   aaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGAc                                                                                    <  2:1810386/139‑1 (MQ=255)                                                                 atTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGc                                                                      >  1:324449/1‑139 (MQ=255)                                                                   tGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTGAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGa                                                                    <  1:245977/139‑1 (MQ=255)                                                                   tGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGa                                                                    >  1:567528/1‑139 (MQ=255)                                                                   tGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGa                                                                    <  2:1159478/139‑1 (MQ=255)                                                                    gCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                                   >  2:214812/1‑107 (MQ=255)                                                                    gCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                                   <  1:214812/107‑1 (MQ=255)                                                                     cTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCg                                                                          <  2:595149/131‑1 (MQ=255)                                                                     cTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCg                                                                          >  1:595149/1‑131 (MQ=255)                                                                     cTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAg                                                                       <  1:1880902/134‑1 (MQ=255)                                                                     cTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAg                                                                       >  2:1880902/1‑134 (MQ=255)                                                                      tGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATgg                                                                 <  2:324449/139‑1 (MQ=255)                                                                        aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCATCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                               >  1:267275/1‑139 (MQ=255)                                                                        aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                               >  1:1842950/1‑139 (MQ=255)                                                                        aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                               >  1:1217575/1‑139 (MQ=255)                                                                        aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                               >  1:471913/1‑139 (MQ=255)                                                                        aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                               >  2:698636/1‑139 (MQ=255)                                                                          ccTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGc                                                                                                                   >  1:977062/1‑85 (MQ=255)                                                                          ccTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGc                                                                                                                   <  2:977062/85‑1 (MQ=255)                                                                          ccTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTcc                                                                                                    <  1:1396860/100‑1 (MQ=255)                                                                          ccTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTcc                                                                                                    >  2:1396860/1‑100 (MQ=255)                                                                          ccTCGGTCTGGCTGCGCTGGTTCGCCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                               >  1:441500/1‑137 (MQ=255)                                                                          ccTCGGTCTGGCTGCGCTGGTTCGCCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                               <  2:441500/137‑1 (MQ=255)                                                                                tCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTtctc                                                       <  1:148489/139‑1 (MQ=255)                                                                                tCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTtctc                                                       <  2:1842950/139‑1 (MQ=255)                                                                                tCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTtctc                                                       <  1:1858850/139‑1 (MQ=255)                                                                                   ggCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTTGCAGCGATGGTATCTTTctcctc                                                    >  2:1600728/1‑139 (MQ=255)                                                                                    gttGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGAATCttt                                                          >  1:1126552/3‑132 (MQ=255)                                                                                    gCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGtggtggt                                                                                                           <  1:221566/83‑2 (MQ=255)                                                                                    gCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGtggtggt                                                                                                           >  2:221566/1‑82 (MQ=255)                                                                                    gCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGAATCttt                                                          <  2:1126552/132‑1 (MQ=255)                                                                                          cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGc                                                                         <  1:1785823/111‑1 (MQ=255)                                                                                          cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGc                                                                         >  2:1785823/1‑111 (MQ=255)                                                                                          cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                             >  1:335008/1‑139 (MQ=255)                                                                                          cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                             >  2:1376961/1‑139 (MQ=255)                                                                                          cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                             <  1:776280/139‑1 (MQ=255)                                                                                           tGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGc                                            <  1:1376961/139‑1 (MQ=255)                                                                                              ttCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCAcc                                                                                                                                      <  1:901148/46‑1 (MQ=255)                                                                                              ttCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCAcc                                                                                                                                      >  2:901148/1‑46 (MQ=255)                                                                                              ttCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCtgt                                         >  2:735954/1‑139 (MQ=255)                                                                                                       ttACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGa                                                                    <  1:914396/103‑1 (MQ=255)                                                                                                       ttACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGa                                                                    >  2:914396/1‑103 (MQ=255)                                                                                                       ttACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGgct                                <  2:471913/139‑1 (MQ=255)                                                                                                         aCCGTCCATTCTGATGTTCTAAGCCGTGAATCACCCTTAGTGTTGTTAGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCAGGCTGTGCTATGgctgc                              <  1:1600728/139‑1 (MQ=255)                                                                                                             tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCgtg                                                                                                                <  2:557555/53‑1 (MQ=255)                                                                                                             tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCgtg                                                                                                                >  1:557555/1‑53 (MQ=255)                                                                                                             tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGAc                                                                                    <  2:1561373/81‑1 (MQ=255)                                                                                                             tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGAc                                                                                    >  1:1561373/1‑81 (MQ=255)                                                                                                             tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAg                                                                       <  1:1104226/94‑1 (MQ=255)                                                                                                             tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAg                                                                       >  2:1104226/1‑94 (MQ=255)                                                                                                             tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGtt                          >  1:1539610/1‑139 (MQ=255)                                                                                                             tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGtt                          <  1:226016/139‑1 (MQ=255)                                                                                                               cATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCa                        >  1:1653861/1‑139 (MQ=255)                                                                                                                 ttgtGATGTTCGACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCg                                                                                                                  >  1:1090850/4‑47 (MQ=38)                                                                                                                 ttCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCg                                                                                                                  <  2:1090850/47‑1 (MQ=255)                                                                                                                        gTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                           <  2:1761633/47‑1 (MQ=255)                                                                                                                        gTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                           >  1:1761633/1‑47 (MQ=255)                                                                                                                           cTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTGGTCGCAGCGa                                                                    <  2:1311883/83‑1 (MQ=255)                                                                                                                           cTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTGGTCGCAGCGa                                                                    >  1:1311883/1‑83 (MQ=255)                                                                                                                           cTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTg            <  2:267275/139‑1 (MQ=255)                                                                                                                              cGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCAc         >  2:822556/1‑139 (MQ=255)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                           <  1:1751015/36‑1 (MQ=38)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                           >  2:1751015/1‑36 (MQ=38)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                         <  1:1827730/38‑1 (MQ=255)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                         >  2:1827730/1‑38 (MQ=255)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGt                                                                                                       >  2:1284213/1‑40 (MQ=39)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGt                                                                                                       <  1:1284213/40‑1 (MQ=255)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGa                                                                                         <  2:406019/54‑1 (MQ=255)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGa                                                                                         >  1:406019/1‑54 (MQ=255)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACgg                                                                                      >  2:467898/1‑57 (MQ=255)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCttt                                                 >  2:17820/1‑94 (MQ=255)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCttt                                                 <  1:17820/94‑1 (MQ=255)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                             <  1:1431632/98‑1 (MQ=255)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                             >  2:1431632/1‑98 (MQ=255)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTCAg    >  1:573484/1‑139 (MQ=255)                                                                                                                                   tGAATCACCCTTAATGTTGTTGGTAGGCGTGGTGGCCGGTTCCGTACTCTATGACgg                                                                                      <  1:467898/57‑1 (MQ=255)                                                                                                                                     aaTCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTCAGGc  <  2:1539610/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        CGTCTCTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATTACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTCAGGC  >  NC_000913/3340453‑3340723

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